Spectroscopic and Density Functional Theory ... - ACS Publications

Jan 3, 2017 - 11 software. To ensure confidence in the deconvolution, constraints of a consistent Gaussian to Lorentzian ratio were maintained on all ...
6 downloads 29 Views 3MB Size
Subscriber access provided by GAZI UNIV

Article

Spectroscopic and Density Functional Theory Characterization of Common Lithium Salt Solvates in Carbonate Electrolytes for Lithium Batteries Navid Chapman, Oleg Borodin, Taeho Yoon, Cao Cuong Nguyen, and Brett L Lucht J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b12234 • Publication Date (Web): 03 Jan 2017 Downloaded from http://pubs.acs.org on January 4, 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 48

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

Spectroscopic and Density Functional Theory Characterization of Common Lithium Salt Solvates in Carbonate Electrolytes for Lithium Batteries Navid Chapman,1 Oleg Borodin,2 Taeho Yoon,1 Cao Cuong Nguyen,1 and Brett L. Lucht*,1

1

Department of Chemistry, University of Rhode Island, Kingston, Rhode Island 02881, USA

2.

Electrochemistry Branch, Sensor and Electron Devices Directorate, Power and Energy Division, U.S. Army Research Laboratory Adelphi, Maryland, 20783, USA e-mail: Brett Lucht: [email protected], Oleg Borodin: [email protected]

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

Page 2 of 48

Abstract The structure and composition of lithium ion solvation spheres of electrolyte solutions composed of common lithium salts (LiTFSI, LiPF6, LiBF4, and LiClO4) dissolved in aprotic polar linear and cyclic carbonate solvents (propylene carbonate (PC) or dimethyl carbonate (DMC)) has been investigated via a combination of FTIR,

13

C NMR spectroscopy and density

functional theory (DFT). Results from the two different spectroscopic methods are in strong agreement with each other and with predictions from quantum chemistry calculations. The coordination of the carbonyl oxygen of the solvents to the lithium cation is observed by IR spectroscopy. The ratio of coordinated to uncoordinated PC and DMC has been used to determine solvent coordination numbers which range from 2 to 5 depending on salt, solvent, and concentration. The relative stability of the lithium-anion solvates were examined using DFT employing the cluster – continuum approach including changes to the intensity and frequency of the IR bands along with the populations of the cis-cis and cis-trans conformers of DMC in the lithium ion solvation shell. Solvent coordination is dependent upon the nature of the salt. Weakly associating salts, LiTFSI, LiPF6, and LiClO4, dissociate to a similar degree with LiPF6 being the most dissociated, while LiBF4 had significantly less dissociation in both solvents. This investigation provides significant insight into the solution structure commonly used LIB electrolytes over wide range of salt concentrations.

2 ACS Paragon Plus Environment

Page 3 of 48

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 There is significant need to further improve lithium Ion Batteries (LIBs) due to the increased global demand for the portable energy storage devices.1-3 While the electrolyte is a critical component of the LIB, our current understanding of the solution structures of these electrolytes is incomplete. Further investigation is important since the electrolyte solution structure plays a critical role in the performance of electrochemical systems.4 These fundamental physicochemical properties affect many aspects of performance such as electrochemical stability window, working temperature range, and ionic conductivity.5-11 Determining the solution structure and relative binding affinities of different solvents is also of interest since it would provide insight into electrolyte reduction and solid electrolyte interface (SEI) formation on the on the anode surface.12 The lithium ion solvation structure and aggregation state often determine the electrolyte components that undergo reduction and to generate the SEI.13-16 Salt aggregation can lead to a significant change in the reduction potential of the electrolyte.13, 17-18 Binary mixtures of cyclic and linear carbonates were first implemented by Tarascon and Guyomard to optimize the properties of the electrolyte since the cyclic carbonates have high permittivity and high viscosity while the linear carbonates have low viscosity and low permittivity.19 However, investigations of the solution structures of cyclic and linear carbonates with different lithium salts are limited.11,

20-27

Using LiPF6 as a salt, Seo et al. previously

examined cation-solvent coordination via IR and NMR spectroscopies.11 This study expands on Seo’s work by examining solvation and ion pairing with several different salts. Spectroscopic methods including IR, Raman, and NMR can be used to indirectly observe solvent coordination by detecting changes in the local electronic environment of the ions or solvent molecules.

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

Page 4 of 48

In this report, IR and NMR spectroscopy have been used to investigate the cation solvent interactions resulting from lithium salts dissolved in carbonate solvents. Four lithium salts (lithium bis(trifluoromethane)sulfonamide (LiTFSI), LiPF6, LiBF4, LiClO4) and two carbonate solvents (Propylene carbonate (PC) and dimethyl carbonate (DMC)) were investigated over a wide concentration range from a solvent to Li ratio of 20:1 to 4:1.28 Density functional theory (DFT) and the more computationally expensive wave function-based quantum chemistry (QC) methods were also used to examine the relative stability of the lithium solvates and changes in the IR spectra due to complexation of the carbonyl group of the solvent to the Li+ cation. Populations of two most stable DMC conformers in lithium solvates, cis-cis and cis-trans, have also been examined. In this study, previous QC27,

29-36

and simulation20,

24, 29, 37-47

studies of Li+

solvation are built up and utilize the cluster-continuum approach which explicitly considers solvent and anion in the first coordination shell of Li+, while solvent molecules outside of the first solvation shell are described using an implicit solvation model.21, 48

Experimental Section PC was used as-received (MP Biomedicals). DMC was used as-received (Sigma-Aldrich, anhydrous). LiTFSI, LiPF6, LiBF4, and LiClO4 were also used as-received (Sigma-Aldrich). The four lithium salts (LiTFSI, LiPF6, LiBF4, LiClO4) were chosen due to the interest in these salts in lithium ion battery electrolytes. The salts readily dissolve and dissociate in non-aqueous polar aprotic electrolyte solvents, have high ionic mobility, are non-reactive with electrolyte solvents, and have relatively good thermal stability.3, 49-51 The salt-to-solvent ratios (1:20 to 1:4) have been selected to incorporate a full range of salt concentrations with the conventional concentration of ~1 M in the center of the range. Solutions were prepared in an air-free environment and analyzed

4 ACS Paragon Plus Environment

Page 5 of 48

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

shortly after preparation at room temperature to minimize contamination by water.3, 52 Constant temperature was maintained.4, 20 Samples were prepared in an inert atmosphere (N2) glovebox ( LiTFSI > LiClO4 > LiBF4. This trend is in good agreement with the experimentally observed trends for the Li+ solvation numbers. Importantly, the cluster-continuum DFT calculations not only predict the correct trends in the lithium-anion CIP formation, but also yield a reasonable estimate of the fraction of SSIP in dilute solution. Specifically, application of the Boltzmann factors (eq. 3) suggests that in PC LiBF4 is quite aggregated with only 55% SSIP, while LiClO4 has around 10% CIP. Both LiPF6 and LiTFSI are largely dissociated. DFT calculations also show that both mono-dentate and bidentate LiTFSI binding will be present with the monodentate binding dominating. In contrast to the CIP free energies with the PCM model, application of the SMD model predicts higher 22 ACS Paragon Plus Environment

Page 23 of 48

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

energies for CIP formation vs. SSIP. Thus, it will yield more dissociated electrolytes and higher dissociation for LiBF4 than LiClO4. When energies are corrected by subtracting the PBE/631+(d,p)-G4MP2 energy differences calculated using SMD implicit solvent for ion pairs as shown in Table 2, the binding energies for the SMD solvation model start to follow the experimental tends. Another important observation is that when a PC was explicitly included in the first solvation shell of Li+ (Table 4), the difference for the CIP formation reaction from SSIP is much smaller between the SMD and PCM solvation models compared to the differences observed for the CIP formation reaction when all solvent was treated implicitly. Relative stability of the Li+(PC)5 vs Li+(PC)4 solvates were estimated from the PBE/631+G(d,p), PCM(PC) cluster-continuum calculations. Free energy of the Li+(PC)5 solvate was found to be 9 kcal/mol higher than the free energy of the Li+(PC)4 solvate indicating that five PC coordinated Li+ cation SSIP are not expected in the PC-based electrolyte solutions. Interestingly, the energy of the Li+(PC)5 solvate was found to be 0.6 kcal/mol lower than Li+(PC)4 solvate showing that entropic penalty destabilizes higher solvation number complexes in accord with the recent study of the Li+(EC)n-based complexes.73 Thus, DFT calculation support the Li+ solvation numbers from Raman measurements62 at 20:1 EC:Li around four for the strongly dissociating salts (LiTFSI, LiPF6) instead of the higher solvation number of five for the PC-based electrolyte obtained from IR spectra for PC:Li=20 for the same highly dissociating salts.

Table 4. The contact ion pair formation reaction free energy at 298.15K (∆G) for the reaction Li(PC)4+Anion -→

LiAnion(PC)3+PC from PBE/6-31+G(d,p) calculations using PCM and SMD

solvation model with PC parameters. Solvation model: energy

PCM

SMD

SMD

PCM(corr.)

SMD(corr.)a

none, gas-phase

Solvation model:

PCM

PCM

SMD

PCM

SMD

none, gas-phase

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

Page 24 of 48

geometry optimization LiTFSI(mono)

2.2

1.4

LiTFSI(bidentate)+3PC

2.8

5.3

LiPF6

4.1

5.4

5.6

LiClO4

1.4

2.3

LiBF4

0.1

3.4

a

3.4

2.6

3.8

-52.9

2.6

4.1

-54.4

2.1

2.1

2.8

-59.5

2.8

0.0

2.7

-62.4

3.2

free energies were corrected using the G4MP2-PBE/6-31+G(d,p) energy differences from Table 2 for

ion pairs in SMD implicit solvent.

Lithium Solvates in Dimethyl Carbonate. Stability of CIP and SSIP Li+ solvates in DMC was

examined

by

calculating

free

energies

for

the

Li(DMCcc)4+Anion

->

LiAnion(PC)3+DMCcc reaction using PBE/6-31+G(d,p) DFT with SMD and PCM solvation models as shown in Table 5 and correspond to geometries shown in Figure 9. The build-in parameters for DiButylEther (ε=3.0473) were utilized to model DMC as implicit solvent. The CIP stability vs. SSIP obtained from PBE/6-31+G(d,p) DFT calculations was also corrected for the differences between PBE and G4MP2 binding free energies in implicit solvent shown in Table 2. We find that while absolute values for the CIP vs. SSIP stability are different for solvates calculated using SMD and PCM solvation models, the trends for the lithium salts are similar. Specifically, LiTFSI and LiPF6 have much smaller CIP formation energies than LiClO4 and LiBF4. LiTFSI monodentate and bidentate binding in CIP are also similar. After correction for the differences between G4MP2 and PBE/6-31+G(d,p), the solvate of LiBF4 has the highest CIP formation energies indicating that it is expected to have the lowest number of DMC in the first solvation shell of Li+ which agrees with IR data shown in Figure 4. Overall, the propensity for CIP formation has the following order: LiBF4 > LiClO4 > LiPF6 > LiTFSI. This agrees with the reverse trend for the Li+ solvation numbers extracted from IR measurements for dilute solutions. High negative free energies for CIP solvates vs. SSIP in DMC indicate that a 24 ACS Paragon Plus Environment

Page 25 of 48

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

negligible amount of SSIP is expected even in dilute solutions of these salts, while experiments indicate DMC solvation numbers approaching four in dilute solutions which are indicative of some SSIP formation. At least two factors might influence the solvate stabilities. First, previous work indicated that the quadrupole moment of DMC is expected to significantly contribute to the solvation and might not be adequately accounted for in the implicit solvation model.77 Second, the presence of the high dipole moment DMC(cis-trans) conformer denoted as DMCct could influence the Li+ solvate stabilities.

Figure 9. Optimized contact ion pair geometries for the Li+(DMCcc)3Anion and Li+(DMCcc)4 solvates from PBE/6-31+G(d,p) DFT calculations with PCM(ε=3) implicit solved model surrounding the solvates.

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

Page 26 of 48

Table 5. The contact ion pair formation reaction energy (∆E) 0 K and free energy (∆G) at 298.15K for the reaction Li(DMCcc)4+Anion → LiAnion(PC)3+DMCcc from PBE/6-31+G(d,p) calculations using PCM and SMD(ε=3.05) solvation model. SMD ∆G

a

PCM ∆G(corrected)a ∆G

∆G(corrected)a

LiTFSI (monodentate) -24.3

-23.9

-17.3 -16.9

LiTFSI (bidentate)

-23.3

-22.9

-17.6 -17.2

LiPF6

-23.0

-24.5

-16.8 -18.3

LiClO4

-28.7

-28.0

-21.3 -20.6

LiBF4

-28.3

-28.4

-21.5 -21.6

free energies were corrected using the G4MP2-PBE/6-31+G(d,p) energy differences from Table 2 for

ion pairs in SMD implicit solvent.

In order to further explore sensitivity of the solvation energy to the size of the Li+ solvation shell in the cluster – continuum approach, energy and free energy for adding an additional DMCcc was calculated corresponding to the reaction Li+(DMCcc)5 → Li+(DMCcc)4 + DMCcc in implicit solvent as shown in Figure S3. The reaction energy was found to be -2.9 kcal/mol, while reaction free energy is 11.9 kcal/mol. Such a high free energy difference between the Li+(DMCcc)5 and Li+(DMCcc)4 description indicates a high sensitivity of the solvation free energy on the Li+ shell description. Of course, BOMD or force field-based MD simulations with an explicit representation of all solvent would alleviate this drawback but they are either significantly more computationally expensive or rely on the accuracy and adequacy of the force field used in simulations. We have also explored an alternative way to estimate the stability of CIP vs. SSIP solvates in DMC by directly comparing their stability from LiAnion(DMCcc)4 calculations as shown in Figure 10. We stress that while we explored a number of initial configurations for each

26 ACS Paragon Plus Environment

Page 27 of 48

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

CIP and SSIP solvate and utilized only the most energetically stable configurations for our analysis, the structures are not guaranteed to be global minima. They were confirmed to be local minima by frequency calculations with the exception of the SSIP for LiPF6(DMC)4 shown in Figure 10a which has one imaginary frequency. The later structure has the Li-P distance constrained, otherwise it would converge to the CIP geometry shown in Figure 10b as a result of geometry optimization. The calculated free energies of CIP vs. SSIP are less negative from this analysis than the CIP formation energies listed in Table 5. Note that in the SSIPs in Figure 10, the Li+ and Anion- interact through the first solvation shell of the low dielectric constant DMC solvent, while calculations summarized in Table 5 assume that the Li+(DMC)4 complex is at infinite separation from the anion, which is solvated only via PCM(ε=3) model.

27 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 28 of 48

Figure 10. Optimized geometries, relative energies (∆E) and free energies (∆G) of CIP vs. SSIP for the LiAnion(DMC)4 solvates from PBE/6-31+G(d,p) calculations with PCM(ε=3).

Furthermore, our cluster-continuum calculations assume that the DMC solvent stays in the cis-cis conformation and that the dielectric constant of DMC liquid ε≈3.0 is adequate. 28 ACS Paragon Plus Environment

Page 29 of 48

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

Previous MD simulations and experiments suggested that population of the high dipole moment DMC cis-trans conformer increases when DMC is bond to a Li+ cation.64 The appearance of this high dipole moment conformer near Li+ would suggest that the local dielectric constant is expected to be significantly higher than bulk DMC dielectric constant. In order to further characterize populations of DMCcc and DMCct in the first solvation shell of Li+, energies of the LiClO4(DMCcc)m(DMCct)n, n+m=3 solvates were calculated using DFT as shown in Figure 11. Because most salts form CIP in DMC, these solvates are representative complexes for the investigation of the DMCct vs. DMCcc population in the first solvation shell of Li+. The dispersion correction was recently shown to influence the relative stability of the Li+ solvates with high and low dipole moment carbonates,

78

therefore, additional calculations were

performed using PBE+GD3BJ D3 version of Grimme’s dispersion with Becke-Johnson damping. DFT calculations shown in Figure 11 demonstrate that the most energetically favorable solvates have a mixture of DMCcc and DMCct in the first solvation shell of Li+ because the free energies of all solvates as shown in Figure 11 differ by less than 4 kcal/mol. When the Boltzmann factor is applied to estimate the solvate populations based upon free energies of the LiClO4(DMCcc)n(DMCct)m, n+m=3 solvates calculated using PBE-GD3BJ functional, we find that approximately a quarter of the DMC in this CIP solvate is in the DMCct conformation. DFT calculations further support the evidence that a substantial fraction of DMCct is present in the first coordination shell of Li+ as obtained in the neutron diffraction with isotopic substitution experiments on DMC-LiClO4 electrolytes.79 In IR experiments, both DMCcc and DMCct contribute to the same carbonyl band position as discussed below, while yielding different position of the Raman active band around 913 cm-1.

29 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 30 of 48

Figure 11. Optimized geometries and relative free energies for LiClO4(DMCcc)m(DMCct)n, m+n=3 solvates from PBE/6-31+G(d,p), PCM(ε=3) calculations. Numbers in parentheses are from DFT calculations with the GD3BJ dispersion correction.

DFT Study of Carbonyl Shift in IR Spectra in Lithium Solvates We follow our previous studies on other Li+ solvates10,

21, 51, 80

and systematically

investigate the influence of Li+ binding on the PC IR vibrational bands in the (Li+)n(PC)m solvates with n=1,3,4 complexes shown in Figure 12 and the PC3(LiBF4) complexes as shown in Figure 8f. Unlike the previous gas-phase QC study of the influence of Li+ binding on the vibrational signature of the solvent, this work incorporates the solvation effects via both the implicit solvent model and explicit solvents. The Li+ solvate size was systematically increased 30 ACS Paragon Plus Environment

Page 31 of 48

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

from 1 PC to 4 PC. Most calculations were performed using B3LYP or M05-2X functionals and a compact 6-31+G(d,p) basis set. A number of additional calculations B3LYP, PBE functionals with the same (6-31+G(d,p)) and much larger aug-cc-pvTz (abbreviated as Tz) basis sets were performed in order to examine the sensitivity of the band shifts and changes of the IR intensity to the choice of density functional and basis set.

Figure 12. Optimized geometries of Li+PC, Li+(PC)3 and Li+(PC)4 from B3LYP/6-31+G(d,p) optimization with PCM(PC).

Shifts of the C=O vibrational band and change of IR intensity of PC upon Li+ complexation from B3LYP calculations are shown in Table 6 for the Li+PC, Li+(PC)3 and Li+(PC)4 clusters surrounding PCM(PC) implicit solvent. Addition of the Li+ cation yields a red shift of the C=O vibrational band by 20-28 cm-1 for the representative larger clusters Li+(PC)3 and Li+(PC)4. This shift is in good agreement with IR spectra. DFT calculations also predict that the Li+ complexation of PC would increase IR intensity of the band by 7-16%. Examination of the LiAnion(PC)3 CIP complexes shown in Figure 10 yielded the C=O band shift between 15 and 25 cm-1, while the intensity of the band increased by 5-8% in PBE/6-31+G(d,p), PCM(PC) calculations.

31 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 32 of 48

Table 6. Shifts (δν) and intensity ratios (I/Ibulk) of the C=O vibrational band of PC upon Li+ complexation from DFT calculations with PCM(PC). PBE/

DFT method

B3LYP/6-31+G(d,p

B3LYP/Tz

Solvate

Li(PC)

Li(PC)3

Li(PC)4

Li(PC)

Li(PC)4

shift δν (cm )

-13

-28

-20

-26

-19, -25

intensity ratio (I/Ibulk)

1.15

1.08

1.09

1.16

1.07,1.07

-1

631+G(d,p)

Analysis of the shifts in the LiAnion(DMC)3 complexes shown in Figure 9 suggests that the carbonyl band shifts by 18-32 cm-1 upon Li+ complexation with the band intensity increasing from 30% to 64% depending on the solvate. A broad range of intensity increases indicates a source of potential uncertainty in extracting the Li+ solvation numbers in DMC from analysis of IR spectra. Moreover, the assumption of the unchanged activity is going to result in overestimation of the number of solvating DMC by 30-60%.

Conclusion Two spectroscopic methods and DFT calculations were utilized to investigate the solvation sphere of the lithium cation. The study was conducted with four lithium salts (LiTFSI, LiPF6, LiBF4, and LiClO4) and two common carbonate solvents (PC and DMC) used in lithiumion battery electrolytes. Solutions were prepared with specific salt-to-solvent mole ratios and the quantities of coordinated and uncoordinated solvent were calculated based on the relative areas of carbonyl IR absorptions characteristic of the coordinated and uncoordinated solvents. Three of the salts, LiTFSI, LiPF6, and LiClO4, are found to readily dissociate in PC at low salt concentrations, while dissociation of LiBF4 occurs to a lesser extent. As the amount of solvent

32 ACS Paragon Plus Environment

Page 33 of 48

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

per salt is decreased, the solvation numbers steadily decrease to ~3 in the most dissociative salts in PC and to ~2 for LiBF4. The solvation numbers for LiBF4 are significantly lower. At moderate concentrations ionic transport in highly dissociated salts LiTFSI, LiPF6, and to a lesser degree for LiClO4 salt, is expected to occur via diffusion of the dissociated ions, while diffusion of charged clusters is expected for LiBF4 in PC. Due to much smaller solvation numbers in DMC compared to PC, all salts are expected to be strongly aggregating with ion transport occurring via diffusion of charged aggregates or structural diffusion via anion – cation dissociation. Ionic aggregation also has ramifications to electrolyte electrochemical stability, LiPF6 and LiTFSI salt aggregation has been shown to increase electrolyte reduction potential due to anion defluorination upon aggregate reduction,

13, 17

while aggregation of LiBF4 is not expected

increase electrolyte reduction potential because the defluorination of LiBF4 aggregates upon reduction occurs at low potentials lower than solvent reduction.13 The solvation structures determined by IR spectroscopy are supported via

13

C NMR

spectroscopy of the carbonyl resonance of the carbonate solvents as a function of salt concentration. The 13C NMR resonance of the carbonyl carbon shifts upon incorporation of the lithium salt due to solvent coordination and rapid solvent exchange. The change in the chemical shift correlates well with the fraction of coordinated solvent determined by IR spectroscopy providing further support for the solvation numbers. DFT calculations using the cluster-continuum approach yield stability of the Li+ solvates in good agreement with the experimentally observed trends for the Li+ solvation numbers. While both mono- and bi-dentate Li+ coordination to TFSI- were predicted in DFT calculations, the LiPF6, LiBF4 and LiClO4 salts favor the monodendtate binding to Li+. Investigation of the Li+(Anion)(DMC)3 solvates showed that both the DMC(cis-cis) and DMC(cis-trans) conformers

33 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 34 of 48

are expected to be present, with an approximate ratio of 3:1. DFT calculations predict shifts of the C=O bond upon Li+ complexation in reasonable agreement with experimental observations. The ratio of the Li+ coordinated PC C=O band intensity to the PC C=O non-coordinated band intensity is predicted to be 1.08-1.16, while this ratio is significantly higher for DMC, 1.3-1.6.

Supporting Information Optimized geometries of Li-anion and Li-solvent configurations, Thermodynamic values obtained by DFT calculations such as binding energy, reaction energy, and free energy

Acknowledgement The authors gratefully acknowledge funding from Department of Energy Office of Basic Energy Sciences EPSCoR Implementation award (DE–SC0007074)

34 ACS Paragon Plus Environment

Page 35 of 48

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

Table of Contents Graphic:

35 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 36 of 48

References 1.

Bockris, J. O. M.; Reddy, A. K., Modern Electrochemistry 2b: Electrodics in Chemistry,

Engineering, Biology and Environmental Science; Springer Science & Business Media, 2001; Vol. 2. 2.

Lucht, B. L.; Markmaitree, T.; Yang, L., Lithium–Ion Batteries. In Encyclopedia of

Inorganic and Bioinorganic Chemistry, John Wiley & Sons, Ltd: 2011. 3.

Xu, K., Electrolytes and Interphases in Li-Ion Batteries and Beyond. Chem. Rev. 2014,

114, 11503-11618. 4.

Seo, D. M.; Borodin, O.; Han, S.-D.; Boyle, P. D.; Henderson, W. A., Electrolyte

Solvation and Ionic Association II. Acetonitrile-Lithium Salt Mixtures: Highly Dissociated Salts. J. Electrochem. Soc. 2012, 159, A1489-A1500. 5.

Barthel, J.; Buchner, R.; Wismeth, E., Ftir Spectroscopy of Ion Solvation of LiClO4 and

LiSCN in Acetonitrile, Benzonitrile, and Propylene Carbonate. J. Solution. Chem. 2000, 29, 937954. 6.

Aroca, R.; Nazri, M.; Nazri, G. A.; Camargo, A. J.; Trsic, M., Vibrational Spectra and

Ion-Pair Properties of Lithium Hexafluorophosphate in Ethylene Carbonate Based MixedSolvent Systems for Lithium Batteries. Journal of Solution Chemistry 2000, 29, 1047-1060. 7.

Burba, C. M.; Frech, R., Spectroscopic Measurements of Ionic Association in Solutions

of LiPF6. J. Phys. Chem. B 2005, 109, 15161-15164. 8.

Takeuchi, M.; Kameda, Y.; Umebayashi, Y.; Ogawa, S.; Sonoda, T.; Ishiguro, S.-i.;

Fujita, M.; Sano, M., Ion–Ion Interactions of LiPF6 and LiBF4 in Propylene Carbonate Solutions. J. Moleq. Liquids 2009, 148, 99-108.

36 ACS Paragon Plus Environment

Page 37 of 48

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

9.

Nie, M.; Abraham, D. P.; Seo, D. M.; Chen, Y.; Bose, A.; Lucht, B. L., Role of Solution

Structure in Solid Electrolyte Interphase Formation on Graphite with LiPF6 in Propylene Carbonate. J. Phys. Chem. C 2013, 117, 25381-25389. 10.

Grande, L.; von Zamory, J.; Koch, S. L.; Kalhoff, J.; Paillard, E.; Passerini, S.,

Homogeneous Lithium Electrodeposition with Pyrrolidinium-Based Ionic Liquid Electrolytes. ACS Applied Materials & Interfaces 2015, 7, 5950-5958. 11.

Seo, D. M.; Reininger, S.; Kutcher, M.; Redmond, K.; Euler, W. B.; Lucht, B. L., Role of

Mixed Solvation and Ion Pairing in the Solution Structure of Lithium Ion Battery Electrolytes. J. Phys. Chem. C 2015, 119, 14038-14046. 12.

Xu, K.; Lam, Y.; Zhang, S. S.; Jow, T. R.; Curtis, T. B., Solvation Sheath of Li+ in

Nonaqueous Electrolytes and Its Implication of Graphite/Electrolyte Interface Chemistry. J. Phys. Chem. C 2007, 111, 7411-7421. 13.

Delp, S. A.; Borodin, O.; Olguin, M.; Eisner, C. G.; Allen, J. L.; Jow, T. R., Importance

of Reduction and Oxidation Stability of High Voltage Electrolytes and Additives. Electrochimica Acta 2016, 209, 498-510. 14.

von Wald Cresce, A.; Borodin, O.; Xu, K., Correlating Li+ Solvation Sheath Structure

with Interphasial Chemistry on Graphite. J. Phys. Chem. C 2012, 116, 26111-26117. 15.

von Cresce, A.; Xu, K., Preferential Solvation of Li+ Directs Formation of Interphase on

Graphitic Anode. Electrochem. and Solid-State Lett. 2011, 14, A154-A156. 16.

Xu, K.; von Cresce, A., Interfacing Electrolytes with Electrodes in Li Ion Batteries. J

Mater Chem 2011, 21, 9849-9864. 17.

Suo, L.; Borodin, O.; Sun, W.; Fan, X.; Yang, Wang, F.; Gao, T.; Ma, Z.; Schroeder, M.;

von Cresce, A.; Russell, S. M.; Armand, M.; Angell, A.; Xu K.; Wang, C.; Advanced High-

37 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 38 of 48

Voltage Aqueous Lithium-Ion Battery Enabled by “Water-in-Bisalt” Electrolyte. Angew. Chem. Int. Ed. 2016, 55, 7136-7141. 18.

Kim, H.; Wu, F.; Lee, J. T.; Nitta, N.; Lin, H.-T.; Oschatz, M.; Cho, W. I.; Kaskel, S.;

Borodin, O.; Yushin, G., In Situ Formation of Protective Coatings on Sulfur Cathodes in Lithium Batteries with Lifsi-Based Organic Electrolytes. Adv. Energ. Mater. 2015, 5, 1401792. 19.

Tarascon, J. M.; Guyomard, D., New Electrolyte Compositions Stable over the 0 to 5 V

Voltage Range and Compatible with the Li1+Xmn2o4/Carbon Li-Ion Cells. Solid State Ionics 1994, 69, 293-305. 20.

Soetens, J. C.; Millot, C.; Maigret, B., Molecular Dynamics Simulation of Li+BF4- in

Ethylene Carbonate, Propylene Carbonate, and Dimethyl Carbonate Solvents. J. Phys. Chem. A 1998, 102, 1055-1061. 21.

Borodin, O.; Olguin, M.; Ganesh, P.; Kent, P. R. C.; 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. 22.

von Wald Cresce, A., et al., Anion Solvation in Carbonate-Based Electrolytes. J. Phys.

Chem. C 2015, 119, 27255-27264. 23.

Seo, D. M.; Boyle, P. D.; Sommer, R. D.; Daubert, J. S.; Borodin, O.; Henderson, W. A.,

Solvate Structures and Spectroscopic Characterization of Litfsi Electrolytes. J. Phys. Chem. B 2014, 118, 13601-13608. 24.

Han, S. D.; Yun, S. H.; Borodin, O.; Seo, D. M.; Sommer, R. D.; Young, V. G.;

Henderson, W. A., Solvate Structures and Computational/Spectroscopic Characterization of Lipf6 Electrolytes. J. Phys. Chem. C 2015, 119, 8492-8500.

38 ACS Paragon Plus Environment

Page 39 of 48

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

25.

Smith, J. W.; Lam, R. K.; Sheardy, A. T.; Shih, O.; Rizzuto, A. M.; Borodin, O.; Harris,

S. J.; Prendergast, D.; Saykally, R. J., X-Ray Absorption Spectroscopy of LiBF4 in Propylene Carbonate: A Model Lithium Ion Battery Electrolyte. Phys. Chem. Chem. Phys. 2014, 16, 2356823575. 26.

Bogle, X.; Vazquez, R.; Greenbaum, S.; Cresce, A. v. W.; Xu, K., Understanding Li+–

Solvent Interaction in Nonaqueous Carbonate Electrolytes with

17

O NMR. J. Phys. Chem. Lett.

2013, 4, 1664-1668. 27.

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.; Xu, K.; Solvation Behavior of CarbonateBased Electrolytes in Sodium Ion Batteries. Phys. Chem. Chem. Phys. 2017, 19, 574-586. 28.

Takeuchi, M.; Matubayasi, N.; Kameda, Y.; Minofar, B.; Ishiguro, S.-i.; Umebayashi, Y.,

Free-Energy and Structural Analysis of Ion Solvation and Contact Ion-Pair Formation of Li+ with BF4– and PF6– in Water and Carbonate Solvents. J. Phys. Chem. B 2012, 116, 6476-6487. 29.

Li, T.; Balbuena, P. B., Theoretical Studies of Lithium Perchlorate in Ethylene

Carbonate, Propylene Carbonate, and Their Mixtures. J. Electrochem. Soc. 1999, 146, 36133622. 30.

Johansson, P., Electronic Structure Calculations on Lithium Battery Electrolyte Salts.

Phys Chem Chem Phys 2007, 9, 1493-1498. 31.

Tachikawa, H., Mechanism of Dissolution of a Lithium Salt in an Electrolytic Solvent in

a Lithium Ion Secondary Battery: A Direct Ab Initio Molecular Dynamics (Aimd) Study. Chemphyschem 2014, 15, 1604-1610.

39 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

32.

Page 40 of 48

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, 45024515. 33.

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. 34.

Popov, S. E.; Nikiforov, A. E.; Bushkova, O. V.; Zhukovskii, V. M., Quantum-Chemical

Investigation of Ionic Association of Lithium Salts LiXF6 (X = As, P). Russian Journal of Electrochemistry 2005, 41, 476-484. 35.

Eilmes, A.; Kubisiak, P., Stability of Ion Triplets in Ionic Liquid/Lithium Salt Solutions:

Insights from Implicit and Explicit Solvent Models and Molecular Dynamics Simulations. J. Comput. Chem. 2015, n/a-n/a. 36.

Malliakas, C. D.; Leung, K.; Pupek, K. Z.; Shkrob, I. A.; Abraham, D. P., Spontaneous

Aggregation of Lithium Ion Coordination Polymers in Fluorinated Electrolytes for High-Voltage Batteries. Phys. Chem. Chem. Phys. 2016. 37.

Masia, M.; Rey, R., Computational Study of Gamma-Butyrolactone and Li+/Gamma-

Butyrolactone in Gas and Liquid Phases. J. Phys. Chem. B 2004, 108, 17992-18002. 38.

Masia, M.; Probst, M.; Rey, R., Ethylene Carbonate-Li+: A Theoretical Study of

Structural and Vibrational Properties in Gas and Liquid Phases. J. Phys. Chem. B 2004, 108, 2016-2027. 39.

Silva, L. B.; Freitas, L. C. G., Structural and Thermodynamic Properties of Liquid

Ethylene Carbonate and Propylene Carbonate by Monte Carlo Simulations. J. Mol. Struct.Theochem 2007, 806, 23-34.

40 ACS Paragon Plus Environment

Page 41 of 48

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

40.

Borodin, O.; Smith, G. D., Litfsi Structure and Transport in Ethylene Carbonate from

Molecular Dynamics Simulations. J. Phys. Chem. B 2006, 110, 4971-4977. 41.

Borodin, O.; Smith, G. D., Development of Many-Body Polarizable Force Fields for Li-

Battery Components: 1. Ether, Alkane, and Carbonate-Based Solvents. J. Phys. Chem. B 2006, 110, 6279-6292. 42.

Borodin, O.; Smith, G. D., Development of Many-Body Polarizable Force Fields for Li-

Battery Applications: 2. LiTFSI-Doped Oligoether, Polyether, and Carbonate-Based Electrolytes. J. Phys. Chem. B 2006, 110, 6293-6299. 43.

Borodin, O.; Smith, G. D., Li+ Transport Mechanism in Oligo(Ethylene Oxide)S

Compared to Carbonates. J. Solution. Chem. 2007, 36, 803-813. 44.

Tasaki, K.; Nakamura, S., Computer Simulation of LiPF6 Salt Association in Li-Ion

Battery Electrolyte in the Presence of an Anion Trapping Agent. J. Electrochem. Soc. 2001, 148, A984-A988. 45.

Newman, J.; Thomas, K. E.; Hafezi, H.; Wheeler, D. R., Modeling of Lithium-Ion

Batteries. J. Power Sources 2003, 119, 838-843. 46.

McOwen, D. W.; Seo, D. M.; Borodin, O.; Vatamanu, J.; Boyle, P. D.; Henderson, W.

A., Concentrated Electrolytes: Decrypting Electrolyte Properties and Reassessing Al Corrosion Mechanisms. Ener.& Env. Sci. 2014, 7, 416-426. 47.

Leung, K.; Leenheer, A. J., How Voltage Drops Are Manifested by Lithium Ion

Configurations at Interfaces and in Thin Films on Battery Electrodes. J. Phys. Chem. C 2015. 48.

Korsun, O. M.; Kalugin, O. N.; Fritsky, I. O.; Prezhdo, O. V., Ion Association in Aprotic

Solvents for Lithium Ion Batteries Requires Discrete-Continuum Approach: Lithium

41 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 42 of 48

Bis(Oxalato)Borate in Ethylene Carbonate Based Mixtures. J. Phys. Chem. C 2016, 120, 1654516552. 49.

Kondo, K.; Sano, M.; Hiwara, A.; Omi, T.; Fujita, M.; Kuwae, A.; Iida, M.; Mogi, K.;

Yokoyama, H., Conductivity and Solvation of Li+ Ions of LiPF6 in Propylene Carbonate Solutions. J. Phys. Chem. B 2000, 104, 5040-5044. 50.

Xu, K., Nonaqueous Liquid Electrolytes for Lithium-Based Rechargeable Batteries.

Chem. Rev. 2004, 104, 4303-4418. 51.

Seo, D. M.; Borodin, O.; Han, S.-D.; Ly, Q.; Boyle, P. D.; Henderson, W. A., Electrolyte

Solvation and Ionic Association. J. Electrochem. Soc. 2012, 159, A553-A565. 52.

Campion, C. L.; Li, W.; Lucht, B. L., Thermal Decomposition of LiPF6-Based

Electrolytes for Lithium-Ion Batteries. J. Electrochem. Soc. 2005, 152, A2327-A2334. 53.

Bogle, X.; Vazquez, R.; Greenbaum, S.; Cresce, A. v. W.; Xu, K., Understanding Li+–

Solvent Interaction in Nonaqueous Carbonate Electrolytes with

17

O NMR. J. Phys. Chem. Lett.

2013, 4, 1664-1668. 54.

Reichardt, C.; Welton, T., Solvents and Solvent Effects in Organic Chemistry; John Wiley

& Sons, 2011. 55.

Prassides, K., Mixed Valency Systems: Applications in Chemistry, Physics and Biology;

Springer Science & Business Media, 2012; Vol. 343. 56.

Yanase, S.; Oi, T., Solvation of Lithium Ion in Organic Electrolyte Solutions and Its

Isotopie Reduced Partition Function Ratios Studied by Ab Initio Molecular Orbital Method. Journal of Nuclear Science and Technology 2002, 39, 1060-1064.

42 ACS Paragon Plus Environment

Page 43 of 48

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

57.

Yuan, K.; Bian, H.; Shen, Y.; Jiang, B.; Li, J.; Zhang, Y.; Chen, H.; Zheng, J.,

Coordination Number of Li+ in Nonaqueous Electrolyte Solutions Determined by Molecular Rotational Measurements. J. Phys. Chem. B 2014, 118, 3689-3695. 58.

Ue, M.; Mori, S., Mobility and Ionic Association of Lithium Salts in a Propylene

Carbonate‐Ethyl Methyl Carbonate Mixed Solvent. J. Electrochem. Soc. 1995, 142, 2577-2581. 59.

Nishikawa, K.; Fukunaka, Y.; Sakka, T.; Ogata, Y. H.; Selman, J. R., Measurement of

LiClO4 Diffusion Coefficient in Propylene Carbonate by Moiré Pattern. J. Electrochem. Soc. 2006, 153, A830-A834. 60.

Nishida, T.; Nishikawa, K.; Fukunaka, Y., Diffusivity Measurement of LiPF6, LiTFSI,

LiBF4 in PC. ECS Transactions 2008, 6, 1-14. 61.

Henderson, W. A., Glyme−Lithium Salt Phase Behavior. J. Phys. Chem. B 2006, 110,

13177-13183. 62.

Allen, J. L.; Borodin, O.; Seo, D. M.; Henderson, W. A., Combined Quantum

Chemical/Raman Spectroscopic Analyses of Li+ Cation Solvation: Cyclic Carbonate Solvents— Ethylene Carbonate and Propylene Carbonate. J. Power Sources 2014, 267, 821-830. 63.

Hayamizu, K.; Aihara, Y.; Arai, S.; Martinez, C. G., Pulse-Gradient Spin-Echo H-1, Li-7,

and F-19 Nmr Diffusion and Ionic Conductivity Measurements of 14 Organic Electrolytes Containing LiN(SO2CF3)(2). J. Phys. Chem. B 1999, 103, 519-524. 64.

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.

43 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

65.

Page 44 of 48

Marenich, A. V.; Cramer, C. J.; Truhlar, D. G., Universal Solvation Model Based on

Solute Electron Density and on a Continuum Model of the Solvent Defined by the Bulk Dielectric Constant and Atomic Surface Tensions. J. Phys. Chem. B 2009, 113, 6378-6396. 66.

Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman,

J. R.; Scalmani, G., Gaussian 09, Revision C1. 2010. 67.

Bauschlicher, C. W.; Haskins, J. B.; Bucholz, E. W.; Lawson, J. W.; Borodin, O.,

Structure and Energetics of Li+–(BF4–)N, Li+–(FSI–)N, and Li+–(TFSI–)N: Ab Initio and Polarizable Force Field Approaches. J. Phys. Chem. B 2014, 118, 10785-10794. 68.

Bryantsev, V. S., Calculation of Solvation Free Energies of Li+ and O2 − Ions and

Neutral Lithium–Oxygen Compounds in Acetonitrile Using Mixed Cluster/Continuum Models. Theor. Chem. Acc. 2012, 131, 1250. 69.

Seo, D. M.; Boyle, P. D.; Borodin, O.; Henderson, W. A., Li+ Cation Coordination by

Acetonitrile-Insights from Crystallography. RSC Advances 2012, 2, 8014-8019. 70.

Borodin, O.; Behl, W.; Jow, T. R., Oxidative Stability and Initial Decomposition

Reactions of Carbonate, Sulfone, and Alkyl Phosphate-Based Electrolytes. J. Phys. Chem. C 2013, 117, 8661-8682. 71.

Borodin, O., Molecular Modeling of Electrolytes. In Electrolytes for Lithium and

Lithium-Ion Batteries, Jow, T. R.; Xu, K.; Borodin, O.; Ue, M., Eds. Springer New York: 2014; Vol. 58, pp 371-401. 72.

Ding, W.; Lei, X.; Ouyang, C., Coordination of Lithium Ion with Ethylene Carbonate

Electrolyte Solvent: A Computational Study. Int. J. Quant. Chem. 2016, 116, 97-102. 73.

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

44 ACS Paragon Plus Environment

Page 45 of 48

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

Theory with the Hybrid Solvation Model and Free Energy Correction in Solution. Phys. Chem. Chem. Phys. 2016, 18, 23607-23612. 74.

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. 75.

Borodin, O.; Smith, G. D., Litfsi Structure and Transport in Ethylene Carbonate from

Molecular Dynamics Simulations. J. Phys. Chem. B 2006, 110, 4971-4977. 76.

Borodin, O.; Smith, G. D.; Jaffe, R. L., Ab Initio Quantum Chemistry and Molecular

Dynamics Simulations Studies of LiPF6/Poly(Ethylene Oxide) Interactions J. Comput. Chem. 2001, 22, 641-654. 77.

Barnes, T. A.; Kaminski, J. W.; Borodin, O.; Miller, T. F., Ab Initio Characterization of

the Electrochemical Stability and Solvation Properties of Condensed-Phase Ethylene Carbonate and Dimethyl Carbonate Mixtures. J. Phys. Chem. C 2015, 119, 3865-3880. 78.

Grimme, S.; Ehrlich, S.; Goerigk, L., Effect of the Damping Function in Dispersion

Corrected Density Functional Theory. J. Comput. Chem. 2011, 32, 1456-1465. 79.

Kameda, Y.; Saito, S.; Umebayashi, Y.; Fujii, K.; Amo, Y.; Usuki, T., Local Structure of

Li+ in Concentrated LiPF6–Dimethyl Carbonate Solutions. J. Mol. Liq. 2016, 217, 17-22. 80.

Allen, J. L.; Borodin, O.; Seo, D. M.; Henderson, W. A., Combined Quantum

Chemical/Raman Spectroscopic Analyses of Li+ Cation Solvation: Cyclic Carbonate Solvents— Ethylene Carbonate and Propylene Carbonate. Journal of Power Sources 2014, 267, 821-830.

45 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

338x190mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 46 of 48

Page 47 of 48

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

338x190mm (300 x 300 DPI)

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

587x577mm (96 x 96 DPI)

ACS Paragon Plus Environment

Page 48 of 48