Solvation and Dynamics of Sodium and Potassium in Ethylene

J. Phys. Chem. C , 2017, 121 (40), pp 21913–21920. DOI: 10.1021/acs.jpcc.7b06457. Publication Date (Web): September 18, 2017. Copyright © 2017 Amer...
1 downloads 6 Views 999KB Size
Subscriber access provided by University of Sussex Library

Article

Solvation and Dynamics of Sodium and Potassium in Ethylene Carbonate from Ab Initio Molecular Dynamics Simulations Tuan Anh Pham, Kyoung E Kweon, Amit Samanta, Vincenzo Lordi, and John E. Pask J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b06457 • Publication Date (Web): 18 Sep 2017 Downloaded from http://pubs.acs.org on September 23, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Solvation and Dynamics of Sodium and Potassium in Ethylene Carbonate from Ab Initio Molecular Dynamics Simulations Tuan Anh Pham,∗,† Kyoung E. Kweon,† Amit Samanta,‡ Vincenzo Lordi,∗,† and John E. Pask∗,‡ Quantum Simulations Group, Materials Science Division, Lawrence Livermore National Laboratory, Livermore, CA 94550, USA, and Physics Division, Lawrence Livermore National Laboratory, Livermore, CA 94550, USA E-mail: [email protected]; [email protected]; [email protected] Phone: 925-423-6501. Fax: 925-423-5733



To whom correspondence should be addressed Quantum Simulations Group, Materials Science Division, Lawrence Livermore National Laboratory, Livermore, CA 94550, USA ‡ Physics Division, Lawrence Livermore National Laboratory, Livermore, CA 94550, USA †

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract The development of sodium and potassium batteries offers a promising way to meet the scaling and cost challenges of energy storage. However, compared to Li+ , several intrinsic properties of Na+ and K+ , including their solvation and dynamics in typical organic electrolytes utilized in battery applications, are less well-understood. Here, we report a systematic investigation of Na+ and K+ in ethylene carbonate (EC) using first-principles molecular dynamics simulations. Our simulations reveal significant differences in the solvation structure and dynamical properties of Na+ and K+ compared to Li+ . We find that, in contrast to Li+ which exhibits a well-defined first solvation shell, the larger Na+ and K+ ions show more disordered and flexible solvation structures. These differences in solvation were found to significantly influence the ion dynamics, leading to larger diffusion coefficients of Na+ and K+ compared to Li+ . Our simulations also reveal a clear and interesting analog in the behavior of the ions in EC and aqueous environments, particularly in the specific ion effects on the solvent dynamics. This work provides fundamental understanding of the intrinsic properties of Na+ and K+ in organic electrolytes, which may ultimately influence the intercalation mechanism at the electrode-electrolyte interface and therefore battery performance, lifetime, and safety.

2

ACS Paragon Plus Environment

Page 2 of 26

Page 3 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Introduction Li-ion secondary batteries (LIBs) remain the most popular rechargeable batteries utilized in electronic devices due to their high energy and power densities. However, to meet the scaling and cost challenges of energy storage, the development of alternative energy storage technologies beyond traditional LIBs is necessary. 1,2 This is partly due to the limited lithium reserves in the Earth’s crust, that could lead to LIBs technological unsustainability as the demand of lithium usage in LIBs for new technologies has been steady growing. In this regard, sodium 3–7 and potassium 8–10 ion batteries have been considered as a promising way to mitigate the potential shortage of lithium due to their greater abundance and low cost. While significant efforts are ongoing to optimize host materials for both sodium 11–13 and potassium 14–16 ion batteries, several fundamental properties of Na+ and K+ in battery electrolytes remain less well-understood than those of Li+ . Of particular importance are the ion solvation structure and transport properties that play a key role in the intercalation/deintercalation mechanism and therefore battery performance. 17,18 In this regard, some aspects of the ion’s solvation have been recently addressed by both experimental and theoretical studies. Specifically, electrochemical impedance spectroscopy experiments show that Na+ insertion and desolvation processes in propylene carbonate and tetraglyme are much faster than Li+ , while yielding an activation for Na+ desolvation more than two times smaller than that of Li+ . 19 This implies that Na+ generally exhibits a weaker interaction with the solvents considered compared to Li+ . Such a conclusion is consistent with a recent combined theoryexperimental study 20 and density functional theory (DFT)-based calculations of Na+ in various organic electrolytes. 21,22 Studies of solvation properties of K+ in organic electrolytes are rather rare. The most recent DFT calculations using gas phase models indicate that compared to Li+ and Na+ , K+ displays the lowest interaction energy with the solvents considered. 23 These theoretical studies of Na+ and K+ , however, were based on static gas phase models 20–23 that neither account for temperature effects nor yield dynamical properties of the ions. On the other hand, existing molecular dynamics (MD) simulations of, e.g., Na+ 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

in organic electrolytes, have relied on classical force fields, which were derived based on specific model systems. 24 In this regard, first-principles simulations, which do not require a priori assumptions, are particularly valuable for providing a microscopic description of the solvation and dynamics of Na+ and K+ in organic electrolytes. Here, we present a detailed investigation of the solvation and transport properties of Na+ and K+ in ethylene carbonate (EC) using first-principles molecular dynamics (FPMD) simulations based on DFT. Our simulations show significant differences in the solvation structure and dynamical properties of Na+ and K+ compared to Li+ , which are directly related to the variations in their intrinsic electronic properties. In addition, our results reveal analogous behavior of these ions in EC solvent compared to aqueous environments, particularly regarding specific ion effects on solvent dynamics. The current study builds on our previous work 25,26 that examined the relationship between the structural and dynamical properties of Li+ in organic electrolytes, which is critical for optimizing the charge/discharge process in LIBs. The remainder of the paper is organized as follows. First, we describe our computational methods, including the construction of the electrolyte model. Then, we discuss our results regarding the solvation structures, dynamics, and electronic structures of Na+ and K+ in EC solvent, along with comparisons to corresponding results for Li+ . Finally, we discuss our main conclusions.

Methods The electrolytes were modeled by cubic supercells consisting of one cation (Li+ , Na+ , and K+ ), a counter anion (PF− 6 ), and 63 EC molecules. The simulation models consist of 638 total atoms in a supercell 19.283 ˚ A on a side, yielding a density of 1.32 g/cm3 and an ion concentration of 0.23 M. A low ion concentration was chosen in order to investigate the solvation structure of single ions. For the same reason, configurations were generated with

4

ACS Paragon Plus Environment

Page 4 of 26

Page 5 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

− the cation and PF− 6 dissociated, and we verified that PF6 is never associated with the cation

during the course of our simulations. 25 The FPMD simulations were carried out using the Born-Oppenheimer approximation with the VASP code, 27,28 where the electronic ground-state wave functions were optimized at each ionic step. The interatomic forces were computed using DFT with the projector augmented wave (PAW) method 29,30 and the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation for the exchange-correlation functional. 31 The electronic wave-functions were expanded in a plane-wave basis set truncated at a cutoff energy of 450 eV, and the Brillouin zone was sampled at the Γ-point. We carried out the FPMD simulations in the NVT canonical ensemble using a Nose-Hoover thermostat 32,33 with a frequency of ∼1000 cm−1 , which corresponds to a period of ∼63 fs, and a time step of 0.5 fs. Following our previous work, 25 we used a simulation temperature of 330 K to mimic an intermediate battery operating temperature and to ensure EC was not frozen (Tmelt = 310 K). For each system, equilibration runs were carried out for 7 ps, and we collected statistics over 25 ps production runs.

Results and Discussion Conformation of the First Ion Solvation Shell Our initial examination of the ion solvation structures is based on the radial distribution functions (RDF) between the solvated ions and oxygen atoms of the EC molecules, gXO (r), where X=Li+ , Na+ , or K+ . We focus on the average ion–oxygen distance, rX , that is relevant for characterizing the ion solvation. As shown in Fig. 1a, the position of the gXO (r) first A, 2.34 ˚ A, and 1.95 ˚ A maximum follows the order rK+ > rNa+ > rLi+ , yielding values of 2.80 ˚ for K+ , Na+ , and Li+ , respectively. This ordering indicates an increase in the size of the first solvation shell from Li+ to K+ and reflects the trend in the ion radius among the three ions. 34 We also report in Fig. 1b the RDFs between the ions and either carbonyl oxygen 5

ACS Paragon Plus Environment

The Journal of Physical Chemistry

10

gXO (r)

8

+

Li + Na + K

6 4 2 0

(a)

1

2

4

3

5

7

6

8

r (Å) 30 15

25 gXO (r)

+

20

gXO (r)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 26

Na + K

10

5

15 0 1

2

4

3

5 0

5

6

7

8

r (Å)

10

Li-OC Li-OE

(b) 1

2

3

4

5

6

7

8

r (Å)

Figure 1: (a) Ion–oxygen radial distribution functions, gXO (r), for solvated ions in EC, obtained from first-principles molecular dynamics simulations. Black, red, and blue lines indicate results for X=Li+ , Na+ , and K+ , respectively. (b) Radial distribution functions between Li+ and either the carbonyl oxygen atoms (OC , solid line) or the ether oxygen atoms (OE , dashed line) of EC molecules. The corresponding results for Na+ and K+ are shown in the inset. (OC ) or ether oxygen (OE ) atoms of EC molecules, which are discussed in more detail below. We then examine the distributions of the oxygen coordination number, nXO , in the first ion solvation shell to further describe the solvation geometries. Specifically, we find that nXO increases with the size of the ion solvation shell, yielding average values of 4.0, 5.7, and 7.6 for Li+ , Na+ , and K+ , respectively. Typical solvation structures of the three ions are depicted in Fig. 2. For Li+ (Fig. 2a), the solvation exhibits the well-known tetrahedral arrangement, as discussed in more detail in our previous study. 25 However, the result is significantly different for Na+ and K+ (Figs. 2b and Figs. 2c, respectively). For Na+ , we find the geometrical arrangement is intermediate between trigonal bipyramidal and square pyramidal. In contrast 6

ACS Paragon Plus Environment

Page 7 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(c)

(b)

(a)

Na+ Li+

K+

Figure 2: Typical solvation structures of (a) Li+ , (b) Na+ , and (c) K+ in EC. The solvation structures of Li+ and Na+ , respectively, can be described as tetrahedral and distorted trigonal bipyramid (or square pyramidal), while the solvation structure of K+ is less well-defined. Thick solid and transparent lines denote bonds to carbonyl and ether EC oxygen atoms, respectively. to Li+ and Na+ , the solvation structure of K+ is not as well characterized by a definite geometry, due to a large number of coordinating oxygen atoms. The representative snapshots presented in Fig. 2 for the ion solvation also hint at the notable differences between the solvation structure of Na+ and K+ compared to that of Li+ , namely, that the first solvation shell of Li+ involves only carbonyl oxygens of EC molecules, whereas those of Na+ and K+ consist of both carbonyl and ether oxygens. We quantified this feature by decomposing the RDFs for the carbonyl (OC ) and ether (OE ) oxygen atoms of the EC solvent separately, as shown in Fig. 1b. The calculated partial RDFs provide clear evidence that the first solvation shell of Li+ consists only of carbonyl oxygens, and moreover that these carbonyl oxygens rarely exchange with the rest of the liquid, as indicated by a well-defined boundary located at the first gLi+ O (r) minimum (Fig. 1a). This is in contrast to Na+ and K+ , where gXO (r) (X=Na+ or K+ ) shows a much less well-defined boundary between the first and second solvation shells, while also having notable contributions of the ether oxygens of EC molecules to the first solvation shells (inset, Fig. 1b). Together, these results suggest a much more fluid interface between the first and second solvation shells of Na+ and K+ compared to that of Li+ . We notice a general trend of increased coordination number from Li+ to Na+ to K+ , with further insight provided by examining the probability distributions of the ion–oxygen 7

ACS Paragon Plus Environment

The Journal of Physical Chemistry

100

+

Li

80

P(n)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

60 +

Na 40

+

K

20 0

2

4

6

8

10

Coordination number (n)

Figure 3: Histograms of the oxygen coordination number in the first solvation shell around Li+ (black), Na+ (red) and K+ (blue) ions in EC. The first minima in the corresponding gXO (r) were used as distance cutoffs for determination of the first solvation shells. coordination numbers, shown in Fig. 3. We see that nLi+ O exhibits a sharp and narrow distribution, showing a nearly constant value (of four) for the duration of the simulation. On the other hand, Na+ and K+ prefer an oxygen coordination number of six and eight, respectively, but also exhibit much broader distributions of nXO (K+ broader than Na+ ). These broad nXO distributions imply that the solvation shells of Na+ and K+ are rather flexible and are characterized by frequent exchange of individual oxygen atoms of the EC molecules between the first and second ion solvation shells. This conclusion is further supported by the analysis of the residence time (τ ) of oxygen atoms of EC molecules in the first ion solvation shell. We estimated τ by fitting an exponential of the form e−(t/τ )β to the time correlation function P (t) = hH(t) · H(0)i, where H(t) = 1 if a given EC oxygen atom is within the first solvation shell and 0 otherwise. 35,36 For Li+ , the residence time is well beyond the length of the simulations (∼130 ps), whereas for Na+ or K+ , τ ∼ 22-23 ps. Collectively, analyses of the time evolution of the ion first solvation shell support the interpretation of a more rigid, well-defined solvation structure of Li+ , in contrast to more flexible solvation of Na+ and K+ . This difference is directly related to the variation in the solvation energy and transport of the ions, as discussed below. We also examine the variation in the ion solvation shell due to the formation of a contact ion pair, which is particularly relevant in the high salt concentration regime. We focused on 8

ACS Paragon Plus Environment

Page 8 of 26

Page 9 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Na+ , for which the ion solvation in carbonate-based electrolytes has been more extensively studied. 20 In particular, we carried out a relatively short FPMD simulation of 5 ps, starting from an initial configuration with the cation and anion associated. This configuration was generated based on the final snapshot of a 30 ps simulation of Li+ and PF+ 6 ions being associated, 25 by replacing Li+ with Na+ . Analyses of the trajectory show that the ion pair remains associated for the entire simulation. In addition, we find that the Na+ ion first solvation shell is coordinated by ∼4.5 oxygen atoms and the anion, simply showing a reduction of ∼1 oxygen atom compared to Na+ dissociated from the counter-ion. To the best of our knowledge, the only experimental study of the Na+ solvation number in the literature was performed for linear carbonate solvents. 20 However, the ion solvation can be significantly different in linear and cyclic carbonate solvents, 25 making direct comparison with experiment difficult presently. We also compared the thermodynamic stabilities of electrolytes with associated and dissociated ion pairs by computing the average relative energies over the trajectories. Our simulations show that configurations with an associated ion pair are slightly favorable by ∼0.1 eV compared to those with a dissociated ion pair. We note that the stability of the ion pair depends on the solvent, 25 while accurate estimate of the kinetics of ion pairing requires extensive free-energy calculations, e.g., with respect to cation-anion separation; 37 however, these topics are beyond the scope of the present study. Finally, it is of interest to compare the present results with the existing literature. For A is consistent with both previous FPMD 38,39 and Li+ in EC, the value of rLi+ = 1.95 ˚ classical simulations, 36,40 which report a range between 1.96 and 2.0 ˚ A. Comparison of Na+ and K+ results with the existing literature is less straightforward due to a lack of available A simulation data for these ions in the liquid phase. Nevertheless, the calculated rNa+ = 2.34 ˚ is consistent with a previous DFT study of Na+ –EC clusters in the gas phase that reports a value of 2.36 ˚ A. 22 On the other hand, our conclusion on the solvation structure of Na+ differs significantly from a recent classical simulation that utilized the CHARMM General Force field. 24 Specifically, Ref. 24 yields a notably larger Na+ –EC carbonyl oxygen distance

9

ACS Paragon Plus Environment

The Journal of Physical Chemistry

6

∆Esol (eV)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 26

5

4 +

Li + Na + K

3 0

20

40

60

80

100

Configuration

Figure 4: Solvation energies (∆Esol ) estimated from the binding energy of the X + (EC)n clusters, where n is the number of EC molecules in the first ion solvation shell, and X ≡ Li+ (circle), Na+ (square) or K+ (triangle). The results were obtained for 100 structural configurations extracted at equal time intervals from FPMD simulations of each electrolyte. A, while showing a clear boundary between the first and second solvation shells of rNa+ = 2.5 ˚ in the gNa+ O (r). These results indicate that, unlike simulations of Li+ , larger discrepancy between classical and DFT-based simulations remains to be addressed for Na+ in EC. We also note that, to the best of our knowledge, the current study presents the first investigation of the solvation of K+ in organic cabonate-based electrolytes, which can also be utilized for benchmarking classical simulations.

Solvation Energy The ion solvation energy is directly related to the strength of the ion–EC interaction that governs the solvation structure. In addition, it plays an important role in determining the energetics of ion intercalation from the liquid electrolyte into the battery electrodes, since desolvation is required for intercalation. 17,18,41 Here, we rely on cluster calculations to estimate the ion solvation energies (∆Esol ). Specifically, ∆Esol was approximated from the binding energy of the X + (EC)n clusters, where X ≡ Li+ , Na+ or K+ , and n is the number of EC molecules in the first ion solvation shell. In particular,

∆Esol = EX + + E(EC)n − EX + (EC)n , 10

ACS Paragon Plus Environment

(1)

Page 11 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

where EX + (EC)n , E(EC)n , and EX + are the energies of the X + (EC)n and (EC)n clusters, and the ion X + in vacuum, respectively. In these calculations, configurations of the X + (EC)n clusters were extracted directly from the FPMD simulations, therefore temperature and dynamical effects are implicitly included. For each ion, we report in Fig. 4 the solvation energies computed for 100 configurations of X + (EC)n extracted at equal time intervals (250 fs) from the corresponding FPMD simulations. Focusing first on Li+ for which the solvation energy has been widely discussed in the literature, we obtained an average value of ∆Esol = 5.85 eV at the PBE level of theory, which is in good agreement with results of 5.5-6.2 eV reported in previous studies using similar gas phase calculations. 42–44 More interestingly, our result is also in reasonable agreement with a value of 5.2 eV obtained from more computationally intensive free energy calculations using MD simulations performed at the same level of theory. 45,46 This observation suggests that the ion solvation energy is largely determined by the Coulombic interactions between the ion and the EC molecules that constitute the first ion solvation shell, and thus that cluster calculations provide a good estimate for the ion solvation energy in EC. When compared to Li+ , the solvation energies of Na+ and K+ are notably smaller. We obtain values of 4.76 eV and 4.12 eV for Na+ and K+ , respectively, which are 19% and 30% smaller than the solvation energy estimated for Li+ . In addition, we find that K+ yields a smaller solvation energy than Na+ ; however, the difference (about 0.5 eV) is much less significant than the variation between the solvation energy of K+ or Na+ with respect to that of Li+ , indicating only a slightly weaker K+ solvation compared to Na+ . These results are consistent with the generally less rigid solvation structures of Na+ and K+ , compared to the strongly tetrahedral solvation of Li+ , discussed above. The results in Fig. 4 also show significantly larger fluctuations of ∆Esol for Na+ and K+ compared to Li+ . Quantitatively, the standard deviation of ∆Esol is 0.43 eV and 0.37 eV for Na+ and K+ , respectively, as compared to the corresponding value of 0.28 eV for Li+ . The larger standard deviations are associated with greater variations of the oxygen coordination

11

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 26

Lone pair

(a)

Li +

(b)

Na+

(c)

K+

Figure 5: The location of Wannier centers (yellow spheres) is shown for representative snapshots of the first solvation shell of (a) Li+ , (b) Na+ , and (c) K+ . Carbon, oxygen and hydrogen atoms are represented by blue, red and white spheres, respectively, whereas the cations are represented by transparent gray spheres. number of Na+ and K+ , which stems from the intrinsic flexibility of the ion solvation structures discussed above. Thus, both the structural analysis and solvation energy calculations collectively show that Li+ yields a rigid solvation structure, whereas those of Na+ and K+ are much more flexible, with K+ having the most flexible ion solvation of the three ions. We note that, in the context of energetics for ion intercalation, smaller desolvation energies of Na+ and K+ compared to Li+ indicate that these ions can move more easily into electrodes from the electrolyte.

Electronic Structure To further understand the origin of the intrinsic rigidity/flexibility of the different ion solvations, we explore the differences in electronic structure of each of the solvated ions. We employ maximally localized Wannier function (MLWF) analysis 47,48 as a convenient way to examine the ion–solvent interaction. With this analysis, the ground-state charge density is transformed into a local-orbital basis so that the associated positions of Wannier centers (WCs) can provide an intuitive picture of the ion polarization and bonding character. The positions of WCs around each ion, relative to the ion center, is related to the degree of localization of the ion’s electrons around the nucleus. Figure 5 shows the WCs computed for representative configurations of the ion solvation shells of Li+ , Na+ , and K+ . For Li+ , we observe that the WC of the 1s2 electrons is located at the center of the ion, in contrast to Na+ and K+ for which the four WCs associated with the 2s2 2p6 and 3s2 3p6 electrons 12

ACS Paragon Plus Environment

Page 13 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(a): Li+ MLWFs

(b): K+ MLWFs

(c): MLWFs of an Oxygen lone pair near Li +

(d): MLWFs of an Oxygen lone pair near K +

Figure 6: Isosurfaces of the maximally-localized Wannier functions (MLWFs) associated with semicore electrons of solvated (a) Li+ and (b) K+ ions in EC. Also shown are isosurfaces of the MLWFs associated with a representative oxygen lone pair of a solvating EC molecule in each case: (c) for Li+ and (d) for K+ . Positive and negative values are indicated by the yellow and blue surfaces, respectively. An identical isosurface value of 3.0 au was used in all cases. We observe in the case of Li+ that both the Li 1s2 electron and the EC oxygen lone pair electrons are drawn close to the Li+ ion, while for K+ (and Na+ , not shown) more delocalization around the ion is apparent. for Na+ and K+ , respectively, are offset from and less tightly bound to the ion center. The delocalization of the ion electrons is quantified by the ensemble-averaged scalar distance between the WCs and the ion center, dbX + −W , computed from 100 configurations extracted from each FPMD trajectory. The results reveal that the WC of Li+ is located at the ion center (dbLi+ −W = 0.00 ˚ A), while dbNa+ −W = 0.30 ˚ A and dbK+ −W = 0.35 ˚ A. The offset WCs for Na+ and K+ indicate weaker electron localization, with K+ being the weakest, while Li+ shows strong localization. This difference in the electron localization is also reflected in the shape and extent of the MLWFs for these electrons, as shown representatively in Fig. 6a and Fig. 6b for Li+ and K+ , respectively. Specifically, we observe that the MLWF of Li+ is much more tightly bound than those of K+ , consistent with the computed dbX + −W values. We also observe that the solvating oxygen orbitals protrude closer to the Li+ center than to K+ (or 13

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Na+ , not shown), as discussed further below. The above analyses reflect the Coulombic attraction between the ion nucleus and its electrons, which is further related to the interaction strength between the ion and the valence electrons of the solvating EC molecules. This interaction strength was examined by considering the MLWFs associated with the carbonyl oxygen lone pair electrons of the solvating EC molecules of Li+ and K+ , as shown in Fig. 6c and Fig. 6d, respectively. These orbitals were analyzed due to the preferred interaction of the solvated ions with the EC oxygen lone pairs that are oriented toward the ions. We find that the WLWFs of the EC oxygen lone pair show a clear tendency to share electrons with Li+ , in contrast to K+ (and Na+ ). Accordingly, our analysis implies that the origin of the weaker solvation of Na+ and K+ by EC compared to Li+ can be traced to the delocalization of Na+ and K+ electrons that results in weaker interactions between the ions and EC oxygen lone pair electrons. This ordering in the delocalization is also consistent with the electron binding energies of Li, Na, and K. For example, the first ionization potentials of the species follow the ordering: Li (5.39 eV) > Na (5.14 eV) > K (4.34 eV). 49 We note that the strength of the ion solvation is determined not only by the intrinsic properties of the ion, but also the solvent types, as discussed in our previous study for Li+ . 25

Ion Dynamics The strength of the ion solvation shell discussed above is directly related to the transport properties of the ion. For example, our previous studies have shown that the strongly solvated Li+ diffuses much more slowly than the weakly solvated PF− 6 ion in carbonate solvents, including EC; 25,26 in particular, the dynamics of the former is limited by the caging effect, in contrast to the simpler Brownian-like motion of the latter. A similar observation was made in Ref. 38 when comparing propylene-carbonate (PC) and EC solvents, where it was shown that the weaker solvation of Li+ in PC results in a relatively faster ion diffusion. In addition, examination of Li+ in different solvents, including EC, ethyl methyl carbonate (EMC), and a 14

ACS Paragon Plus Environment

Page 14 of 26

Page 15 of 26

12 10 +

Li + Na + K

8

2

MSD (Å )

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

6 4 2 0 0

2

4

6

8

10

Time (ps)

Figure 7: Mean-square displacement (MSD) versus time computed for Li+ (black), Na+ (red) and K+ (blue) ions in EC. mixture of EC and EMC, revealed that the diffusivity of Li+ varies appreciably depending on the solvent interaction strength, i.e., the ion diffusion coefficient is smaller in solvents where Li+ is more tightly solvated. 25 Similarly, it has been demonstrated that mixing EC with solvents having lower dielectric constants, such as dimethyl carbonate (DMC), weakens the Li+ solvation shell and leads to faster diffusion compared to the pure solvent. 50 Collectively, these studies suggest that the transport of Li+ in EC (and other organic solvents) may be significantly different from that of Na+ and K+ , which exhibit significantly weaker solvation shells. The transport properties of Li+ , Na+ and K+ in EC were investigated by comparing their mean square displacements (MSD) versus time. As shown in Fig. 7, the calculated MSDs clearly indicate that Li+ diffuses much more slowly than Na+ and K+ . Specifically, the diffusion coefficients were estimated from the slopes of the MSD plots as 4.3 ± 0.7, 12.1 ± 2.0, and 11.2 ± 0.8 cm2 /µs for Li+ , Na+ and K+ , respectively. As a result, our simulations indicate that the diffusion coefficient of Li+ is about 3 times smaller than those of either Na+ or K+ , which can be directly related to a significantly stronger solvation shell of Li+ . On the other hand, although the K+ solvation is (slightly) weaker than Na+ , the computed diffusion coefficients are not significantly different, implying that the Na+ and K+ transport are more dictated by the bulk solvent properties, due to the lack of the strong ion–solvent

15

ACS Paragon Plus Environment

The Journal of Physical Chemistry

12 10

+

Li + Na + K

8

2

MSD (Å )

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 26

6 4 2 0 0

2

4

6

8

10

Time (ps)

Figure 8: Mean-square displacement (MSD) computed for all EC molecules in the electrolyte with Li+ (black), Na+ (red) and K+ (blue) ions. The MSD computed for EC molecules only in the first solvation shell of Li+ is presented by the dashed line. interaction. (We also note the close quantitative similarity of the electron delocalization for Na+ and K+ discussed in the previous section.) The relatively large uncertainties in the estimated diffusion coefficients also implies that significantly more statistics are required for a better quantitative description of the ion diffusion coefficients; however, the distinction between Li+ and both Na+ and K+ is clear. Reducing these errors significantly further to allow more precise distinction between the diffusion coefficients of Na+ and K+ in FPMD simulations is nontrivial, since it would require many long (at least several tens of ps or more) and independent simulations. The ion dynamics is clearly governed by the nature of the ion–solvent interaction; conversely, the dynamics of EC molecules may also be affected by the presence of the ion. The effect of the ions on the solvent molecules is analyzed with MSDs computed for the EC molecules in the presence of the ions, as shown in Fig. 8. The results show a slight suppression of the EC diffusion in the electrolyte containing Li+ compared to Na+ and K+ . In particular, we obtained values of 14.4 ± 0.6, 16.0 ± 0.5, and 15.9 ± 0.2 cm2 /µs for the diffusion coefficients of EC molecules in electrolytes containing Li+ , Na+ and K+ , respectively. To understand the origin of these specific ion effects, we investigated the dynamics of EC molecules only in the first ion solvation shell, which are directly influenced by the ion, and we focus on Li+ that exhibits the strongest ion solvation. The result is shown by 16

ACS Paragon Plus Environment

Page 17 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

the dashed line in Fig. 8, indicating a much slower diffusion of EC molecules in the first solvation shell compared to bulk EC. The computed diffusion coefficient of the first solvation shell EC molecules is 4.9 ± 0.4 cm2 /µs, which is about three times smaller than the overall diffusion coefficients of bulk EC molecules. These results provide clear evidence that the low diffusivity of the EC molecules near Li+ is responsible for the slightly lower diffusivity of EC molecules in the electrolyte with Li+ compared to those with Na+ or K+ , and further highlights the specific ion effects in the organic electrolytes. We note that previous work 26,51 describes a diffusion mechanism for Li+ in EC that includes the four solvent molecules as a complex, which is reflected in the results presented here. The weaker, less-well-structured solvation of Na+ and K+ thus contributes to significantly higher diffusivity of those ions. Finally, we relate and compare the present findings with existing studies of these ions in aqueous solutions. For the latter, it is known that the dynamics of water molecules in salt solutions deviates from that of pure liquid water, and the concepts of structure-breaking and structure-making described by the empirical Jones-Dole viscosity coefficient, B, have been used for a qualitative explanation of the ion effects on the water dynamics. 52 Specifically, ions with B > 0 reflect a strong hydration and are classified as “structure-making”, while B < 0 denotes “structure-breaking.” In this context, Li+ is classified as a strong structure-maker (B=0.14), Na+ is a weak structure-maker (B=0.06), and K+ is a weak structure-breaker (B = −0.01) in water. Accordingly, the ordering of the solvation strength is the same for ions in water and EC, i.e., Li+ > Na+ > K+ . In addition, it has been shown experimentally that the water diffusion coefficient in aqueous solutions decreases with increasing LiCl or NaCl concentrations, with the strongest reduction observed for LiCl, whereas the variation is less significant for the KCl solution. 53,54 This implies that the diffusion of water molecules is slowest for solutions containing Li+ , indicating analagous ion dynamical behavior in the aqueous solution and EC. Finally, it is also worth pointing out that the oxygen coordination number in the ion first solvation shell is quite similar for water and EC, particularly for Li+ and Na+ , which yield values of 4.0 and 5.3 in water, respectively, 55,56 while K+ ex-

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

hibits a slightly smaller oxygen coordination number (6.0-6.1) in water compared to EC (see Fig. 3). 56,57

Conclusions To conclude, we investigate the solvation and transport properties of Na+ and K+ in the EC solvent, compared to Li+ , using first-principles molecular dynamics simulations based on density functional theory. Our simulations reveal significant differences in the solvation structure and dynamical properties of Na+ and K+ compared to Li+ . We find that, in contrast to Li+ , which exhibits a well-defined first solvation shell with a tetrahedral arrangement of carbonyl oxygen atoms of EC molecules, the larger Na+ and K+ ions show more disordered and flexible solvation structures. Specifically, Na+ and K+ exhibit increased average coordination numbers of 5.7 and 7.6, respectively, which are mostly comprised of carbonyl oxygens but also involve a non-negligible fraction of ether oxygens. The less-ordered solvation structures for the larger ions is reflected by weaker solvation energies; specifically, we obtain values of 5.85, 4.72, and 4.12 eV for Li+ , Na+ , and K+ , respectively. Using the analysis based on maximally localized Wannier functions, we trace the origin of this behavior to the increased delocalization of Na+ and K+ electrons that results in weaker interactions between the ions and EC oxygen lone-pair electrons. Our simulations imply a complex interplay between the ion dynamics and ion solvation strength in the EC solvent. Specifically, we find that weaker-solvated Na+ and K+ ions show a much faster diffusion than Li+ , yielding a diffusion coefficient about three times higher than that of Li+ . On the other hand, although Na+ exhibits a slightly stronger solvation structure than K+ , the dynamics of Na+ and K+ are not significantly different, indicating that the inhibited diffusivity of Li+ stems from its strong interaction with the solvent EC molecules, whereas the Na+ and K+ transport is more likely indicative of the bulk solvent due to the lack of strong ion–solvent interactions. Our study also reveals specific ion effects

18

ACS Paragon Plus Environment

Page 18 of 26

Page 19 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

on the dynamics of the solvent; particularly, we demonstrate that EC molecules diffuse slower in the Li+ first solvation shell, which is in turn dictated by the strong ion–solvent interaction. Furthermore, the strong ion–solvent interactions of Li+ reflects the “structure-maker” nature of Li+ in the EC solvent, indicating an analog to the behavior of solvated ions in aqueous electrolytes. Finally, our results indicate that the use of K+ or Na+ instead of Li+ in an organic electrolyte may benefit the charge/discharge rates in energy storage devices, due to their fast transport. In the context of energetics for ion intercalation, the smaller solvation energy of K+ and Na+ compared to that of Li+ implies easier migration to the electrode from the electrolyte. We note that ion solvation and dynamics, and the energy barrier for ion intercalation/de-intercalation, also depend on the solvent types as well as chemistry of the electrode/electrolyte interface, which are topics of future study. Work is also in progress to compute the kinetics of ion pairings, and to understand the dynamics of the ion pairs in different solvents.

Acknowledgement This work was performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344. Support for this work was provided through Scientific Discovery through Advanced Computing (SciDAC) program funded by the U.S. Department of Energy, Office of Science, Advanced Scientific Computing Research and Basic Energy Sciences. T.A.P acknowledges the support from the Lawrence Fellowship Program. Computational resources were from the Lawrence Livermore National Laboratory Institutional Computing Grand Challenge Program.

References (1) Palomares, V.; Serras, P.; Villaluenga, I.; Hueso, K. B.; Carretero-Gonz´alez, J.; Rojo, T. Na-ion batteries, recent advances and present challenges to become low cost energy 19

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

storage systems. Energy Environ. Sci. 2012, 5, 5884–5901. (2) Palomares, V.; Casas-Cabanas, M.; Castillo-Mart´ınez, E.; Han, M. H.; Rojo, T. Update on Na-based battery materials. A growing research path. Energy Environ. Sci. 2013, 6, 2312–2337. (3) Ellis, B. L.; Nazar, L. F. Sodium and sodium-ion energy storage batteries. Curr. Opin. Solid State Mater. Sci. 2012, 16, 168–177. (4) Slater, M. D.; Kim, D.; Lee, E.; Johnson, C. S. Sodium-ion batteries. Adv. Funct. Mater 2013, 23, 947–958. (5) Yabuuchi, N.; Kubota, K.; Dahbi, M.; Komaba, S. Research development on sodium-ion batteries. Chem. Rev. 2014, 114, 11636–11682. (6) Sawicki, M.; Shaw, L. L. Advances and challenges of sodium ion batteries as post lithium ion batteries. RSC Adv. 2015, 5, 53129–53154. (7) Hwang, J.-Y.; Myung, S.-T.; Sun, Y.-K. Sodium-ion batteries: present and future. Chem. Soc. Rev. 2017, 46, 3529–3614. (8) Eftekhari, A. Potassium secondary cell based on Prussian blue cathode. J. Power Sources 2004, 126, 221–228. (9) Eftekhari, A.; Jian, Z.; Ji, X. Potassium secondary batteries. ACS Appl. Mater. Interfaces 2017, 9, 4404–4419. (10) Wang, Y.; Chen, R.; Chen, T.; Lv, H.; Zhu, G.; Ma, L.; Wang, C.; Jin, Z.; Liu, J. Emerging non-lithium ion batteries. Energy Storage Mater. 2016, 4, 103–129. (11) Kim, S.-W.; Seo, D.-H.; Ma, X.; Ceder, G.; Kang, K. Electrode materials for rechargeable sodium-ion batteries: potential alternatives to current lithium-ion batteries. Adv. Energy Mater 2012, 2, 710–721. 20

ACS Paragon Plus Environment

Page 20 of 26

Page 21 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(12) Han, M. H.; Gonzalo, E.; Singh, G.; Rojo, T. A comprehensive review of sodium layered oxides: powerful cathodes for Na-ion batteries. Energy Environ. Sci. 2015, 8, 81–102. (13) Wang, L. P.; Yu, L.; Wang, X.; Srinivasan, M.; Xu, Z. J. Recent developments in electrode materials for sodium-ion batteries. J. Mat. Chem. A 2015, 3, 9353–9378. (14) Padigi, P.; Thiebes, J.; Swan, M.; Goncher, G.; Evans, D.; Solanki, R. Prussian green: a high rate capacity cathode for potassium ion batteries. Electrochim. Acta 2015, 166, 32–39. (15) Jian, Z.; Luo, W.; Ji, X. Carbon electrodes for K-ion batteries. J. Am. Chem. Soc 2015, 137, 11566–11569. (16) Zhang, C.; Xu, Y.; Zhou, M.; Liang, L.; Dong, H.; Wu, M.; Yang, Y.; Lei, Y. Potassium Prussian blue nanoparticles: A low-cost cathode material for potassium-ion batteries. Adv. Funct. Mater 2016, 27, 1604307. (17) Xu, K.; von Cresce, A.; Lee, U. Differentiating contributions to “ion transfer” barrier from interphasial resistance and Li+ desolvation at electrolyte/graphite interface. Langmuir 2010, 26, 11538–11543. (18) Xu, K.; von Cresce, A. Interfacing electrolytes with electrodes in Li ion batteries. J. Mater. Chem. 2011, 21, 9849–9864. (19) Browning, K. L.; Sacci, R. L.; Veith, G. M. Energetics of Na+ transport through the electrode/cathode interface in single solvent electrolytes. J. Electrochem. Soc 2017, 164, A580–A586. (20) Cresce, A. V.; Russell, S. M.; Borodin, O.; Allen, J. A.; Schroeder, M. A.; Dai, M.; Peng, J.; Gobet, M. P.; Greenbaum, S. G.; Rogers, R. E. et al. Solvation behavior of carbonate-based electrolytes in sodium ion batteries. Phys. Chem. Chem. Phys. 2017, 19, 574–586. 21

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(21) Okoshi, M.; Yamada, Y.; Yamada, A.; Nakai, H. Theoretical analysis on de-solvation of lithium, sodium, and magnesium cations to organic electrolyte solvents. J. Electrochem. Soc 2013, 160, A2160–A2165. (22) Shakourian-Fard, M.; Kamath, G.; Smith, K.; Xiong, H.; Sankaranarayanan, S. K. Trends in Na-ion solvation with alkyl-carbonate electrolytes for sodium-ion batteries: insights from first-principles calculations. J. Phys. Chem. C 2015, 119, 22747–22759. (23) Okoshi, M.; Yamada, Y.; Komaba, S.; Yamada, A.; Nakai, H. Theoretical analysis of interactions between potassium ions and organic electrolyte solvents: A comparison with lithium, sodium, and magnesium ions. J. Electrochem. Soc 2017, 164, A54–A60. (24) Kamath, G.; Cutler, R. W.; Deshmukh, S. A.; Shakourian-Fard, M.; Parrish, R.; Huether, J.; Butt, D. P.; Xiong, H.; Sankaranarayanan, S. K. In silico based rankorder determination and experiments on nonaqueous electrolytes for sodium ion battery applications. J. Phys. Chem. C 2014, 118, 13406–13416. (25) Ong, M. T.; Verners, O.; Draeger, E. W.; van Duin, A. C.; Lordi, V.; Pask, J. E. Lithium ion solvation and diffusion in bulk organic electrolytes from first-principles and classical reactive molecular dynamics. J. Phys. Chem. B 2015, 119, 1535–1545. (26) Ong, M. T.; Bhatia, H.; Gyulassy, A. G.; Draeger, E. W.; Pascucci, V.; Bremer, P.T.; Lordi, V.; Pask, J. E. Complex ion dynamics in carbonate lithium-ion battery electrolytes. J. Phys. Chem. C 2017, 121, 6589–6595. (27) Kresse, G.; Furthm¨ uller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169–11186. (28) Kresse, G.; Furthm¨ uller, 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. (29) Bl¨ochl, P. E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953. 22

ACS Paragon Plus Environment

Page 22 of 26

Page 23 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(30) Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmentedwave method. Phys. Rev. B 1999, 59, 1758. (31) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865. (32) Nos´e, S. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 1984, 81, 511–519. (33) Hoover, W. G. Canonical dynamics: equilibrium phase-space distributions. Phys. Rev. A 1985, 31, 1695. (34) Nightingale Jr, E. Phenomenological theory of ion solvation. Effective radii of hydrated ions. J. Phys. Chem. 1959, 63, 1381–1387. (35) Borodin, O.; Smith, G. D. LiTFSI structure and transport in ethylene carbonate from molecular dynamics simulations. J. Phys. Chem. B 2006, 110, 4971–4977. (36) Borodin, O.; Smith, G. D. Quantum chemistry and molecular dynamics simulation study of dimethyl carbonate: ethylene carbonate electrolytes doped with LiPF6 . J. Phys. Chem. B 2009, 113, 1763–1776. (37) Dang, L. X.; Chang, T.-M. Rate theory of solvent exchange and kinetics of Li+ − BF− 4 /PF6 ion pairs in acetonitrile. J. Chem. Phys. 2016, 145, 094502.

(38) Ganesh, P.; Jiang, D.-e.; Kent, P. Accurate static and dynamic properties of liquid electrolytes for Li-ion batteries from ab initio molecular dynamics. J. Phys. Chem. B 2011, 115, 3085–3090. (39) Bhatt, M. D.; Cho, M.; Cho, K. Density functional theory calculations and ab initio molecular dynamics simulations for diffusion of Li+ within liquid ethylene carbonate. Modell. Simul. Mater. Sci. Eng. 2012, 20, 065004.

23

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(40) Li, T.; Balbuena, P. B. Theoretical studies of lithium perchlorate in ethylene carbonate, propylene carbonate, and their mixtures. J. Electrochem. Soc 1999, 146, 3613–3622. (41) Xu, K.; von Wald Cresce, A. Li+ -solvation/desolvation dictates interphasial processes on graphitic anode in Li ion cells. J. Mater. Res. 2012, 27, 2327–2341. (42) Bhatt, M. D.; Cho, M.; Cho, K. Conduction of Li+ cations in ethylene carbonate (EC) and propylene carbonate (PC): comparative studies using density functional theory. J. Solid State Electrochem. 2012, 16, 435–441. (43) Skarmoutsos, I.; Ponnuchamy, V.; Vetere, V.; Mossa, S. Li+ solvation in pure, binary, and ternary mixtures of organic carbonate electrolytes. J. Phys. Chem. C 2015, 119, 4502–4515. (44) Cui, W.; Lansac, Y.; Lee, H.; Hong, S.-T.; Jang, Y. H. Lithium ion solvation by ethylene carbonates in lithium-ion battery electrolytes, revisited by density functional theory with the hybrid solvation model and free energy correction in solution. Phys. Chem. Chem. Phys. 2016, 18, 23607–23612. (45) Leung, K.; Tenney, C. M. Toward first principles prediction of voltage dependences of electrolyte/electrolyte interfacial processes in lithium ion batteries. J. Phys. Chem. C 2013, 117, 24224–24235. (46) Li, Y.; Leung, K.; Qi, Y. Computational exploration of the Li-electrode— electrolyte interface in the presence of a nanometer thick solid-electrolyte interphase layer. Acc. Chem. Res. 2016, 49, 2363–2370. (47) Mostofi, A. A.; Yates, J. R.; Lee, Y.-S.; Souza, I.; Vanderbilt, D.; Marzari, N. wannier90: A tool for obtaining maximally-localised Wannier functions. Comput. Phys. Commun. 2008, 178, 685–699.

24

ACS Paragon Plus Environment

Page 24 of 26

Page 25 of 26

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(48) Marzari, N.; Mostofi, A. A.; Yates, J. R.; Souza, I.; Vanderbilt, D. Maximally localized Wannier functions: Theory and applications. Rev. Mod. Phys. 2012, 84, 1419. (49) Kramida, A.; Yu. Ralchenko,; Reader, J.; and NIST ASD Team, NIST Atomic Spectra Database (version 5.3). National Institute of Standards and Technology, Gaithersburg, MD, 2015, http://physics.nist.gov/asd, (accessed August 28, 2017). (50) Borodin, O.; Olguin, M.; Ganesh, P.; Kent, P. R.; Allen, J. L.; Henderson, W. A. Competitive lithium solvation of linear and cyclic carbonates from quantum chemistry. Phys. Chem. Chem. Phys. 2016, 18, 164–175. (51) Xu, K. Nonaqueous liquid electrolytes for lithium-based rechargeable batteries. Chem. Rev. 2004, 104, 4303–4418. (52) Marcus, Y. Effect of ions on the structure of water: structure making and breaking. Chem. Rev. 2009, 109, 1346–1370. (53) M¨ uller, K.; Hertz, H. A parameter as an indicator for water-water association in solutions of strong electrolytes. J. Chem. Phys. 1996, 100, 1256–1265. (54) Ben Ishai, P.; Mamontov, E.; Nickels, J. D.; Sokolov, A. P. Influence of Ions on Water Diffusion–A Neutron Scattering Study. J. Phys. Chem. B 2013, 117, 7724–7728. (55) Ikeda, T.; Boero, M.; Terakura, K. Hydration of alkali ions from first principles molecular dynamics revisited. J. Chem. Phys. 2007, 126, 01B611. (56) Pham, T. A.; Ogitsu, T.; Lau, E. Y.; Schwegler, E. Structure and dynamics of aqueous solutions from PBE-based first-principles molecular dynamics simulations. J. Chem. Phys. 2016, 145, 154501. (57) Pham, T. A.; Mortuza, S. G.; Wood, B. C.; Lau, E. Y.; Ogitsu, T.; Buchsbaum, S. F.; Siwy, Z. S.; Fornasiero, F.; Schwegler, E. Salt solutions in carbon nanotubes: The role of cation- π interactions. J. Phys. Chem. C 2016, 120, 7332–7338. 25

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Graphical TOC Entry

Na+ Li+

K+

Increasing Ion Coordination Number

26

ACS Paragon Plus Environment

Page 26 of 26