* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Download PDF - ACS Publications - American Chemical Society
Survey
Document related concepts
Thermodynamics wikipedia , lookup
2-Norbornyl cation wikipedia , lookup
Homoaromaticity wikipedia , lookup
Ultraviolet–visible spectroscopy wikipedia , lookup
Equilibrium chemistry wikipedia , lookup
Stability constants of complexes wikipedia , lookup
Ionic compound wikipedia , lookup
Ionic liquid wikipedia , lookup
History of electrochemistry wikipedia , lookup
Determination of equilibrium constants wikipedia , lookup
Electrolysis of water wikipedia , lookup
Transcript
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes. Research Article http://pubs.acs.org/journal/acscii Transport Properties of Polysulfide Species in Lithium−Sulfur Battery Electrolytes: Coupling of Experiment and Theory M. Safari,† C. Y. Kwok, and L. F. Nazar* Department of Chemistry, University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada S Supporting Information * ABSTRACT: A comprehensive experimental and theoretical analysis of the isothermal transport of species for the two model ternary-electrolytes with LiTFSILi2S4/dioxolane (DOL)-dimethoxyethane (DME) and LiTFSI-Li2S6/DOL-DME formulations is presented. An unambiguous picture of the polysulfide’s mobility is set forth after a detailed investigation of the macroscopic transference number and diffusion coefficients. The new findings of incongruent diffusion for Li2S4 species and high significance of cross-term diffusion coefficients reformulate a fledgling view of the prevalent redox-shuttle phenomena. The practical implications of this complex mechanism are discussed in detail. ■ INTRODUCTION Successful penetration of lithium−sulfur (Li−S) technology into the rechargeable battery market could serve the thirst for higher energy density battery systems in the electric-transport sector. These batteries that utilize lithium metal as the negative electrode and sulfur as positive electrode benefit from high theoretical specific capacity and energy density compared to lithium-ion batteries, coupled with low cost.1−3 However, before their practical realization, there are many hurdles to overcome such as the insulating character of the active material end members (sulfur and Li2S) on discharge and charge;4 the corresponding large volume change;5 lithium dendrite formation6 during cell operation; andmost importantly dissolution of the intermediate lithium polysulfide redox species into the electrolyte.7 Unlike conventional Li-ion battery technologies, in a Li-sulfur battery, formation of a stable anode/ electrolyte interface is mainly challenged by the intermediate species produced at the cathode. Understanding the underlying mechanisms of operation in Li−S cells especially remains a major obstacle to their continued improvement.8 Extensive research is currently ongoing to minimize lithium polysulfide shuttles in Li−S batteries and to eliminate all its concomitant risks to the battery state of health. Much effort has been expanded in attempts to confine the soluble polysulfides via modification of cathode and cell design.9−11 Some of these approaches include physical entrapment12,13 or chemical adsorption14,15 of polysulfides within the cathode network. Moreover, electrochemical and operando spectroscopic studies reveal that the exact composition of these intermediates is solvent dependent and can be influenced by a series of chemical disproportionation reactions in the electrolyte.16,17 Lu et al. used a rotating-ring disk electrode technique to probe the kinetics of the electrochemical reactions in Li−S cells,18 while other researchers19−21 focused more on the dominant soluble polysulfide species. These are effectively Li2S6 and Li2S4 in solvents of intermediate polarity such as dioxolane/dimethoxyethane (1:1 DOL/DME), which serves as the classical Li−S © 2016 American Chemical Society battery electrolyte, where lithium bis-trifluoromethane-sulfonimide (LiTFSI) is the conducting salt. Many studies show that Li2S8 is unstable toward disproportionation to these species, and Li2S2 is either unstable and/or insoluble.20 In highly polar media such as dimethyl amide (DMA) and dimethyl sulfoxide (DMSO), the hexa-sulfide and tetra-sulfide dianions undergo cleavage to their monoanion radicals. In such media, S3−• the product of S6−2 cleavage is especially prominent and is strongly favored in dilute solutions. This gives rise to the characteristic deep blue color of such solutions.22,23 However, cleavage in DOL/DME is not favored at room temperature, and hence the concentration of the trisulfide radical is neglible.24,25 Despite these many experimental studies to elucidate polysulfide speciation, very limited modeling work on Li−S batteries has been conducted that could provide an understanding of polysulfide transport properties, and hence their behavior in the cell. Only a few seminal and important studies have been reported,7,26,27 contrary to lithium ion battery systems.28−31 An accurate picture of species transport in the electrolyte is not only of fundamental importance, but is a prerequisite for any meaningful macroscopic modeling/optimization of the cell performance in a Li−S battery. 32−34 In this regard, experimental measure of a minimum number of the iontransference numbers and diffusion coefficients for neutral combination of the ions, together with the ionic conductivity of the electrolyte is required. Here, for the first time, we determine a complete set of transport properties for the binary electrolytes composed of Li2S6, Li2S4 (and the electrolyte salt, LiTFSI) dissolved in DOL/DME (1:1). These measurements are performed according to the original approach developed by Ma and Newman.35 The theoretical/experimental frames provided by the thermodynamics of irreversible processes are further applied to estimate the transport properties for the more Received: June 9, 2016 Published: July 28, 2016 560 DOI: 10.1021/acscentsci.6b00169 ACS Cent. Sci. 2016, 2, 560−568 Research Article ACS Central Science with a common cation (C) and A1 and A2 as two anions can be represented as the following: practical ternary cases where the polysulfide and lithium salt are both present in the electrolyte.36 ■ RESULTS The restricted diffusion (RD), transference-polarization (TP), and concentration-cell (CC) experiments were all carried out on 1 mL of the electrolyte confined inside a Teflon (PFA) Swagelok cell. A schematic of this symmetric cell is presented in Figure 1 where a 12 mm PFA tube with an internal diameter of C ϑ1,cA1ϑ1,a ↔ ϑ1,cCz1 + ϑ1,aA1z 2 (1) C ϑ2,cA2 ϑ2,a ↔ ϑ2,cCz1 + ϑ2,aA2z3 (2) where, ϑi,c and ϑi,a are the stoichiometries of the cation and anions in the ith solute (i = 1, 2), and the charge numbers of the common cation and two anions are represented by z1, z2, and z3, respectively. Isothermal transport of species in such a system is described by two independent transference numbers, four diffusion coefficients for neutral combinations of ions (out of which three are independent) and the ionic conductivity of the electrolyte. These measurable macroscopic properties are expressed with the aid of nine fundamental transport coefficients (l) where only six are independent according to the Onsager reciprocal relations.36 In a solvent-fixed frame (SF), hereafter represented by 0 in the superscript, the flux equations for the ions (i.e., J0C, J0A1, J0A2) in one dimension are t10 0 I − (ϑ1,c+ 022 + ϑ2,c+ 32 )∇μ12 Fz1 JC0 = J10 = Figure 1. Schematic of the Swagelok cell used in this study for the electrochemical characterization of the electrolytes. 0 − (ϑ1,c+ 023 + ϑ2,c+ 33 )∇μ13 10 mm is positioned between the two lithium electrodes and filled with the electrolyte. In the concentration-cell (CC) experiments, the internal PFA tube is divided in two compartments separated with a 3 mm ultrafine porous glass frit to avoid mixing among the two compartments. The binary and ternary electrolytes with Li2Sn (n = 6, 4) and LiTFSI as solutes were prepared in a 1,3 dioxolane/1,2 dimethoxyethane (DOL/DME) solvent (1−1 volume ratio) over a concentration range of 0.05−1.1 M in solute. Additional details are provided in Experimental Methods (see below). Development of the Theoretical Framework. In a Li−S battery, the concentrated binary electrolyte (i.e., a primary lithium salt dissolved in a solvent) upon assembly turns into an electrolyte with a single cation (Li+) and multiple anions (i.e., polysulfides in addition to the primary lithium-salt anion) over the course of battery operation.20 Hence, in a Li−S battery and in the presence of n independent series of polysulfides, (n + 3)(n + 2) transport properties are required to describe the 2 isothermal transport of species in the electrolyte. This means that in the simplest case where only one of the high-order lithium-polysulfides (Li2Sn, 4 ≤ n ≤ 8) is dominant, one is left with a ternary electrolyte and six transport properties to be measured. Here, our target is to provide consistent estimates of diffusion coefficients and transference numbers for a typical electrolyte used in Li−S cells. To do so, first, an appropriate formulation of species transport is presented for a general ternary electrolyte composed of a primary lithium salt, a single lithium-polysulfide (Li2Sn), and a solvent. Further, the two approaches of Miller37,38 and Ma et al.35 are combined and adopted to the ternary electrolytes of LiTFSI and Li2S6 (Li2S4) dissolved in DOL/DME (namely, the binary coefficients were determined with the aid of three orthogonal experiments as set forth by Ma et al., and the ternary data are estimated according to Miller et al.). Fundamental Transport Coefficients in a Ternary Electrolyte with a Common Cation. A ternary electrolyte (3) 0 JA1 = J20 = t 20 0 0 I − ϑ1,a+ 22 ∇μ12 − ϑ1,a+ 23 ∇μ13 Fz 2 (4) 0 JA2 = J30 = t30 0 0 I − ϑ2,a+ 32 ∇μ12 − ϑ2,a+ 33 ∇μ13 Fz 3 (5) where I is the total current density and F is the Faraday constant. In these flux equations, ∇μ12 and ∇μ13 are the electrochemical potential gradients of the neutral combination of the common cation (denoted as subscript 1) with the first (2) and second (3) anions, respectively, and are related to their concentration counterparts (i.e., ∇c12, ∇c13) according to the following equation:39 ⎡ ∂ ln f±1i ϑ1,cϑi − 1,c ∇μ1i = RT ⎢(ϑi − 1,c + ϑi − 1,a) + ⎢⎣ ∂ c12 ϑ1,cc12 + ϑ2,cc13 + + ϑi − 1,a c1i ⎤ ⎡ ∂ ln f±1i ⎥ ⎢ δ1, i − 1 ∇c12 + RT (ϑi − 1,c + ϑi − 1,a) ⎥⎦ ⎢⎣ ∂ c13 ϑ2,cϑi − 1,c ϑ1,cc12 + ϑ2,cc13 + ϑi − 1,a c1i ⎤ δ2, i − 1⎥∇c13 , (i = 2, 3) ⎥⎦ (6) where δ is the Kronecker delta function, R is the universal gas constant, and T is temperature. The molar concentrations and mean-molar ionic activity coefficients of the two solutes (i.e., Cϑ1,cA1ϑ1,a and Cϑ2,cA2ϑ2,a) are represented by (c12, c13) and 13 ( f 12 ± , f ± ), respectively. The transference numbers (t0) and fundamental diffusion coefficients (+ 0 ) are expressed with the aid of more fundamental transport coefficients (l0) set forth by the thermodynamics of irreversible processes according to eqs 7 and 8, respectively.36−38 k=3 ti0 = 561 F 2zi ∑k = 1 zklik0 κ , i = 1, 2, 3 (7) DOI: 10.1021/acscentsci.6b00169 ACS Cent. Sci. 2016, 2, 560−568 Research Article ACS Central Science +ij0 = F2 κ 3 3 ∑∑ k=1 w=1 N = |ϑ1,az 2c12 + ϑ2,az 3c13| = z1(ϑ1,cc12 + ϑ2,cc13) z k zw 0 − liw0 lkj0 ), i , j = 2, 3 (lij0lkw ϑi − 1,aϑj − 1,a In eqs 13−16, x12 and x13 are the equivalent fractions and are defined as the following: (8) where, in eq 7 and 8, κ is the electrolyte conductivity: 3 κ = F2 ∑ (9) k=1 w=1 3 ∑ Dij∇c1j = ∑ +ij∇μ1j , i = 2, 3 j=2 (10) j=2 (18) N1i = |ϑi − 1,azic1i| , i = 2, 3 (19) Transference numbers, conductivity, and fundamental diffusion coefficients (+ij0 ) for the ternary system (eqs 1 and 2) are readily available (eqs 7−9) after estimation of l0ij. Derivatives of mean-molar ionic activity coefficients ( f12 ± and f13 ± ) are required, however, in order to obtain practical diffusion coefficients (D ij0 ). Providing that these thermodynamic parameters are known for the binary constituents ( f12,b and ± f13,b ), the Scatchard’s neutral-electrolyte description can be used ± 42 13 to estimate f12 ± and f± in the ternary system: However, the transference numbers and diffusion coefficients are dependent on the chosen frame of reference. Transference numbers are experimentally measured in a solvent-fixed (SF) frame, whereas the practical diffusion coefficients are usually measured in a volume-fixed frame of reference (VF). Practical diffusion coefficients (D) and fundamental ones (+ ) are related according to eq 10 3 N1i , i = 2, 3 N x1i = 3 ∑ zkzwlkw0 ln f±1i It is more convenient to set forth the transport equations in SF and further use the following equation to transform the practical diffusion coefficients between the two frames of reference (i.e., SF and VF):40 ∑ (δki − c1iV1̅ k)Dkj0 , i , j = 2, 3 k=2 V1̅ i = 3 ∂ρ ∂c1i ∂ρ ρ − ∑ j = 2 c1j ∂c (11) − 0 0,b 0,b l11 = x12(l++ )12, N + x13(l++ )13, N (14) 0,b l10i = x1i(l+− )1i , N , i = 2, 3 (15) 0 0 0 0.5 l 23 = x12x13(l12 l13) 0,b l0,b + −, l+ +, (21) where c0 and are the solvent concentrations in the ternary and binary electrolytes, respectively, ρ0 is the solvent density, and M0 is molar mass of the solvent. In eqs 20 and 21, all the binary parameters are evaluated at the total ionic-strength of the ternary electrolyte and y defines the ionic-strength fraction42 of each solute in the ternary electrolyte according to where M12 and M13 are the molar masses of Cϑ1,cA1ϑ1,a and Cϑ2,cA2ϑ2,a, respectively, and ρ is the density of the ternary electrolyte. Note that in Miller’s original treatment of ternary electrolytes with a common anion,38 he distinguished between the chemical and electrical potentials in the entropy production term. This separation is unnecessary,39 though similar results are obtained when the thermodynamic force is defined in terms of the electrochemical potentials. Estimation of Ternary Transport Properties from Experimental Binary Data. The six fundamental transport coefficients (l0ij) for the system defined by eqs 1 and 2 are estimated based on the measured coefficients (l0,b ij ) for the two isolated binary systems, defined by eqs 1 and 2, and application of the following correlations after Miller:38 (13) cM 1 ln 0 0 , i = 2, 3 z1zi ρ0 (cb0)1i (12) 0,b lii0 = x1i(l −− )1i , N , i = 2, 3 (20) b y + δi3 (c0b)13 M 0 ⎤ 1 ⎡ y12 + δi2 (c0 )12 M 0 ⎢ ⎥ ln ln + 13 z1z 3 2 ⎢⎣ z1z 2 ρ0 ρ0 ⎥⎦ βi = , i = 2, 3 1j 1i ,b ⎡ ln f 12,b ln f±13,b ⎤ 1 ln f± 1 ± ⎥ + β, + ⎢y12 + y13 i z1z 2 z1z 3 ⎥⎦ 2 z1zi 2 ⎢⎣ i = 2, 3 where V̅ 12 and V̅ 13 are the partial molar-volumes of Cϑ1,cA1ϑ1,a and Cϑ2,cA2ϑ2,a, respectively, and are experimentally measured with the aid of the following equation:41 M1i − = z1zi 3 Dijv = (17) y1i = c1i ϑi − 1, cz1(z1 − zi) , 3 ∑ j = 2 c1j ϑj − 1,cz1(z1 − zj) i = 2, 3 (22) Fundamental Transport Coefficients for the Binary Electrolytes of LiTFSI and Li2Sn (n = 6, 4) in DOL/DME. 0,b Here, we are interested at experimental measure of l0,b + −, l+ +, and 0,b l− − for the three binary systems described by eqs 23−25 Li TFSI ↔ Li+1 + TFSI−1, (solvent: DOL/DME) (23) Li 2S6 ↔ 2Li+1 + S6−2 , (solvent: DOL/DME) (24) Li 2S4 ↔ 2Li+1 + S4 −2 , (solvent: DOL/DME) (25) It is straightforward to show that in a binary electrolyte, the three fundamental transport coefficients can be expressed as a function of three macroscopic and measurable transport V,b properties (i.e., κb, t0,b i , and D ) in eq 26 lij0,b = (16) κ bti0,bt j0,b F 2zizj + ( −1)δij ϑ+ϑ−z+z −c bDV,b , i , j = +, − ⎛ ∂ln f±b ⎞ b b ⎟c V ̅ RTzizj(ϑ+ + ϑ−)⎜1 + ∂ln c b ⎠ 0 0 ⎝ l0,b −− where and are the cation−anion, cation−cation, and anion−anion interaction coefficients in the binary electrolytes, respectively. These binary coefficients are evaluated at a total (N) number of equivalents in the ternary electrolyte: (26) 562 DOI: 10.1021/acscentsci.6b00169 ACS Cent. Sci. 2016, 2, 560−568 Research Article ACS Central Science Figure 2. Transport coefficients: (a) practical diffusion coefficient in a volume-fixed frame of reference, (b) transference number of lithium ions in a solvent-fixed frame of reference (c) electrolyte ionic conductivityin the binary electrolytes of LiTFSI-DOL/DME (dash-dot-○-), Li2S4-DOL/ DME (dot-△-), and Li2S6-DOL/DME (dash-*-) at T = 25 °C. Markers are experimental data and lines are shape-preserving fits only for improved visualization. where V̅ 0b and c 0b are the partial molar-volume and concentration of the solvent, respectively, and cb is the solute concentration in the binary electrolyte. The value κb is readily attained with the application of ACimpedance spectroscopy. The high frequency intercept on the x (real)-axis in a Nyquist plot is a sufficiently accurate measure of the ohmic resistance of the electrolyte. This resistance, together with the known geometry of the cell, provides us with the specific ionic conductivity of the electrolyte which is independent of the frame of reference used during the measurement. However, the experimental measurement of t0,b i and DV,b is not trivial. In the following, the elegant approach originally developed by Ma et al.35 and lately improved by Hafezi and Newman43 is applied to the case of binary electrolytes investigated in this study (eqs 23−25). Here, and for the sake of simplicity, the solvent system (DOL/DME) is treated as a single component. The transference number of lithium ions (t0,b + ) and the mean-molar ionic activity coefficient of the solutes ( fb±) are coupled together according to eq 27 in the absence of current (I = 0)41 ∇Φ = different concentrations of a given solute at the left (cbl ) and right (cbr ) compartments. Transference-polarization (TP) tests are performed to compensate for the dearth of data regarding the thermodynamic factor of the solute (which could be obtained independently via vapor pressure or isopiestic measurements with great effort). In TP experiments, a short pulse of current is applied to an equilibrated symmetric cell made of the lithium electrodes and the binary electrolyte under study at a given initial concentration (cb∞). The difference of solute concentration (Δcb) between the vicinity of the two electrodes created by the current pulse is43 Δc b = 4(1 − t+0,b)I D∞V,bt p F ϑ+z+D∞V,b π τ (28) where tp is the pulse duration and DV,b ∞ is the diffusion coefficient of solute in the electrolyte at cb = cb∞ and in VF. τ is a dimensionless time and is defined by eq 29: τ= ⎛ ∂ ln f±b ⎞ (ϑ+ + ϑ−)RT ⎜ ⎟(1 − t 0,b)∇ln c b 1 + + b⎟ ⎜ ϑ+z+F ln c ∂ ⎝ ⎠ tp t + t − tp (29) In eq 29, t is the time elapsed from the beginning of the TP experiment and extends to the relaxation period to the point where the cell retrieves its equilibrium (Δcb = 0). During the TP experiment, this is the potential across the cell and not the concentration difference which is recorded. For sufficiently small perturbations (I√tp → 0), Δcb can be evaluated from eq 30: (27) where Φ is a well-defined potential in the electrolyte with respect to a lithium reference electrode. Application of eq 27 to the results from a series of concentration-cell experiments provides one with the product ∂ ln f b ± ) and the transference of thermodynamic factor (1 + ∂ln c b 0,b number of the anion (i.e., 1 − t+ ). In the concentration-cell experiments, the potential difference is recorded (I = 0) between the two lithium electrodes (ΔΦ) in a cell filled with Δc b = 563 ΔΦc∞b ( d(ΔΦ) ) d ln c b CC,c b ∞ (30) DOI: 10.1021/acscentsci.6b00169 ACS Cent. Sci. 2016, 2, 560−568 Research Article ACS Central Science Table 1. Calculated Partial Molar Volumes of Solvent (i.e., DOL/DME) and Solutes in the Binary and Ternary Mixtures of Li2S4, Li2S6, and LiTFSIa component a solution V̅ Li2S4 (cm /mol) DOL/DME-Li2S4 DOL/DME-Li2S6 DOL/DME-LiTFSI DOL/DME-Li2S4−LiTFSI DOL/DME-Li2S6−LiTFSI 44 3 V̅ Li2S6 (cm /mol) 3 V̅ LiTFSI (cm3/mol) V̅ DOL/DME (cm3/mol) 134 128 125 84 84 83 84 83 65 38 75 solution density (g/cm3) ρ ρ ρ ρ ρ = = = = = 0.0996c13 + 0.9647 0.1438c13 + 0.9638 0.157c12 + 0.968 0.1076c13 + 0.1638c12 + 0.9644 0.1336c13 + 0.1656c12 + 0.9689 Calculations are based on the solution densities measured at 25 °C. Individual solute concentrations are between 0.1 and 1 M. Figure 3. (a) Thermodynamic factor and (b) mean-molar ionic activity coefficients of solutes in the binary electrolytes of LiTFSI-DOL/DME (dashdot-○-), Li2S4-DOL/DME (dot-△-), and Li2S6-DOL/DME (dash-*-) at T = 25 °C. Markers are experimental data and lines represent the empirical fits or processed data. which a linear trend is found toward equilibrium. The slope of the linear tail (sD) is related to DV,b ∞ according to eq 33 where the denominator on the right-hand side of the equality is the slope of a ΔΦ−ln cb plot evaluated at cb = cb∞ and is accessible through the data analysis of the concentration-cell experiments. Accordingly, combination of eqs 28 and 30 results in the following correlation between ΔΦ and τ in TP experiments ⎡⎛ 4(1 − t+0,b) d(ΔΦ) ⎞ ΔΦ = ⎢⎜ ⎟ ⎢⎣⎝ d ln c b ⎠CC, c b F ϑ+z+c∞b D∞V,b ∞ ⎤ D∞V,b I t p ⎥τ ⎥⎦ π −sD = F ϑ+z+c∞b πD∞V,b ( 4 d(ΔΦ) d ln c b ) b CC, c∞ (33) where L is the distance between the two electrodes. ■ DISCUSSION It is important to note that the application of eqs 27−33 for the case of polysulfide systems (eqs 24 and 25) is subject to two main assumptions. First, any disproportionation reactions cannot be explicitly accounted for owing to the complexity involved (and lack of experimental data on their kinetics), and second, the flux of polysulfide anions (i.e., S4−2 and S6−2) is taken to be zero at both electrodes. A more comprehensive analysis would be possible only with a priori access to the kinetic parameters of the polysulfides in the disproportionation and side reactions at the bulk of electrolyte and surface of the lithium electrode, respectively. We believe these are both reasonable assumptions. Under dilute conditions, and in a solvent of intermediate polarity such as DOL/DME, disproportionation reactions at static Li/S ratios as used here are not expected to be high.45 The macroscopic transport coefficients for the three binary electrolyte systems (i.e., eqs 23−25) at 25 °C and for solute concentrations in the range of 0.1−1 M are presented in Figure 2a−c. The diffusion coefficient of the lithium salt and lithium polysulfides features a decreasing trend upon the concentration increase, and the average DV,b is the same order of magnitude (31) The slope of ΔΦ−τ plots are calculated for a series of TP tests with different values of I√tp. The quantity inside the bracket in eq 31 is then plotted against I√tp and the resulting slope represented by m is used for unequivocal resolution of transference number: t+0,b = 1 − π 2 V,b D∞ L2 m (32) An a priori knowledge of DV,b ∞ enables us to calculate both the thermodynamic factor and the transference number at the same concentration of solute by the application of eq 27 and eq 32, respectively. Here, DV,b ∞ is obtained with the application of the method of the restricted diffusion (RD) to a similar cell used in TP experiment.44 Contrary to the TP experiment where tp is in the order of few minutes, in RD, the cell is held under polarization for a very long time (i.e., few hours). The natural logarithm of the potential relaxation is plotted against time of 564 DOI: 10.1021/acscentsci.6b00169 ACS Cent. Sci. 2016, 2, 560−568 Research Article ACS Central Science Figure 4. Lithium transference numbers at T = 25 °C, in the ternary electrolytes (a) LiTFSI-Li2S4 and (d) LiTFSI-Li2S6 and their ratio with respect to the corresponding values in a binary electrolyte of LiTFSI at the same concentration of LiTFSI as in those of ternary (c) LiTFSI-Li2S4 and (f) and LiTFSI-Li2S6. Polysulfide transference numbers in the ternary electrolytes of (b) LiTFSI-Li2S4 and (e) LiTFSI-Li2S6. (3.5 × 10−10 < (DV,b)ave < 6.5 × 10−10 m2/s) for the three solutes (Figure 2a). LiTFSI and Li2S6 have the highest and lowest DV,b over the whole concentration range, respectively. The transference number of lithium ions in LiTFSI-DOL/DME is significantly higher than those of polysulfide systems and falls in a range of 0.47−0.57. The fraction of current (in the absence of a concentration gradient) carried by Li+ in Li2S6-DOL/DME is in the range of 0.001−0.2, whereas the negative transference numbers clearly distinguish the Li2S4-DOL/DME among the three binary systems, i.e., −0.39 < t0,b + < −0.32 (Figure 2b). These small transference numbers impose a very steep gradient of solute concentrations to the polysulfide binary systems under current control and the Li2S4 species is clearly the most vulnerable. The ionic conductivity (Figure 2c) of the electrolyte in LiTFSI system is well above 1 mS/cm over the whole concentration range and increases to 14.7 mS/cm at 1 M. The average ionic conductivity in Li2S4-DOL/DME and Li2S6DOL/DME systems are less than 20% and 29% of the corresponding value in LiTFSI system, respectively. Table 1 summarizes the measured electrolyte density and partial molar volumes of the solute in the binary and ternary electrolytes. The thermodynamic factor (1 + ∂ ln f ±b ∂ ln c b Scatchard’s neutral-electrolyte description (i.e., eqs 20−22) to the binary data (i.e., Figure 3b) predicts a similar situation (Figure S1) for the ternary electrolytes of LiTFSI with Li2Sn (n = 6, 4). The presence of polysulfides in the ternary electrolytes of LiTFSI is concurrent with a substantial change in the mobility of lithium ions. Figure 4a,b and Figure 4d,e present the transference numbers of lithium ions and polysulfides in the ternary electrolytes of LiTFSI-Li2S4-DOL/DME and LiTFSILi2S6-DOL/DME, respectively. Analogous to the trends found for the binary electrolytes of lithium-polysulfides (Figure 2b), the lithium ions have a higher transference number in the ternary electrolyte of Li2S6 (Figure 4d) in comparison to the Li2S4 ternary system (Figure 4a). The lithium transference is as low as −0.2 for extremely low concentrations of LITFSI and polysulfide concentrations above 0.1 M. The contribution of polysulfide anions in the ionic current is not negligible (Figure 4b,e) even at very low concentrations of Li2Sn (n = 4, 6). A better illustration of the polysulfide’s detrimental impact on the transference of lithium ions is presented in Figure 4c,f. In these figures, the ratio of ternary and binary transference numbers of lithium ions for the same concentration of LiTFSI is presented for the ternary electrolytes based on Li2S4 (Figure 4c) and Li2S6 (Figure 4f). Introduction of the Li2Sn (n = 4, 6) into the binary electrolyte of LiTFSI is concomitant with up to 20% and 80% of reduction in the transference number of lithium ions for polysulfide concentrations as low as 0.01 and 0.1 M, respectively. At lithium salt concentrations below 0.5 M, the onset of polysulfide-induced decline in the cationic transference number of the LiTFSI binary electrolytes stands at a very low concentration of polysulfide species, e.g., 20% decrease in t0Li+ for a 0.3 M LiTFSI is observed as early as polysulfide concentration exceeds 0.005 M. The analysis of diffusion in a ternary electrolyte is not trivial since a solute diffuses not only as a response to concentration ) of the three solutes in DOL/DME are presented in Figure 3a. The thermodynamic factor increases with concentration for the three binary systems and significantly deviates from the value expected for an ideal solution (i.e., 1) at lower concentrations. The mean-molar ionic activity coefficients are available by integration of the curves in Figure 3a and with respect to the secondary state of reference. Figure 3b presents these activity coefficients for the three binary systems. The activity coefficients are far from unity for the three electrolyte cases and the highest degree of nonideality belongs to the polysulfide solutions. Higher gradients of potential in these electrolytes are further expected as a direct consequence of their extreme nonideality. Application of the 565 DOI: 10.1021/acscentsci.6b00169 ACS Cent. Sci. 2016, 2, 560−568 Research Article ACS Central Science Figure 5. Practical diffusion coefficients of LiTFSI (D022) at T = 25 °C in ternary electrolytes of (a) LiTFSI-Li2S4 and (c) LiTFSI-Li2S6 and those of polysulfides (D033) in (b) LiTFSI-Li2S4 and (d) LiTFSI-Li2S6. 0 0 of D23 and D32 presented in Figure S2 of Supporting Information. Very high values of D032 (Figure S2b,d) indicate that a subtle gradient in the concentration of LiTFSI has a significant impact on the diffusion of polysulfide species (i.e., both Li2S4 and Li2S6). Accordingly, higher rates of cell discharge and charge are counterintuitively in favor of polysulfide retention and crossover, respectively. This might suggest that the relaxation periods after slow rate of discharge and high rate of charge are detrimental to the battery state of health. In contrast to the other two transport properties and not surprisingly, the ionic conductivity of the electrolyte improves upon polysulfide build-up: especially for lower concentrations of LiTFSI (Figure S3). While it is risky to speculate on the underlying difference in the behavior of the two polysulfides that this analysis clearly points to, we note that under the relatively dilute conditions employed here (<1 M polysulfide) we can expect disproportionation to occur, although it would be less compared to more concentrated solutions. Thus, S62− is probably less prone to disproportionation in DOL/DMEgenerally being subject to cleavage to S3−• only in EDP solvents25whereas Li2S4 is likely to be the end polysulfide that precedes precipitation of Li2S via disproportionation, as it represents the polysulfide species at the inflection point in a theoretical discharge profile.20,24−26,45 Computational studies may shed light on the complex dynamics of this system. Further analysis of the system in electron-pair donor (EPD) solvents such as DMSO and DMA where the trisulfide ion dominates the kinetics would be expected to yield very different results, and will be the subject of future studies. gradient of its own, but also to that of the second solute (eqs 3−5). The practical coefficients characterizing the diffusion rate of LiTFSI and Li2Sn (n = 4, 6) under their own concentration gradient, i.e., D022 and D033, respectively, are presented in Figure 5. LiTFSI has slightly higher diffusion coefficients (D022) in the presence of Li2S4 (Figure 5a) rather than Li2S6 (Figure 5c). The 0 diffusion behavior (D33 ) of Li2S4 and Li2S6 in ternary electrolytes of LiTFSI is significantly different. In LiTFSILi2S6-DOL/DME ternary electrolyte, the diffusion coefficient of Li2S6 (D033 in Figure 5d) is superior to that of LiTFSI (D022 in Figure 5c), up to 1 order of magnitude. However, the incongruent diffusion of Li2S4 (D033 in Figure 5b) is dominant for polysulfide concentrations lower than 0.01 M and higher than 0.1 M over the full range and cLiTFSI < 0.3 M, respectively. ■ CONCLUSIONS Providing D033 < 0, i.e., incongruent diffusion, the Li2S4 species diffuse against their own concentration gradient. This phenomenon is in favor of polysulfide retention at the cathode and against polysulfide crossover toward the lithium anode during cell discharge. The above-mentioned D033 behavior suggests that the kinetics of Li2S6 to Li2S4 reduction has a pivotal role in the onset of polysulfide crossover phenomena. A longer residence time of Li2S6 at the cathode during cell discharge, i.e., sluggish kinetics of Li2S6 reduction, is readily followed by the fast crossover of these species toward the anode. On the contrary, under the fast kinetics of Li2S6 reduction, Li2S4 species have a very low tendency to shuttle toward the anode during cell discharge. Our diffusion analysis reveals the significant share of Li2S6 in the shuttle mechanism and further aging of the lithium−sulfur batteries with a similar electrolyte formulation (i.e., LiTFSI-DOL/DME). The importance of cross-term diffusion coefficients cannot be overlooked. This is clearly evident from the magnitude orders ■ EXPERIMENTAL METHODS All experiments were conducted at 25 °C inside an argon-filled glovebox with oxygen and water contents below 3 ppm. A VMP3 566 DOI: 10.1021/acscentsci.6b00169 ACS Cent. Sci. 2016, 2, 560−568 ACS Central Science ■ potentiostat was used for the AC-impedance measurements and the galvanostatic polarization studies, together with relaxation tests performed with the aid of a Macpile potentiostat. The density of the electrolytes was measured using an oscillating U-tube digital density meter (DMA 35, Anton-Paar). To prepare the polysulfides, precise stoichiometric amounts of lithium (as the stabilized powder; FMC Corporation), sulfur, and/or LiTFSI together with the solvent were mixed inside a magnetically stirred glass vial for 4 days. The electrolytes were filtered to remove any leftover reactants. The 1.27 cm2 circular lithium electrodes were punched out from a lithium foil with a thickness of 150 μm. All the cells were subjected to three formation cycles after assembly in order to partially stabilize the surface of lithium electrodes. A formation cycle consisted of a charge and discharge segment at 0.05 mA and continued for 0.5 h in each segment. The galvanostatic polarization pulses were 0.1−2 mA in magnitude and continued for 20 min and 10 h in the experiments involved in the measure of the transference number and diffusion coefficient, respectively. Relaxation was recorded until equilibrium was reached both during concentration cell studies and after galvanostatic polarization periods in TP and RD. The AC-impedance spectra were obtained over a frequency range of 200 kHz−50 mHz and with a potential amplitude of 15 mV. The ionic conductivity of three ternary samples (i.e., (0.2 M LiTFSI, 0.2 M Li2S6), (0.4 M TFSI, 0.1 M Li2S6), (0.3 M TFSI, 0.2 M Li2S6)) were measured in order to check the accuracy of the prediction methodology proposed in this work. The predictions overestimate the experimental data with a maximum error of 20%, which is reasonable considering the experimental protocol. A representative set of data from the restricted diffusion, transference polarization, and concentration-cell experiments are included in the Supporting Information for the Li2S6 electrolytes. Figures S4 and S5 present the data from the TP and RD experiments, respectively, for the binary electrolyte of Li2S6 at 0.5M. The results of the concentration-cell experiments are summarized in Figure S6, for the concentration cells at five different concentrations of Li2S6, i.e., 0.1, 0.3, 0.5, 0.7, and 1 M. ■ SYMBOLS AND NOMENCLATURE ρ density c concentration + , D diffusion coefficient I binary/ternary fundamental coefficients L distance between two electrodes Φ potential vs Li t transference number and time ϑ stochiometric coefficient F Faraday constant z charge number π Pi number m slope in transference-polarization experiment tp pulse duration τ dimensionless time I geometrical-surface current density κ ionic conductivity of the electrolyte y ionic-strength fraction of the solute in a ternary electrolyte M molar mass N number of equivalents f± mean-molar ionic activity coefficient V partial molar volume μ electrochemical potential Superscripts 0 solvent-fixed frame (SF) b binary electrolyte V volume-fixed frame (VF) Subscripts a anion c cation ∞ electrolyte bulk at infinite time or distance from the electrode + cation in the binary electrolyte − anion in the binary electrolyte 0 solvent (i.e., DOL/DME) 1 cation in the ternary electrolyte with common cation (Li+) 2 anion of the primary solute in the ternary electrolyte with common cation (i.e., TFSI−) 3 anion of the secondary solute in the ternary electrolyte with common cation (i.e., S6−2/S4−2) 12 first solute in the ternary electrolyte with common cation (i.e., LiTFSI) 13 second solute in the ternary electrolyte with common cation (i.e., Li2S6/Li2S4) ASSOCIATED CONTENT S Supporting Information * The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acscentsci.6b00169. Experimental results and estimations for activity coefficients, cross-term diffusion coefficients, and ionic conductivity of the ternary electrolytes, together with a representative set of raw data for the measurement of binary coefficients (PDF) ■ Research Article ■ AUTHOR INFORMATION Corresponding Author *E-mail: [email protected]. REFERENCES (1) Evers, S.; Nazar, L. F. New approaches for high energy density lithium-sulfur battery cathodes. Acc. Chem. Res. 2013, 46, 1135−1141. (2) Manthiram, A.; Fu, Y.; Chung, S. H.; Zu, C.; Su, Y. S. Rechargeable lithium-sulfur batteries. Chem. Rev. 2014, 114, 11751− 11787. (3) Erickson, E. M.; Markevich, E. M.; Salitra, G.; Sharon, D.; Hirshberg, D.; de la Llave, E.; Shterenberg, I.; Rozenman, A.; Frimer, A.; Aurbach, D. Development of advanced rechargeable batteries: a continuous challenge in the choice of suitable electrolyte solution. J. Electrochem. Soc. 2015, 162, A2424−A2438. (4) Lv, D.; Zheng, J.; Li, Q.; Xie, X.; Ferrara, S.; Nie, Z.; Mehdi, L. B.; Browning, N. D.; Zhang, J. G.; Graff, G. L.; Liu, J.; Xiao, J. High energy density lithium-sulfur batteries: challenges of thick sulfur cathodes. Adv. Energy Mater. 2015, 5, 1402290−1402298. (5) Park, K.; Cho, J. H.; Jang, J. H.; Yu, B. C.; de la Hoz, A. T.; Miller, K. M.; Ellison, C. J.; Goodenough, J. B. Trapping lithium polysulfides Present Address † (M.S.) Department of Industrial Engineering, University of Hasselt, Hasselt, Belgium. Notes The authors declare no competing financial interest. ■ ACKNOWLEDGMENTS We gratefully acknowledge the BASF International Scientific Network for Electrochemistry and Batteries for funding. L.F.N. thanks NSERC for a Canada Research Chair, and support via their Discovery Grants program. CYK also thanks the Government of Ontario for Ontario Graduate Scholarship. 567 DOI: 10.1021/acscentsci.6b00169 ACS Cent. Sci. 2016, 2, 560−568 Research Article ACS Central Science of a Li-S battery by forming lithium bonds in a polymer matrix. Energy Environ. Sci. 2015, 8, 2389−2395. (6) Li, W.; Yao, H.; Yan, K.; Zheng, G.; Liang, Z.; Chiang, Y. M.; Cui, Y. The synergetic effect of lithium polysulfide and lithium nitrate to prevent lithium dendrite growth. Nat. Commun. 2015, 6, 7436−7443. (7) Mikhaylik, Y. V.; Akridge, J. R. Polysulfide shuttle study in the Li/ S battery system. J. Electrochem. Soc. 2004, 151, A1969−A1976. (8) Wild, M.; O’Neill, L. O.; Zhang, T.; Purkayastha, R.; Minton, G.; Marinescu, M.; Offer, G. Lithium sulfur batteries, a mechanistic review. Energy Environ. Sci. 2015, 8, 3477−3494. (9) Zhang, S. S.; Read, J. A. A new direction for the performance improvement of rechargeable lithium/sulfur batteries. J. Power Sources 2012, 200, 77−82. (10) Manthiram, A.; Chung, S.-H.; Zu, C. Lithium-Sulfur batteries: progress and prospects. Adv. Mater. 2015, 27, 1980−2006. (11) Pang, Q.; Liang, X.; Kwok, C. Y.; Nazar, L. F. The importance of chemical interactions between sulfur host materials and lithium polysulfides for advanced lithium-sulfur batteries. J. Electrochem. Soc. 2015, 162, A2567−A2576. (12) He, G.; Evers, S.; Liang, X.; Cuisinier, M.; Garsuch, A.; Nazar, L. F. Tailoring porosity in carbon nanospheres for lithium-sulfur battery cathodes. ACS Nano 2013, 7, 10920−10930. (13) Elazari, R.; Salitra, G.; Garsuch, A.; Panchenko, A.; Aurbach, D. Sulfur-impregnated activated carbon fiber cloth as a binder-free cathode for rechargeable Li-S batteries. Adv. Mater. 2011, 23, 5641− 5644. (14) Liang, X.; Garsuch, A.; Nazar, L. F. Sulfur cathodes based on conductive MXene nanosheets for high-performance lithium-sulfur batteries. Angew. Chem., Int. Ed. 2015, 54, 3907−3911. (15) Liang, X.; Kwok, C. Y.; Lodi-Marzano, L. M.; Pang, Q.; Cuisinier, M.; Huang, H.; Hart, C. J.; Houtarde, D.; Kaup, K.; Sommer, H.; Brezesinski, T.; Janek, J.; Nazar, L. F. Tuning transition metal oxide-sulfur interactions for long life lithium sulfur batteries: the ‘Goldilock’ principle. Adv. Energy Mater. 2016, 6, 150163610.1002/ aenm.201501636. (16) Barghamadi, M.; Best, A. S.; Bhatt, A. I.; Hollenkamp, A. F.; Musameh, M.; Rees, R. J.; Rüther, T. Lithium-sulfur batteries - the solution is in the electrolyte, but is the electrolyte a solution? Energy Environ. Sci. 2014, 7, 3902−3920. (17) Xu, R.; Lu, J.; Amine, K. Progress in mechanistic understanding and characterization techniques of Li-S batteries. Adv. Energy Mater. 2015, 5, 150040810.1002/aenm.201500408. (18) Lu, Y.-C.; He, Q.; Gasteiger, H. A. Probing the lithium-sulfur redox reactions: A rotating-ring disk electrode study. J. Phys. Chem. C 2014, 118, 5733−5741. (19) Barchasz, C.; Molton, F.; Duboc, C.; Lepretre, J. C.; Patoux, S.; Alloin, F. Lithium/sulfur cell discharge mechanism: An original approach for intermediate species identification. Anal. Chem. 2012, 84, 3973−3980. (20) Cuisinier, M.; Cabelguen, P.-E.; Evers, S.; He, G.; Kolbeck, M.; Garsuch, A.; Bolin, T.; Balasubramanian, M.; Nazar, L. F. Sulfur speciation in Li-S batteries determined by operando X-ray absorption spectroscopy. J. Phys. Chem. Lett. 2013, 4, 3227−3232. (21) Gorlin, Y.; Siebel, A.; Piana, M.; Huthwelker, T.; Jha, H.; Monsch, G.; Kraus, F.; Gasteiger, H. A.; Tromp, M. Operando characterization of intermediates produced in a lithium-sulfur battery. J. Electrochem. Soc. 2015, 162, A1146−A1155. (22) Chivers, T. Ubiquitous trisulfur radical ion S3−•. Nature 1974, 252, 32−33. (23) Chivers, T.; Elder, P. J. W. Ubiquitous trisulfur radical anion: fundamentals and applications in materials science, electrochemistry, analytical chemistry and geochemistry. Chem. Soc. Rev. 2013, 42, 5996−6005. (24) Cuisinier, M.; Hart, C.; Balasubramanian, M.; Garsuch, A.; Nazar, L. F. Radical or not radical: revisiting lithium-sulfur electrochemistry in non-aqueous electrolytes. Adv. Energy Mater. 2015, 5, 1401801. (25) Vijayakumar, M.; Govind, N.; Walter, E.; Burton, S. D.; Shukla, A.; Devaraj, A.; Xiao, J.; Liu, J.; Wang, C.; Karim, A.; Thevuthasan, S. Molecular structure and stability of dissolved lithium polysulfide species. Phys. Chem. Chem. Phys. 2014, 16, 10923−10932. (26) Kumaresan, K.; Mikhaylik, Y.; White, R. E. A mathematical model for a lithium-sulfur cell. J. Electrochem. Soc. 2008, 155, A576− A582. (27) Hofmann, A. F.; Fronczek, D. N.; Bessler, W. G. Mechanistic modeling of polysulfide shuttle and capacity loss in lithium-sulfur batteries. J. Power Sources 2014, 259, 300−310. (28) Albertus, P.; Christensen, J.; Newman, J. Experiments on and modeling of positive electrodes with multiple active materials for lithium-ion batteries. J. Electrochem. Soc. 2009, 156, A606−A618. (29) Christensen, J.; Newman, J. A mathematical model for the lithium-ion negative electrode solid electrolyte interphase. J. Electrochem. Soc. 2004, 151, A1977−A1988. (30) Christensen, J.; Newman, J. Stress generation and fracture in lithium insertion materials. J. Solid State Electrochem. 2006, 10, 293− 319. (31) Christensen, J.; Newman, J. Cyclable lithium and capacity loss in Li-ion cells. J. Electrochem. Soc. 2005, 152, A818−A829. (32) Zhang, T.; Marinescu, M.; O’Neill, L.; Wild, M.; Offer, G. Modeling the voltage loss mechanisms in lithium-sulfur cells: the importance of electrolyte resistance and precipitation kinetics. Phys. Chem. Chem. Phys. 2015, 17, 22581−22586. (33) Marinescu, M.; Zhang, T.; Offer, G. J. A zero dimensional model of lithium-sulfur batteries during charge and discharge. Phys. Chem. Chem. Phys. 2016, 18, 584−593. (34) Eroglu, D.; Zavadil, K. R.; Gallagher, K. G. Critical link between materials chemistry and cell-level design for high energy density and low cost lithium-sulfur transportation battery. J. Electrochem. Soc. 2015, 162, A982−A990. (35) Ma, Y.; Doyle, M.; Fuller, T. F.; Doef, M. M.; de Jonghe, L. C.; Newman, J. Measurement of a complete set of transport properties for a concentrated solid polymer electrolyte solution. J. Electrochem. Soc. 1995, 142, 1859−1868. (36) de Groot, S. R.; Mazur, P. Non-Equilibrium Thermodynamics; Dover Publications; New York, 1984. (37) Miller, D. G. Application of irreversible thermodynamics to electrolyte solutions. I. determination of ionic transport coefficients lij for isothermal vector transport processes in binary electrolyte systems. J. Phys. Chem. 1966, 70, 2639−2659. (38) Miller, D. G. Application of irreversible thermodynamics to electrolyte solutions. II. Ionic coefficients lij for isothermal vector transport processes in ternary systems. J. Phys. Chem. 1967, 71, 616− 632. (39) Guggenheim, E. A. Thermodynamics, An Advanced Treatment for Chemists and Physicists; Elsevier Science Publishers B. V.; Amsterdam North-Holland, 1967. (40) Kirkwood, J. G.; Baldwin, R. L.; Dunlop, P. J.; Gosting, L. J.; Kegeles, G. Flow equations and frame of reference for isothermal diffusion in liquids. J. Chem. Phys. 1960, 33, 1505−1513. (41) Newman, J.; Thomas-Alyea, K. E. Electrochemical Systems; Wiley Interscience; New Jersey, USA, 2004. (42) Miller, D. G. Activity coefficient derivatives of ternary systems based on Scatchard’s neutral electrolyte description. J. Solution Chem. 2008, 37, 365−375. (43) Hafezi, H.; Newman, J. Verification and analysis of transference number measurements by the galvanostatic polarization method. J. Electrochem. Soc. 2000, 147, 3036−3042. (44) Newman, J.; Chapman, T. W. Restricted diffusion in binary solutions. AIChE J. 1973, 19, 343−348. (45) Cuisinier, M.; Cabelguen, P. E.; Adams, B. D.; Garsuch, A.; Balasubramanian, M.; Nazar, L. F. Unique behavior of nonsolvents for polysulfides in lithium-sulfur batteries. Energy Environ. Sci. 2014, 7, 2697−2705. 568 DOI: 10.1021/acscentsci.6b00169 ACS Cent. Sci. 2016, 2, 560−568