Introduction
Lithium manganese oxide (LiMn
2O
4) is the primary cathode material being exploited for high power and high energy density applications like electric vehicles [
1] and renewable energy storage [
2]. LiMn
2O
4 cathodes possess several advantages, such as low cost, non-toxicity, and simple preparation process [
3] as compared to alternative cathode materials. However, there are well-documented problems with capacity fade of lithium ion batteries containing LiMn
2O
4 due to active material loss, caused by dissolution of divalent manganese from the cathode surface and deposition of dissolved manganese onto the anode [
4]. The LiMn
2O
4 spinel consists of manganese ions in two oxidation states (Mn
3+ and Mn
4+) in equal proportions to give a net manganese oxidation of+ 3.5 [
5]. Conventional wisdom holds that trivalent manganese ions are unstable to acid attack and are lost from the cathode surface into the electrolyte following charge disproportionation through the Hunter reaction [
4]
The dissolution of divalent manganese from the surface of the cathode into the electrolyte leads to a decrease in the effective amount of cathode material available for lithium intercalation. Furthermore, deposition of dissolved manganese onto the anode decreases the Li content of the anode by consuming electrons.
A key to improving the capacity and service life of LiMn
2O
4 based cells is to understand the reactions that occur at the surface of the cathode leading to dissolution of manganese. The cathode-electrolyte interface has only recently been considered as a potential impediment in the long-term viability of Li-ion cells [
6]. Advances in imaging techniques such tunneling electron microscopy have revealed the presence of an extremely thin (2 to 5 nm thickness [
7]) layer of compounds on the cathode surface. Electrochemical impedance spectroscopy measurements have confirmed that this film causes increased cell impedance [
8], implying that the compounds formed at the cathode-electrolyte interface are poor ion conductors, thus limiting the diffusivity of Li
+ through the layer [
9]. The cathode-electrolyte interface is speculated to consist of a layer of electrolyte oxidation products and cathode active material salts formed by reduction of manganese [
10].
Ex situ analysis of the cathode surface by infrared spectroscopy suggests the presence of polyethylene carbonate (PEC), perhaps deposited by oxidative polymerization of ethylene carbonate (EC) in the electrolyte. However, this finding could not be confirmed by X-ray photoelectron spectroscopy due to overlap between the PEC and polyvinylidine fluoride binder peaks [
11], which yielded an inconclusive spectroscopic analysis of the cathode surface. Oxidation of other common electrolyte solvents, such dimethyl carbonate (DMC) and fluoroethylene carbonate (FEC) have been speculated but remains unconfirmed. Chemical characterization of the cathode-electrolyte interface is challenging because the concentrations of the interface compounds are below the detectable limits of nuclear magnetic resonance spectrometry and gas chromatography/mass spectrometry instruments [
11]. Although manganese deposition on the anode surface [
12] provides compelling evidence of manganese dissolution from the cathode, the speciation and reactions responsible for this at the cathode solid-electrolyte interface (SEI) layer are not yet well understood.
The cathode-electrolyte interface was found to have a layered structure similar to the anode-electrolyte interface through argon ion sputtering experiments [
13], but its composition and electrochemical characteristics differ markedly from the anode SEI layer. Experiments have shown that the amount of Mn
2+ present in the electrolyte increases with cycle number [
14], suggesting that the cathode SEI layer formed during the initial cycling does not prevent further decomposition of electrolyte, in contrast to the anode SEI layer. The dynamic behavior of the cathode SEI layer is attributed to organic radicals and ions formed via electrolyte oxidation, which further combine with solvent molecules in a chain reaction [
15]. Hence, the cathode SEI layer can limit battery performance by consuming electrolyte continuously as cycling proceeds. The LiMn
2O
4 surface is also highly susceptible to acid attack by HF produced via hydrolysis of LiPF
6 salt by trace water present in the electrolyte [
16]. While HF is a putative trigger for Mn dissolution, the reaction mechanisms by which this occurs remain unknown. To improve battery performance and cycle life, a better understanding is therefore sought of the reactions that occur at the cathode surface that cause the loss of active material.
Various strategies have been adopted to improve cathode stability under cycling and decrease active material dissolution; these strategies purportedly involve the formation of a passivating cathode SEI layer [
17]. The most popular approach is cathode modification via active material doping or cathode surface coating, where in Mn sites in the LiMn
2O
4 spinel are doped with other transition metals such as Ni [
18], Fe [
18], Co [
19], Mg [
20], Zn[
20], or Al [
21] to reduce the number of Mn
3+ cations that are susceptible to acid attack. Spinel doping with Ni has proven especially successful for increasing the overall electrochemical performance of the cathode and reducing the Mn dissolution rate. However, nickel-doped spinel exhibits capacity loss after prolonged cycling due to surface film formation [
22]. Another somewhat effective modification is the application of a surface coating to protect the cathode active material, either by altering the chemistry at the cathode-electrolyte interface, or by acting as a physical barrier between the active material and electrolyte. Al
2O
3 is the most commonly used coating for transition metal based cathodes as it is known to scavenge HF, reducing the acidity of the cathode surface, and it also reduces electrolyte decomposition at the cathode surface [
23]. Surface coatings, however, have low lithium ion conductivity, necessitating the optimization of coating thickness in order to act as an effective barrier to electrolyte oxidation while still conducting Li ions [
24,
25]. The use of electrolyte additives is another strategy to suppress electrolyte decomposition on the anode surface. Although a few compounds such as vinylene carbonate have been reported as potentially effective additives to reduce electrolyte degradation on LiFePO
4 cathodes [
26], these results could not be replicated on other cathode chemistries [
27].
All of the aforementioned strategies have been designed based on indirect evidence of increased cell performance, without knowing the changes in actual surface chemistry induced by cathode or electrolyte modifications. To devise more efficient cathode improvement strategies, greater insight into the composition and dynamics of the cathode electrolyte interface is essential. Considering the inherent difficulty of characterizing the composition of the cathode electrolyte interface layer using experimental techniques, computational methods are a resource to probe otherwise inaccessible nanoscale phenomena. Ogawa et al. [
28] used quantum chemical molecular dynamics to study the transport of Li ions at the cathode electrolyte interface and determined that Li ions are solvated by EC molecules at the cathode surface. This study however did not consider differences in the chemical composition of the cathode-electrolyte interface and bulk electrolyte. Density functional theory methods have been used to study the oxidation of isolated solvent molecules and salt anions [
29,
30] that do not explicitly interact with the cathode. Leung [
31] conducted the first study on EC decomposition on LiMn
2O
4 surfaces, revealing the generation of surface hydroxyl species (and possibly H
2O molecules) via H abstraction from EC as the primary reaction mechanism on the surface of the cathode. While this study provided invaluable insights into the chemical nature of the cathode surface, it is important to investigate the reactions occurring at the cathode surface in the presence of the complete electrolyte. Among the important research questions yet to be answered are the products of the oxidative decomposition of the solvent molecules; the fate of acidic species (HF) on the cathode surface; and the complete reaction mechanism generating the cathode electrolyte interface layer. In this work, reactive force field based molecular dynamics is therefore applied to investigate the reactions occurring at the LiMn
2O
4 cathode surface and the mechanisms of manganese dissolution.
Methods
The ReaxFF model is employed to conduct molecular dynamics simulation of the surface-electrolyte interface between a lithium manganese oxide cathode and a nonaqueous electrolyte with a force field that allows for changes in bond order and bond length during the course of the simulation. The ReaxFF methodology is described in full in previous work [
32–
34] and was recently extended to model pairwise interactions of manganese [
35] and fluorine [
36] with carbon, hydrogen, oxygen, and lithium. The developed reactive force field is used to investigate the cathode-electrolyte interface layer in electrochemical cells with LiMn
2O
4 cathode, ethylene carbonate/dimethyl carbonate electrolyte, LiPF
6 salt, and 5% fluoroethylene carbonate additive. To probe manganese dissolution due to acid attack, HF is introduced into the simulation cell. Conventional wisdom holds that HF is generated within the battery when LiPF
6 molecules decompose and react with trace water molecules present in the electrolyte according to the following reactions [
16]
The focus of this work is to investigate the fate of the HF molecules on the cathode surface and not simulate the formation reactions of HF. Hence, in the interest of reducing the complexity of the study, the simulations assume the product of LiPF6 decomposition (HF) to be present in the battery without simulating the actual decomposition mechanisms. This assumption overlooks the effects of other LiPF6 decomposition products on the cathode electrolyte interface. The inactive role of these decomposition products has not been proved or speculated but this assumption simplifies the system without compromising the purpose of this study.
The formation of the cathode electrolyte interface layer was simulated using a cathode half-cell composed of a LiMn
2O
4 crystal with its 100 surfaces exposed to the electrolyte consisting of 208 EC molecules, 480 DMC molecules, 20 FEC molecules and 60 LiPF
6 molecules. The LiMn
2O
4 crystal measuring 3.5 nm × 3.5 nm × 3.5 nm is located at the center of the simulation box with dimensions of 10.5 nm × 3.5 nm × 3.5 nm, the rest of the simulation box is filled with the electrolyte as shown in Fig. 1. This configuration helps with maintaining periodic boundary conditions in all three directions without any discontinuities. Bulk electrolyte character is maintained at distances greater than 2 nm from the electrode surfaces even after 2 ns, the maximum time period in this study. The equilibrium electrolyte density in the simulation cell is 1.26 g/mL with 3:7 mol ratio [
37] of ethylene carbonate (EC) and dimethyl carbonate (DMC), 1 mol/L lithium hexafluorophosphate and 5% fluoroethylene carbonate.
The simulation cells were maintained at 330 K under the NVT ensemble using the Nosé Hoover thermostat with a temperature rescaling constant of 10 fs. The starting configuration of the electrolyte molecules was generated using the PACKMOL [
38] code to make sure short range repulsive forces do not disrupt the simulation during the first time step. Each MD simulation was initiated from an energy-minimized structure and was equilibrated to the simulation temperature for 100 ps prior to the production run for 2 ns with a 0.2 fs time step.
The fate of HF on the cathode surface was studied by running a duplicate simulation as described above and introducing HF molecules into the simulation after 1.5 ns of cathode electrolyte interface layer formation. The HF molecules were introduced at random locations within a distance of 10 Å from the surface of the anode. HF molecules are introduced into the cathode electrolyte interface layer at a rate of 0.5 molecules/ps with a total of 10 molecules introduced into the simulation cell. After the introduction of HF molecules, the simulations were run for an additional 500 ps.
The ReaxFF based MD simulations were implemented using the LAMMPS parallel molecular dynamics code [
39]. Molecular analysis of the simulations, which is to identify SEI layer compounds, was performed using a bond order cutoff of 0.2 as implemented in a home built MATLAB code [
32].
Results
This work sets out to investigate the reaction mechanisms that result in the formation of the cathode-electrolyte interface layer. The reaction pathways of electrolyte molecule oxidation will be elucidated, the role of HF in the dissolution reaction will be determined and finally, the overall reaction mechanism on the surface of the cathode will be discussed.
Electrolyte solvent decomposition mechanisms
The ReaxFF based MD simulations of the cathode-half cell reveal a layer, 10 Å in thickness, consisting mainly of electrolyte solvent oxidation products formed on the cathode surface after 2 ns of simulation at 330 K. Figure 2 shows the primary compounds identified on the cathode surface. The ReaxFF model predicts a highly organic cathode-electrolyte interface consisting of compounds very similar to those speculated by other researchers [
11,
40–
42]. At least three different oxidation products of EC were found in the simulation cell after 2 ns depending on the reaction pathway followed. The oxidation of EC also triggers polymerization, yielding large organic molecules as observed by other studies on EC oxidation [
40]. While EC is most susceptible to oxidation on the cathodes surface, oxidation products of DMC and FEC are also present in the layer.
Figure 3 identifies the primary oxidation pathways of the three solvent molecules. The oxidation pathway of the solvents has been found to be very similar to the pathway reported by Leung [
31], proceeding via alkyl hydrogen abstraction and consequent lattice oxygen removal. The surface lattice oxygen atoms take up the hydrogen atoms from the alkyl groups of the solvent molecules resulting in a surface hydroxyl group and an organic radical adsorbed onto the cathode surface. After about 150–600 ps (depending on the reacting solvent molecule; see Fig. 2), the organic radical pulls out a lattice oxygen atom which is made available due to the reduced oxidation state of the surface manganese atoms in the presence of a nearby hydroxyl species. The removal of lattice oxygen results in the formation of organic aldehyde groups, ester groups and polycarbonate species in the case of EC oxidation. DMC and FEC are found to be oxidized to only one type of product with no evidence of polymerization. The organic radicals formed after H abstraction from EC often react with available surface hydroxyl groups resulting in the formation of alcohol groups. Unlike the anode electrolyte interface layer, no inorganic carbonates or organic lithium salts have been predicted by the ReaxFF model in the absence of salt decomposition; this is in good agreement with experimental evidence [
43]. The oxidation reaction time for the solvent molecules ranges from 70 ps to 650 ps and increase in the order of FEC<EC<DMC.
The oxidation reaction pathways for all the electrolyte solvent molecules reveal the formation of surface hydroxyl species. The removal of oxygen from the cathode lattice exposes manganese atoms to the bulk electrolyte, which react with the surface hydroxyl species to form manganese hydroxide that is no longer a part of the lattice structure. The hydroxyl molecules are further, found to form water molecules that are likely to hydrolyze the LiPF
6 salt yielding HF. The formation of hydroxyl species and removal of lattice oxygen atoms result in manganese dissolution from the cathode into the electrolyte. The role of lattice oxygen vacancy formation on manganese dissolution has been speculated by previous experimental work using XAS analyses on cathode surfaces [
44]. ReaxFF predictions of lattice oxygen removal during oxidation of electrolyte solvent molecules confirms the role played by oxygen vacancies on manganese dissolution and provides further insights into the reaction mechanisms that lead to lattice oxygen removal.
Fate of HF on the cathode surface
HF in the lithium ion cells is believed to be a byproduct of LiPF6 hydrolysis due to the presence of trace water molecules. But, the results from this study show that water molecules can be formed on the cathode surface during electrolyte oxidation. This finding implies that even in the absence of water as an impurity, HF can be formed by the hydrolysis of LiPF6 salt via water molecules generated through electrolyte oxidation on the cathode surface. Therefore, HF is an inevitable part of the battery chemistry and considering that HF is a putative trigger for manganese dissolution, its fate on the cathode surface is an essential part of the cathode-electrolyte interface chemistry. To investigate its fate on the cathode surface, HF is introduced into the simulation cell after 1.5 ns. A total of 10 HF molecules are introduced within 10 Å from the cathode surface and the cathode interface layer composition is analyzed after 500 ps.
The introduction of HF into the simulation cell results in a consequent rise in the number of surface hydroxyl species on the cathode. Figure 4 shows the primary compounds identified in the cathode electrolyte interface after 500 ps of simulation in the presence of HF molecules. The number of Mn(OH)
2 molecules also increases in the presence of HF. As can be seen in Fig. 2, the manganese atoms that form Mn(OH)
2 are extracted from the cathode lattice structure on bond formation with the surface hydroxyl species. The Mn-OH species observed in the simulation remain adsorbed on the cathode surface, but it is highly likely that manganese atoms are solvated by the electrolyte compounds leading to the dissolution of manganese at longer timescales. The fluoride anions from HF mostly end up as fluorinated salts of manganese (MnF
2) and lithium (LiF), which are the only inorganic salts present on the cathode surface, as observed experimentally [
45]. Organic radicals formed after alkyl hydrogen abstraction from electrolyte solvent molecules also react with HF to yield fluorinated organic solvents.
The presence of HF results in an increased amount of surface hydroxyl species and the formation of inorganic salts on the surface of the cathode. While the hydroxyl species react with manganese atoms displacing the metal atoms from the cathode lattice, the fluoride anions react with metal atoms to form inorganic salts which may help passivate the cathode surface in the long run. ReaxFF analysis yields insights into the acid attack mechanism occurring on the cathode surface and reveals that scavenging hydroxyl molecules formed on the surface of the cathode may lead to the formation of a stable cathode-electrolyte interface layer.
Cathode-electrolyte interface layer formation reaction mechanism
Based on the ReaxFF MD analysis of the LiMn
2O
4 cathode electrolyte interface layer, the formation reaction scheme can be illustrated as shown in Fig. 5 [
16]. Electrolyte solvent molecules undergo alkyl hydrogen abstraction and oxidation via the removal of lattice oxygen atoms on the cathode surface. This oxidation process results in the formation of organic aldehydes, esters and alcohols that compose the cathode electrolyte interface layer. The oxidation reaction also results in the production of surface hydroxyl species that react with manganese atoms displacing them from the lattice. The hydroxyl species also form water molecules. The manganese hydroxide molecules and water molecules are likely to be released into the electrolyte. The formation of water molecules is likely to contribute to the continuous degradation of the cathode electrolyte interface layer due to salt hydrolysis that yields HF molecules [
46], further explaining the non-passive nature of the cathode electrolyte interface layer. The presence of HF in the electrolyte adds to the formation of surface hydroxyl species and yields metal fluorides. Manganese fluoride, lithium fluoride and fluorinated solvent molecules can be found in the cathode electrolyte interface layer. The components of the cathode electrolyte interface layer identified by this work (excluding the compounds produced due to salt decomposition) are similar to those speculated by previous studies [
13,
47]. This work adds to previously established cathode-electrolyte interface reaction schemes [
48,
49] by providing detailed to reaction mechanism. The ReaxFF analysis also establishes that the formation of hydroxyl species is central to the formation and non-passivating nature of the cathode-electrolyte interface layer.
Conclusions
The formation of the cathode-electrolyte interface layer was investigated using ReaxFF based molecular dynamics. The MD simulations reveal that the cathode-electrolyte interface layer is mainly composed of oxidation products of electrolyte solvent molecules such asaldehydes, esters, alcohols, polycarbonates and organic radicals. The electrolyte solvent molecules undergo oxidation by removing lattice oxygen atoms and producing surface hydroxyl species. In the presence of HF, increased number of hydroxyl species, metal fluorides and fluorinated organic solvent molecules were formed. Manganese dissolution was aided by the removal of lattice oxygen atoms that expose manganese atoms to the bulk electrolyte, and the formation of surface hydroxyl species that react with the exposed manganese atoms and displace the atoms from the cathode lattice.
The cathode-electrolyte surface compounds predicted by the ReaxFF model are in good agreement with the experimentally observed products. The ReaxFF analysis yielded insights into the acid attack mechanism occurring on the cathode surface, establishing that the generation of hydroxyl species is central to the formation and non-passivating nature of the cathode-electrolyte interface layer as well as the issue of manganese dissolution. The results of this work imply that developing strategies to scavenge hydroxyl molecules formed on the surface of the cathode will help in the formation of a stable cathode-electrolyte interface layer. Furthermore, this work also developed a detailed formation reaction pathway for the cathode-electrolyte interface layer with insights into the exact chemical compounds that can be found on the cathode surface. This information will be especially useful to future studies looking to bridge gaps between computational and experimental techniques.
Higher Education Press and Springer-Verlag Berlin Heidelberg