Overview of Lithium Ion Battery Electrode Development
4.1 Candidate Cathode Material with High Voltage for Lithium Ion Batteries
4.1.1 Introduction
- Early Battery Development (Mid-1970s):
- Whittingham proposed TiS2 for non-aqueous secondary batteries due to its high specific and power density.
- TiS2 has good metallic properties and allows reversible lithium intercalation.
- The discharge voltage was around with current densities of .
- Commercialization was limited to button-size cells due to issues with practical non-aqueous batteries and reversible lithium deposition.
- LiCoO2 Introduction (1980):
- Mizushima et al. suggested LiCoO2 with a-NaFeO2 structure as a cathode material.
- LiCoO2 is structurally more stable than TiS2.
- The open-circuit voltage of LiCoO2/Li cell was approximately , about twice that of LiTiS2.
- The theoretical energy density was .
- It took almost 10 years to commercialize LiCoO2 in batteries.
- Current Demands and Challenges:
- Advancements in microelectronics and electric vehicles (EVs, HEVs) require low-cost, environmentally friendly, and thermally stable lithium-ion batteries.
- These batteries need high power and energy density.
- LiMn2O4 as a Potential Cathode Material:
- Lithium manganese spinel (LiMn2O4) is considered attractive due to these demands.
- LiMn2O4 can be cycled at a range but suffers from gradual capacity degradation, especially above .
- Capacity fading is attributed to Mn dissolution into the electrolyte via the disproportionation reaction: .
- Strategies to Improve LiMn2O4 Stability:
- Bittihn et al. explored Co substitution for Mn in LiMn2O4, which improved cycleability and formed LiCoMn2O4.
- Research group systematically studied LiMyMn2-yO4 (M=Co,Cr,Ni) in the range, finding better cycleability than LiMn2O4, even with small substitutions.
- Substitution enhances spinel stability by strengthening M-O bonds over Mn-O bonds in the MO6 octahedra.
- Thermodynamic stability of MO6, including MnO, in spinels was presented using a Born-Haber cycle with a simple ionic model.
- Molecular dynamics (MD) calculations, considering a partially ionic model, further supported the structural stability of spinels.
- Development of 5V Class Cathodes:
- Extensive charge-discharge studies have been conducted using partially substituted spinel oxides in the class.
- For large-scale EV batteries, operating voltage is typically , making the number of single cells critical.
- Higher single-cell discharge voltage reduces the required cell stack number and improves quality control.
- Overcoming Electrolyte Limitations:
- The discovery of highly oxidation-resistant electrolytes has enabled the study of intercalation compounds above .
- Electrolytes like 1M LiPF6 or LiClO4 in EC and DMC or EC and DEC mixtures are used.
- High Voltage Spinel Materials:
- Several research groups have reported transition-metal substituted spinel materials, , (M=Cr³+, Fe³+, Ni²+, Co³+) with plateaus around .
- These spinels typically exhibit two plateaus at around and .
- The plateau corresponds to reduction to , while the plateau is due to to or to reduction.
- In , plateaus are present at around and , corresponding to the and redox couples, respectively.
- LiNi0.5Mn1.5O4 Composition:
- The ideal composition is of special interest because Mn remains at +4 oxidation state.
- This results in only the redox couple being active during charge/discharge.
- However, synthesis difficulties exist, as high calcination temperatures above can partially reduce Mn from to .
- Impurities like NiO and can also form during solid-state reactions.
- Myung et al. synthesized pure using an emulsion drying method.
- Kim et al. synthesized using a molten salt method at and obtained a single cubic phase of after annealing at in air.
- Structural Order-Disorder Effects: These two compositions exhibit distinct differences in charge-discharge properties. These differences appear to come from structural order-disorder, which affect Li diffusion.
4.1.2 Crystal Structure of Spinel Type Phase
- Spinel Structure of LiNi0.5Mn1.5O4:
- Blasse first reported the structure of stoichiometric .
- Gryffroy et al. proposed it was ordered-type (P4332) using ND and IR measurements, supporting Blasse's proposal.
- Similar results were later obtained from Li MAS-NMR measurements.
- Kim et al. determined the structure of and nonstoichiometric by careful XRD and electron diffraction measurements.
- Their results on were consistent with others, while belongs to cubic spinel(Fd3m).
- Synthesis of LiNi0.5Mn1.5O4 in Research Group:
- Stoichiometric was synthesized by mixing , NiO, and with a molar ratio of Li : Ni : Mn = 1.0 : 1.5 : 1.5 in ethanol.
- NiO was obtained by decomposing in air at for 24h.
- was obtained by heating at for 48h to prevent composition deviation due to hygroscopic properties of .
- After evaporating ethanol, condensed materials were heated at in atmosphere for 4 days and slowly cooled to room temperature at a rate of .
- Nonstoichiometric LiNi0.5Mn1.5O4-δ Synthesis:
- Nonstoichiometric was synthesized by heating the stoichiometric compound at under controlled partial pressure of oxygen.
- The pressure was controlled by changing the mixing ratio of and using and - balanced gas.
- δ values were evaluated from the weight change using TG (Seiko Denshi, EXSTAR6000TG/DTA).
- To obtain data for electrochemical and structural properties at various δ values, nonstoichiometric samples were quenched from to by controlling partial pressures of oxygen in the same manner.
- Slight amount of impurity coexists.
- Disordered-type cubic spinel (Fd3m) is observed in the range of δ values between 0 and 0.31.
- Lattice parameter changes with lithium intercalation.
- Ordered has three cubic phases, while disordered forms only two phases.
- Neutron Diffraction (ND) Patterns:
- ND patterns were obtained for (δ=0.03 and 0.21).
- ND patterns of (P4332) and (Fd3m), estimated by diffraction simulation, were included as references.
- Neutron diffraction was carried out at KEK, Neutron Scattering in Tsukuba, Japan, using the TOF method with the Vega ND apparatus.
- Simulation of the diffraction data was performed.
- As δ values increase, the peaks belonging to ordered (α≈0) gradually diminished, supporting electron diffraction data.
- This is attributed to increased disorder due to an increase in defects.
- Defects in Spinel Oxides:
- Defects in spinel oxides has been discussed by many researchers.
- Density measurements of proposed that the defective structures synthesized under specific oxygen partial pressure are metal excess structures rather than oxygen defects.
- Metals in the octahedral 16d sites may migrate to 16c octahedral sites, which are usually vacant in stoichiometric spinels, and the same would happen in .
- Koyama et al. evaluated the formation energy for , including oxygen-vacancy and metal excess types, using first principles plane-wave pseudopotential calculation, which showed that metal excess defects had the lowest formation energy and were consistent with density measurement results.
4.1.3 Electrochemical Properties of LiNi0.5Mn1.5O4 and LiNi0.5Mn1.5O4-δ
- Charge-Discharge Curves:
- First charge-discharge curves for (δ=0-0.31) at current density were obtained.
- Even at δ=0, a small shoulder at around is observed, which increases with δ values.
- This shoulder is due to the redox reaction, indicating small amounts of exist even in stoichiometric .
- The plateau at around corresponds to the redox reaction.
- Similar charge-discharge profiles have been reported in literature.
- The discharge voltage is 15% higher than range cathodes.
- If a proper stable electrolyte enduring oxidation around is developed, could be a good high-voltage cathode candidate.
- Rate Performance:
- Rate performance of (δ=0 and 0.31) shows nonstoichiometric samples have excellent high-rate property.
- Lithium diffusion in disordered spinels is easier than in ordered ones.
- This may depend on the lithium diffusion path.
- Coulombic Potential Calculation: Site-potential energy for lithium migration in orderd and disorded spinels is caluculated using a simple Coulombic potential calculation
4.1.4 Coulombic Potential Calculation of Diffusion Path for Ordered (P4332) and Disordered (Fd3m) Spinels
- Lithium Diffusion Paths:
- In the disordered spinel (Fd3m), lithium at the tetrahedral 8a site moves to a vacant octahedral 16c site, resulting in an 8a-16c diffusion path.
- In the ordered spinel (P4332), the octahedral vacant 16c sites are split into ordered 4a and 12d sites (1:3 ratio), forming 8c-4a and 8c-12d diffusion paths.
- Coulombic Potential Calculation:
- Electrostatic potential calculation uses the Coulombic potential equation:
, where and are ion charges, and is the distance between ions i and j. - Z values are set as: (-2), (+4), (+2), (+1).
- Coulomb potentials for lithium diffusion paths in ordered and disordered spinels were calculated and normalized at lithium sites.
- Electrostatic potential calculation uses the Coulombic potential equation:
- Potential Curves and Diffusion:
- Potential curves for each diffusion path are similar, with a maximum at the vacant site (saddle point).
- The potential at each bottleneck is ordered as: 8c-4a (order) < 8a-16c (disorder) < 8c-12d (order).
- In the ordered spinel, both the easiest (8c-4a) and most difficult (8c-12d) diffusion paths coexist.
- Disordered-type containing 8a-16c enables easier diffusion.
- Data suggests that lithium diffusivity may increase with increasing δ values.
- Relationship to Electrochemical Properties:
- These calculations support the experimental charge-discharge property, where disordered spinels maintain higher capacity than ordered spinels.
- Introducing defects to improves the rate property of charge-discharge.
4.1.5 Conclusions
- and were synthesized via solid-state reactions at .
- ND data confirmed that the stoichiometric spinel oxide had an ordered structure (P4332), while the nonstoichiometric oxides formed a disordered structure (Fd3m).
- The disordered spinel showed better rate property than the ordered spinel.
- Calculated Coulombic potential for diffusion paths suggested that diffusivity is higher in the disordered spinel.
- The calculation results aligned well with the experimental electrochemical properties.
4.2 Electrode/Electrolyte Interfaces in All-Solid-State Lithium Ion Batteries
4.2.1 Introduction
- Background: Solid electrolytes have gained attention for all-solid-state lithium-ion batteries due to their potential advantages.
- Advantages of Solid Electrolytes: Safety (non-flammable compared to liquid electrolytes) and high energy density for small portable and large-scale batteries.
- Focus of the Chapter: Two types of solid electrolytes: inorganic glass electrolytes and polymer electrolytes.
- Inorganic Glass Electrolytes:
- SiS2-Li2S-Li4SiO4 glass electrolyte exhibits high lithium-ion conductivity (up to ), excellent ionic conductivity (lithium ion transport number of 1.0), and a wide electrochemical window.
- Polymer Electrolytes:
- Polyethers complexed with lithium salts are used, but they suffer from low ionic conductivity and transport number due to:
- Low dissociation degree of lithium salts: low concentration of carrier ions.
- Strong interactions between lithium ions and polymer chains (e.g., ether oxygen atoms), leading to low ionic mobility.
- However, polymer electrolytes are flexible and easily make good contact between the electrode and electrolyte.
- Polyethers complexed with lithium salts are used, but they suffer from low ionic conductivity and transport number due to:
- Research Focus: While many studies target improving bulk properties of electrolytes, few focus on charge transfer reactions at the electrode/solid electrolyte interface.
- Importance of Interfacial Reaction Rate: Enhancing the interfacial reaction rate is crucial for fabricating high-power-density batteries, as battery performance is related to charge-transfer rates at the interface and the electrolyte's bulk conductivity.
4.2.2 Control and Charge Transfer Process at Electrode/Electrolyte Interface
- Three Main Topics:
- (i) Factors determining the charge transfer reaction rate at polyether-based electrolyte/electrode interfaces.
- (ii) Enhancement of ionic conductivity and lithium-ion transport number in polymer electrolytes by adding PEG-borate ester (a Lewis acid).
- (iii) Enhancement of the reaction rate by adding PEG-borate ester.
(i) Study of Factors Influencing the Charge Transfer Reaction Rate at the Polyether-Based Electrolyte/Electrode Interface
- Goal: Clarify the reaction rate at interfaces and the ionic conduction mechanism.
- Method: Investigated the electrokinetics of the couple reaction on lithium metal electrodes in poly(ethylene glycol) dimethyl ether (PEGDME)-based electrolytes.
- PEGDME solutions with are used as a model for polymer electrolytes, resembling the amorphous conducting phase in high-molecular-weight PEO.
- Experiment: Measured charge transfer reaction rate using PEGDME (average molecular weight 500, PEGDME500) with varying amounts of higher molecular weight polyether, PEGDME1000.
- Conclusion: Viscosity is a major factor influencing the charge transfer reaction rate, explained by Marcus microscopic theories.
- Technique:
- Double potential step chronoamperometry was used to study lithium deposition and dissolution on the microelectrode.
- Coulombic efficiency of lithium deposition and dissolution was over 90%, due to fresh deposition and immediate dissolution on the microelectrode.
- Analysis:
- Steady-state currents (i) at various overpotentials (E) were analyzed using the Allen-Hickling equation:
- Where F is the Faraday constant, R is the gas constant, T is the absolute temperature, α is the transfer coefficient, and is the exchange current density.
- A linear relationship was observed in Allen-Hickling plots in the range of overpotentials around -150 to 170 mV vs. .
- Exchange current density and transfer coefficient were estimated from the intercept and slope of the Allen-Hickling plots.
- The exchange current densities decreased with increasing amounts of PEGDME1000.
- Steady-state currents (i) at various overpotentials (E) were analyzed using the Allen-Hickling equation:
- Exchange Current Density Equation:
- The exchange current density can be expressed as:
- Where is the activity of lithium ion, is the activity of lithium metal, and is the standard rate constant.
- According to advanced theory:
- A is the pre-exponential factor,
- is the Gibbs activation energy of the reaction,
- is the longitudinal relaxation time of the solvent.
- Since all values of the transfer coefficient are close to 0.5 and is unity, the equation for the exchange current density is simplified to which indicates that the exchange current density depends on the activity of lithium ions, the Gibbs activation energy, and the longitudinal relaxation time of the solvent.
- Raman Spectroscopy:
- Raman spectra were used to investigate the activity of lithium ions.
- The symmetrical stretching mode of () at around was used.
- The band consists of three components:
- Free ions (~).
- Ion pairs (~).
- Multiply aggregated ions (~).
- The band consists of three components:
- The fraction of free ions (dissociated ions) corresponds to the ratio of free ion area to the areas of all components.
- The fractions of free ions were almost the same (around 20%) even with increasing amounts of PEGDME1000, indicating constant lithium ion activity.
- The solvation state of lithium ions in PEGDME should be similar, since ether oxygens act as donor polar groups.
- The ratio of free ions was almost constant.
- Gibbs Activation Energies:
- Gibbs activation energies were calculated from the slope of linear relationship plots. The electrokinetics of the couple reaction involves a desolvation/solvation process where the Gibbs activation energy may depend on desolvation/solvation energy.
- The Gibbs activation energy depends on the solvation state of lithium ions.
- Since lithium ions are incorporated into the helix structures of EO chains, activation energy is independent of the amount of PEGDME1000.
- Where is the longitudinal relaxation time of the solvent.
- Solvent Effects:
- The obtained results clarified that the activity of lithium ions and the Gibbs activation energy do not depend on the ratio of PEGDME500/PEGDME1000 in the electrolyte solutions, while the exchange current densities decreased with increasing amounts of PEGDME1000.
- The Solvent effect on charge transfer kinetics have elucidated that the reaction rate involves the frequency of solvent reorientation in the desolvation/solvation process of solvent dynamics where this frequency is inversely proportional to the longitudinal relaxation time.
- : Debye relaxation time.
- :static dielectric constant.
- : infinite frequency dielectric constant.
- Therefore, , where viscosity is related to the frequency of solvent reorientation in the desolvation/solvation process.
- Viscosity Relationship:
- Viscosity of PEGDME500-based solutions increased with increasing amounts of PEGDME1000.
- The relationships between the exchange current densities and the viscosities of the solutions showed that the exchange current density was inversely proportional to the viscosity.
- Therefore, the viscosity of the solvent is an important factor in determining the exchange current density of the polyether-based electrolytes, where high rates of charge transfer reactions at the interfaces is difficult when using high-viscosity electrolytes and especially polymer electrolytes.
- Viscosity of PEGDME500-based solutions increased with increasing amounts of PEGDME1000.
- Conclusion:
- Exchange current densities in PEGDME500-based electrolytes decreased with increasing amounts of PEGDME1000.
- Raman spectroscopic studies of the electrolytes revealed that the activity of lithium ions was almost the same even with the addition of PEGDME1000.
- Gibbs activation energies for charge-transfer reaction did not change with the amount of PEGDME1000, suggesting that solvation states of lithium ions were similar in the electrolytes to those in PEGDME whose molecular weight range was 500 to 1000.
- The viscosity of the electrolytes increased with increasing amounts of PEGDME1000.
- Inversely proportional relationships were found between the exchange current densities and the viscosities of the electrolytes.
- Therefore, the viscosity of the electrolyte solution is key for the charge transfer reaction rate at the polyether-based electrolyte/electrode interfaces.
- It will be difficult to get high rates of charge transfer at the interfaces when using polymer electrolytes, or other highly viscous electrolytes.
(ii) Enhancement of Ionic Conductivity and Transport Number of Lithium Ions in Polymer Electrolytes by the Addition of PEG-Borate Ester as a Lewis Acid into the Electrolytes
- Objective: To enhance ionic conductivity and transport number of lithium ions in solid polymer electrolytes for practical electrochemical devices.
- Strategy: Enhancement of the dissociation of supporting salts and decreasing the mobility of the counter anions.
- Lewis Acid Addition: By adding Lewis acid compounds one can enhance the transport number of lithium via interacting with Lewis base anions as the Lewis acid compounds trap the counter anions.
- Boron Compounds: Reportedly act as anion receptors due to Lewis acid character.
- Investigation Focus: The effect of the Lewis acid character of borate ester groups, which are fixed to the chains of the matrix polymer, on ionic conductivity and transport number of lithium ions in the polymer electrolytes. The relationships between the concentration of borate ester groups and ionic conductivity or transport number of lithium ions in the polymer electrolytes were examined.
- Polymer Preparation:
- Backbone polymer prepared from polyethylene glycol and boric acid anhydride.
- Polyethylene glycols were with various molecular weights (PEG150, PEG200, PEG400, PEG600, whose average molecular weight was 150, 200, 400, 600, respectively).
- LIN(CF3SO2)2, LiCF3SO3 and LiClO4 were used for supporting salts without further purification.
- Polymer electrolyte composition of polymer electrolyte, for example, prepared from PEG200, boric acid anhydride and LiN(CF 3 SO2)2, is represented as PEG200-B2O3 + LIN(CF3SO2)2.
- Computational Methods
Ab initio Hartree-Fock (HF) self-consistent field molecular orbital calculations and density functional theory (DFT) calculations were performed with Gaussian98.
Calculations for geometry optimizations were carried out at the HF level of the theory using the standard 3-21G basis set. Subsequently, single point calculation to investigate the energies of the optimized geometries was performed by using DFT methods with the B3LYP form for the exchange-correlation function and the 6-311G** basis set. - Glass Transition Temperature:
- of the polymer electrolytes decreased with increasing molecular weight of PEG, resulting in decreasing the concentration of the borate ester groups.
- The mobility of the polymer chains increases due to decreesing concentration of crosslinking points.
Ionic motion increased with descreased crosslinking
- Ionic Conductivity:
- Increase in ionic conductivity of the polymer electrolytes is observed with increasing molecular weight of PEG.
- The highest ionic conductivity was found for the polymer electrolyte prepared from PEG600, to be at and at .
- Free volume Theory: The free volume theory of polymers states which leads to the temperature dependence for segmental motion of polymer chains expressed by the Williams-Landel-Ferry (WLF) relationship:
Where σ(T) and σ(Tg) are the conductivity at temperatures T and Tg, respectively, and C1, C2 are the WLF parameters for the temperature dependence of ionic conductivity. However, the ionic conductivity at TÅ, σ(Tg), is difficult to measure in the present experiments because σ(Tg) is too low to measure with a complex impedance measurement, and therefore, 10 °C was selected as the reference temperature, To. The equation is rewritten as follows:
The parameters in equation (1) are calculated as follows:
WLF Plots:
- The WLF plots for the ionic conductivity of the polymer electrolytes show that [log(σ(T)/σ (To))] is varies linearly with, which indicates that the temperature dependence of the ionic conductivity for the polymer electrolytes follows a WLF-type equation as the temperature dependence of the ionic conductivity for the polymer electrolytes was dominated by that of the segmental motion of polymer chains. Therefore from the WLF plots and Tå values, the increase in ionic conductivity with increasing molecular weight of PEG is due to the increase in ionic mobility, where parameters were found to be close to the universal values of WLF parameters, = 17.4 and = 51.6/K47
I is the current, AV is the applied potential, R₁ is the bulk resistance,
Re is the interface resistance, and 0 and ∞ refer to the initial and steady states, respectively.It was found that the transport number of the polymer electrolytes containing LiCF3SO3 or LiClO4 was higher than that of the electrolyte with LIN(CF3SO2)2, showwing that the Lewis acidity of the borate ester groups in polymer electrolytes affects the anions, ionic conductivity and Tg of the polymer electrolytes containing LiCF3SO3 or LiClO4 were examined.
The HOMO exists around the oxygen atoms, where electron density around the oxygen atoms apart from the boron atoms is higher than that around the oxygen atoms next to the boron atom, where lithium ions interact more strongly with the oxygen atoms that are apart from the boron atom than to other oxygen atoms.
On the other hand, the LUMO exists perpendicular to the BO3 plane, which indicates a location for interaction between the boron atom and a Lewis basic anion.The order of both X and 7 of the anions is CF3SO3 > C104 > N(CF3SO2)2. Therefore, CF3S03 or C10 should interact more strongly with the PEG-borate ester than N(CF3SO2)2 does as The boron atoms of the backbone polymer act as a