Ion Transport in Polymer Electrolytes: Building New Bridges between Experiment and Molecular Simulation.

ConspectusPolymer electrolytes constitute a promising type of material for solid-state batteries. However, one of the bottlenecks for their practical implementation lies in the transport properties, often including restricted Li+ self-diffusion and conductivity and low cationic transference numbers. This calls for a molecular understanding of ion transport in polymer electrolytes in which molecular dynamics (MD) simulation can provide both new physical insights and quantitative predictions. Although efforts have been made in this area and qualitative pictures have emerged, direct and quantitative comparisons between experiment and simulation remain challenging because of the lack of a unified theoretical framework to connect them.In our work, we show that by computing the glass transition temperature (Tg) of the model system and using the normalized inverse temperature 1000/(T - Tg + 50), the Li+ self-diffusion coefficient can be compared quantitatively between MD simulations and experiments. This allows us to disentangle the effects of Tg and the polymer dielectric environment on ion conduction in polymer electrolytes, giving rise to the identification of an optimal solvating environment for fast ion conduction.Unlike Li+ self-diffusion coefficients and ionic conductivity, the transference number, which describes the fraction of current carried by Li+ ions, depends on the boundary conditions or the reference frame (RF). This creates a non-negligible gap when comparing experiment and simulation because the fluxes in the experimental measurements and in the linear response theory used in MD simulation are defined in different RFs. We show that by employing the Onsager theory of ion transport and applying a proper RF transformation, a much better agreement between experiment and simulation can be achieved for the PEO-LiTFSI system. This further allows us to derive the theoretical expression for the Bruce-Vincent transference number in terms of the Onsager coefficients and make a direct comparison to experiments. Since the Bruce-Vincent method is widely used to extract transference numbers from experimental data, this opens the door to calibrating MD simulations via reproducing the Bruce-Vincent transference number and using MD simulations to predict the true transference number.In addition, we also address several open questions here such as the time-scale effects on the ion-pairing phenomenon, the consistency check between different types of experiments, the need for more accurate force fields used in MD simulations, and the extension to multicomponent systems. Overall, this Account focuses on building new bridges between experiment and simulation for quantitative comparison, warnings of pitfalls when comparing apples and oranges, and clarifying misconceptions. From a physical chemistry point of view, it connects to concentrated solution theory and provides a unified theoretical framework that can maximize the power of MD simulations. Therefore, this Account will be useful for the electrochemical energy storage community at large and set examples of how to approach experiments from theory and simulation (and vice versa).

The sensitive aspects of modelling polymer-ceramic composite solid-state electrolytes using molecular dynamics simulations.

Solid-state composite electrolytes have arisen as one of the most promising materials classes for next-generation Li-ion battery technology. These composites mix ceramic and solid-polymer ion conductors with the aim of combining the advantages of each material. The ion-transport mechanisms within such materials, however, remain elusive. This knowledge gap can to a large part be attributed to difficulties in studying processes at the ceramic-polymer interface, which are expected to play a major role in the overall ion transport through the electrolyte. Computational efforts have the potential of providing significant insight into these processes. One of the main challenges to overcome is then to understand how a sufficiently robust model can be constructed in order to provide reliable results. To this end, a series of molecular dynamics simulations are here carried out with a variation of certain structural (surface termination and polymer length) and pair potential (van der Waals parameters and partial charges) models of the Li7La3Zr2O12 (LLZO) poly(ethylene oxide) (PEO) system, in order to test how sensitive the outcome is to each variation. The study shows that the static and dynamic properties of Li-ion are significantly affected by van der Waals parameters as well as the surface terminations, while the thickness of the interfacial region - where the structure-dynamic properties are different as compared to the bulk-like regime - is the same irrespective of the simulation setup.

Transference Number in Polymer Electrolytes: Mind the Reference-Frame Gap.

The transport coefficients, in particular the transference number, of electrolyte solutions are important design parameters for electrochemical energy storage devices. The recent observation of negative transference numbers in PEO-LiTFSI under certain conditions has generated much discussion about its molecular origins, by both experimental and theoretical means. However, one overlooked factor in these efforts is the importance of the reference frame (RF). This creates a non-negligible gap when comparing experiment and simulation because the fluxes in the experimental measurements of transport coefficients and in the linear response theory used in the molecular dynamics simulation are defined in different RFs. In this work, we show that, by applying a proper RF transformation, a much improved agreement between experimental and simulation results can be achieved. Moreover, it is revealed that the anion mass and the anion-anion correlation, rather than ion aggregates, play a crucial role for the reported negative transference numbers.

Improved Battery Cycle Life Prediction Using a Hybrid Data-Driven Model Incorporating Linear Support Vector Regression and Gaussian.

The ability to accurately predict lithium-ion battery life-time already at an early stage of battery usage is critical for ensuring safe operation, accelerating technology development, and enabling battery second-life applications. Many models are unable to effectively predict battery life-time at early cycles due to the complex and nonlinear degrading behavior of lithium-ion batteries. In this study, two hybrid data-driven models, incorporating a traditional linear support vector regression (LSVR) and a Gaussian process regression (GPR), were developed to estimate battery life-time at an early stage, before more severe capacity fading, utilizing a data set of 124 battery cells with lifetimes ranging from 150 to 2300 cycles. Two type of hybrid models, here denoted as A and B, were proposed. For each of the models, we achieved 1.1 % (A) and 1.4 % (B) training error, and similarly, 8.3 % (A) and 8.2 % (B) test error. The two key advantages are that the error percentage is kept below 10 % and that very low error values for the training and test sets were observed when utilizing data from only the first 100 cycles.The proposed method thus appears highly promising for predicting battery life during early cycles.

Understanding the lithiation limits of high-capacity organic battery anodes by atomic charge derivative analysis.

The superlithiation of organic anodes is a promising approach for developing the next generation of sustainable Li-ion batteries with high capacity. However, the lack of fundamental understanding hinders its faster development. Here, a systematic study of the lithiation processes in a set of dicarboxylate-based materials is carried out within the density functional theory formalism. It is demonstrated that a combined analysis of the Li insertion reaction thermodynamics and the conjugated-moiety charge derivative enables establishing the experimentally observed maximum storage, thus allowing an assessment of the structure-function relationships also.

A molecular dynamics study of a fully zwitterionic copolymer/ionic liquid-based electrolyte: Li+ transport mechanisms and ionic interactions.

The development of polymer electrolytes (PEs) is crucial for advancing safe, high-energy density batteries, such as lithium-metal and other beyond lithium-ion chemistries. However, reaching the optimum balance between mechanical stiffness and ionic conductivity is not a straightforward task. Zwitterionic (ZI) gel electrolytes comprising lithium salt and ionic liquid (IL) solutions within a fully ZI polymer network can, in this context, provide useful properties. Although such materials have shown compatibility with lithium metal in batteries, several fundamental structure-dynamic relationships regarding ionic transport and the Li+ coordination environment remain unclear. To better resolve such issues, molecular dynamics simulations were carried out for two IL-based electrolyte systems, N-butyl-N-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide ([BMP][TFSI]) with 1 M LiTFSI salt and a ZI gel electrolyte containing the IL and a ZI copolymer: poly(2-methacryloyloxyethyl phosphorylcholine-co-sulfobetaine vinylimidazole), poly(MPC-co-SBVI). The addition of ZI polymer decreases the [TFSI]- -[Li]+ interactions and increases the IL ion diffusivities, and consequently, the overall ZI gel ionic conductivity. The structural analyses showed a large preference for lithium-ion interactions with the polymer phosphonate groups, while the [TFSI]- anions interact directly with the sulfonate group and the [BMP]+ cations only display secondary interactions with the polymer. In contrast to previous experimental data on the same system, the simulated transference numbers showed smaller [Li]+ contributions to the overall ionic conductivities, mainly due to negatively charged lithium aggregates and the strong lithium-ion interactions in the systems.

Boosting Rechargeable Batteries R&D by Multiscale Modeling: Myth or Reality?

This review addresses concepts, approaches, tools, and outcomes of multiscale modeling used to design and optimize the current and next generation rechargeable battery cells. Different kinds of multiscale models are discussed and demystified with a particular emphasis on methodological aspects. The outcome is compared both to results of other modeling strategies as well as to the vast pool of experimental data available. Finally, the main challenges remaining and future developments are discussed.

The role of coordination strength in solid polymer electrolytes: compositional dependence of transference numbers in the poly(ε-caprolactone)-poly(trimethylene carbonate) system.

Both polyesters and polycarbonates have been proposed as alternatives to polyethers as host materials for future polymer electrolytes for solid-state lithium-ion batteries. While being comparatively similar functional groups, the electron density on the coordinating carbonyl oxygen is different, thereby rendering different coordinating strength towards lithium ions. In this study, the transport properties of poly(ε-caprolactone) and poly(trimethylene carbonate) as well as random copolymers of systematically varied composition of the two have been investigated, in order to better elucidate the role of the coordination strength. The cationic transference number, a property well-connected with the complexing ability of the polymer, was shown to depend almost linearly on the ester content of the copolymer, increasing from 0.49 for the pure poly(ε-caprolactone) to 0.83 for pure poly(trimethylene carbonate). Contradictory to the transference number measurements that suggest a stronger lithium-to-ester coordination, DFT calculations showed that the carbonyl oxygen in the carbonate coordinates more strongly to the lithium ion than that of the ester. FT-IR measurements showed the coordination number to be higher in the polyester system, resulting in a higher total coordination strength and thereby resolving the paradox. This likely originates in properties that are specific of polymeric solvent systems, e.g. steric properties and chain dynamics, which influence the coordination chemistry. These results highlight the complexity in polymeric systems and their ion transport properties in comparison to low-molecular-weight analogues, and how polymer structure and steric effects together affect the coordination strength and transport properties.

Simultaneous Monitoring of Crystalline Active Materials and Resistance Evolution in Lithium-Sulfur Batteries.

Operando X-ray diffraction (XRD) is a valuable tool for studying secondary battery materials as it allows for the direct correlation of electrochemical behavior with structural changes of crystalline active materials. This is especially true for the lithium-sulfur chemistry, in which energy storage capability depends on the complex growth and dissolution kinetics of lithium sulfide (Li2S) and sulfur (S8) during discharge and charge, respectively. In this work, we present a novel development of this method through combining operando XRD with simultaneous and continuous resistance measurement using an intermittent current interruption (ICI) method. We show that a coefficient of diffusion resistance, which reflects the transport properties in the sulfur/carbon composite electrode, can be determined from analysis of each current interruption. Its relationship to the established Warburg impedance model is validated theoretically and experimentally. We also demonstrate for an optimized electrode formulation and cell construction that the diffusion resistance increases sharply at the discharge end point, which is consistent with the blocking of pores in the carbon host matrix. The combination of XRD with ICI allows for a direct correlation of structural changes with not only electrochemical properties but also energy loss processes at a nonequilibrium state and, therefore, is a valuable technique for the study of a wide range of energy storage chemistries.

Interactions and Transport in Highly Concentrated LiTFSI-based Electrolytes.

To elucidate what properties control and practically limit ion transport in highly concentrated electrolytes (HCEs), the viscosity, ionic conductivity, ionicity, and transport numbers were studied for nine model electrolytes and connected to the rate capability in Li-ion battery (LIB) cells. The electrolytes employed the LiTFSI salt in three molar ratio concentrations; 1 : 2, 1 : 4, and 1 : 16 (LiTFSI:X) vs. solvents (X) with different permittivities; tert-butyl methyl ether (MTBE), tetrahydrofuran (THF) and propylene carbonate (PC). While the low polarity MTBE creates liquid electrolytes, ion-pairing limits the ionic conductivity despite extremely low viscosities. For the less concentrated 1 : 16 LiTFSI:MTBE and 1 : 16 LiTFSI:THF electrolytes the ionic diffusivities decrease with increased temperature, a sign of aggregation, but still their ionic conductivities and LIB performance increase. In general, the low ionic conductivity and high viscosity both limit the use of HCEs in LIBs, and no compensating mechanism seems to be present.

Electronic conductivity of polymer electrolytes: electronic charge transport properties of LiTFSI-doped PEO.

The electronic structure of poly(ethyleneoxide) with and without a common electrolyte lithium bis(trifluoromethane)sulfonimide salt is calculated from first principles. Introducing the salt into the polymer electrolyte significantly reduces the band gap, down to 0.6 eV. Thus, this will have a significant impact on the leakage currents in polymer electrolytes used in all-solid-state batteries.

Ion-Conductive and Thermal Properties of a Synergistic Poly(ethylene carbonate)/Poly(trimethylene carbonate) Blend Electrolyte.

Electrolytes comprising poly(ethylene carbonate) (PEC)/poly(trimethylene carbonate) (PTM C) with lithium bis(trifluoromethane sulfonyl)imide (LiTFSI) are prepared by a simple solvent casting method. Although PEC and PTMC have similar chemical structures, they are immiscible and two glass transitions are present in the differential scanning calorimetry (DSC) measurements. Interestingly, these two polymers change to miscible blends with the addition of LiTFSI, and the ionic conductivity increases with increasing lithium salt concentration. The optimum composition of the blend electrolyte is achieved at PEC6 PTMC4 , with a conductivity as high as 10-6 S cm-1 at 50 °C. This value is greater than that for single PEC- and PTMC-based electrolytes. Moreover, the thermal stability of the blend-based electrolytes is improved as compared to PEC-based electrolytes. It is clear that the interaction between CO groups and Li+ gives rise to a compatible amorphous phase of PEC and PTMC.

Ion transport in polycarbonate based solid polymer electrolytes: experimental and computational investigations.

Among the alternative host materials for solid polymer electrolytes (SPEs), polycarbonates have recently shown promising functionality in all-solid-state lithium batteries from ambient to elevated temperatures. While the computational and experimental investigations of ion conduction in conventional polyethers have been extensive, the ion transport in polycarbonates has been much less studied. The present work investigates the ionic transport behavior in SPEs based on poly(trimethylene carbonate) (PTMC) and its co-polymer with ε-caprolactone (CL) via both experimental and computational approaches. FTIR spectra indicated a preferential local coordination between Li(+) and ester carbonyl oxygen atoms in the P(TMC20CL80) co-polymer SPE. Diffusion NMR revealed that the co-polymer SPE also displays higher ion mobilities than PTMC. For both systems, locally oriented polymer domains, a few hundred nanometers in size and with limited connections between them, were inferred from the NMR spin relaxation and diffusion data. Potentiostatic polarization experiments revealed notably higher cationic transference numbers in the polycarbonate based SPEs as compared to conventional polyether based SPEs. In addition, MD simulations provided atomic-scale insight into the structure-dynamics properties, including confirmation of a preferential Li(+)-carbonyl oxygen atom coordination, with a preference in coordination to the ester based monomers. A coupling of the Li-ion dynamics to the polymer chain dynamics was indicated by both simulations and experiments.

Electrolyte decomposition on Li-metal surfaces from first-principles theory.

An important feature in Li batteries is the formation of a solid electrolyte interphase (SEI) on the surface of the anode. This film can have a profound effect on the stability and the performance of the device. In this work, we have employed density functional theory combined with implicit solvation models to study the inner layer of SEI formation from the reduction of common organic carbonate electrolyte solvents (ethylene carbonate, propylene carbonate, dimethyl carbonate, and diethyl carbonate) on a Li metal anode surface. Their stability and electronic structure on the Li surface have been investigated. It is found that the CO producing route is energetically more favorable for ethylene and propylene carbonate decomposition. For the two linear solvents, dimethyl and diethyl carbonates, no significant differences are observed between the two considered reduction pathways. Bader charge analyses indicate that 2 e- reductions take place in the decomposition of all studied solvents. The density of states calculations demonstrate correlations between the degrees of hybridization between the oxygen of adsorbed solvents and the upper Li atoms on the surface with the trend of the solvent adsorption energies.

Polymer-ionic liquid ternary systems for Li-battery electrolytes: Molecular dynamics studies of LiTFSI in a EMIm-TFSI and PEO blend.

This paper presents atomistic molecular dynamics simulation studies of lithium bis(trifluoromethane)sulfonylimide (LiTFSI) in a blend of 1-ethyl-3-methylimidazolium (EMIm)-TFSI and poly(ethylene oxide) (PEO), which is a promising electrolyte material for Li- and Li-ion batteries. Simulations of 100 ns were performed for temperatures between 303 K and 423 K, for a Li:ether oxygen ratio of 1:16, and for PEO chains with 26 EO repeating units. Li(+) coordination and transportation were studied in the ternary electrolyte system, i.e., PEO16LiTFSI⋅1.0 EMImTFSI, by applying three different force field models and are here compared to relevant simulation and experimental data. The force fields generated significantly different results, where a scaled charge model displayed the most reasonable comparisons with previous work and overall consistency. It is generally seen that the Li cations are primarily coordinated to polymer chains and less coupled to TFSI anion. The addition of EMImTFSI in the electrolyte system enhances Li diffusion, associated to the enhanced TFSI dynamics observed when increasing the overall TFSI anion concentration in the polymer matrix.

On the Elusive Crystallography of Lithium-Rich Layered Oxides: Novel Structural Models.

Lithium-rich layered oxides (LRLOs) are one of the most attractive families among future positive electrode materials for the so-called fourth generation of lithium-ion batteries (LIBs). Their electrochemical performance is enabled by the unique ambiguous crystal structure that is still not well understood despite decades of research. In the literature, a clear structural model able to describe their crystallographic features is missing thereby hindering a clear rationalization of the interplay between synthesis, structure, and functional properties. Here, the structure of a specific LRLO, Li1.28 Mn0.54 Ni0.13 Co0.02 Al0.03 O2 , using synchrotron X-ray diffraction (XRD), neutron diffraction (ND), and High-Resolution Transmission Electron Microscopy (HR-TEM), is analyzed. A systematic approach is applied to model diffraction patterns of Li1.28 Mn0.54 Ni0.13 Co0.02 Al0.03 O2 by using the Rietveld refinement method considering the R 3 ¯ $\bar{3}$ m and C2/m unit cells as the prototype structures. Here, the relative ability of a variety of structural models is compared to match the experimental diffraction pattern evaluating the impact of defects and supercells derived from the R 3 ¯ $\bar{3}$ m structure. To summarize, two possible models able to reconcile the description of experimental data are proposed here for the structure of Li1.28 Mn0.54 Ni0.13 Co0.02 Al0.03 O2 : namely a monoclinic C2/m defective lattice (prototype Li2 MnO3 ) and a monoclinic defective supercell derived from the rhombohedral R 3 ¯ $\bar{3}$ m unit cell (prototype LiCoO2 ).

Initial SEI formation in LiBOB-, LiDFOB- and LiBF4-containing PEO electrolytes.

A limiting factor for solid polymer electrolyte (SPE)-based Li-batteries is the functionality of the electrolyte decomposition layer that is spontaneously formed at the Li metal anode. A deeper understanding of this layer will facilitate its improvement. This study investigates three SPEs - polyethylene oxide:lithium tetrafluoroborate (PEO:LiBF4), polyethylene oxide:lithium bis(oxalate)borate (PEO:LiBOB), and polyethylene oxide:lithium difluoro(oxalato)borate (PEO:LiDFOB) - using a combination of electrochemical impedance spectroscopy (EIS), galvanostatic cycling, in situ Li deposition photoelectron spectroscopy (PES), and ab initio molecular dynamics (AIMD) simulations. Through this combination, the cell performance of PEO:LiDFOB can be connected to the initial SPE decomposition at the anode interface. It is found that PEO:LiDFOB had the highest capacity retention, which is correlated to having the least decomposition at the interface. This indicates that the lower SPE decomposition at the interface still creates a more effective decomposition layer, which is capable of preventing further electrolyte decomposition. Moreover, the PES results indicate formation of polyethylene in the SEI in cells based on PEO electrolytes. This is supported by AIMD that shows a polyethylene formation pathway through free-radical polymerization of ethylene.

Do non-coordinating polymers function as host materials for solid polymer electrolytes? The case of PVdF-HFP.

In the search for novel solid polymer electrolytes (SPEs), primarily targeting battery applications, a range of different polymers is currently being explored. In this context, the non-coordinating poly(vinylidene fluoride-co-hexafluoropropylene) (PVdF-HFP) polymer is a frequently utilized system. Considering that PVdF-HFP should be a poor solvent for cation salts, it is counterintuitive that this is a functional host material for SPEs. Here, we do an in-depth study of the salt dissolution properties and ionic conductivity of PVdF-HFP-based electrolytes, using two different fabrication methods and also employing a low-molecular-weight solvent analogue. It is seen that PVdF-HFP is remarkably poor as an SPE host, despite its comparatively high dielectric constant, and that the salt dissolution properties instead are controlled by fluorophilic interactions of the anion with the polymer.

A method for modelling polymer electrolyte decomposition during the Li-nucleation process in Li-metal batteries.

Elucidating the complex degradation pathways and formed decomposition products of the electrolytes in Li-metal batteries remains challenging. So far, computational studies have been dominated by studying the reactions at inert Li-metal surfaces. In contrast, this study combines DFT and AIMD calculations to explore the Li-nucleation process for studying interfacial reactions during Li-plating by introducing Li-atoms close to the metal surface. These Li-atoms were added into the PEO polymer electrolytes in three stages to simulate the spontaneous reactions. It is found that the highly reactive Li-atoms added during the simulated nucleation contribute to PEO decomposition, and the resulting SEI components in this calculation include lithium alkoxide, ethylene, and lithium ethylene complexes. Meanwhile, the analysis of atomic charge provides a reliable guideline for XPS spectrum fitting in these complicated multicomponent systems. This work gives new insights into the Li-nucleation process, and experimental XPS data supporting this computational strategy. The AIMD/DFT approach combined with surface XPS spectra can thus help efficiently screen potential polymer materials for solid-state battery polymer electrolytes.

Rapid determination of solid-state diffusion coefficients in Li-based batteries via intermittent current interruption method.

The galvanostatic intermittent titration technique (GITT) is considered the go-to method for determining the Li+ diffusion coefficients in insertion electrode materials. However, GITT-based methods are either time-consuming, prone to analysis pitfalls or require sophisticated interpretation models. Here, we propose the intermittent current interruption (ICI) method as a reliable, accurate and faster alternative to GITT-based methods. Using Fick's laws, we prove that the ICI method renders the same information as the GITT within a certain duration of time since the current interruption. Via experimental measurements, we also demonstrate that the results from ICI and GITT methods match where the assumption of semi-infinite diffusion applies. Moreover, the benefit of the non-disruptive ICI method to operando materials characterization is exhibited by correlating the continuously monitored diffusion coefficient of Li+ in a LiNi0.8Mn0.1Co0.1O2-based electrode to its structural changes captured by operando X-ray diffraction measurements.