Thermodynamics of Subzero Lithium Plating in Graphite Anodes

Potential

As thermal energy drops, electrochemical equilibrium at the graphite negative electrode shifts substantially. In a standard lithium-ion cell operating at room temperature, lithium ions insert between graphene sheets at thermodynamic potentials between 200 mV and 65 mV versus the reference lithium metal electrode (Li/Li+). Because this window remains above zero volts, solvated lithium ions are prevented from directly reducing into metallic lithium.

At subzero temperatures, however, the thermodynamic equilibrium potential shifts according to the Nernst equation while kinetic resistance increases exponentially. When the total overpotential driven by the charging current exceeds this equilibrium potential, the graphite anode’s operational potential drops below 0 V versus Li/Li+.

Negative anode potentials create the thermodynamic driving force for metallic lithium deposition. Below zero volts against lithium metal, two competing reaction pathways occur at the electrode surface: ion intercalation into the graphite host lattice and phase transformation into solid metallic lithium. Intercalation requires ions to desolvate, cross the solid-electrolyte interphase, and diffuse into carbon particles.

Metallic plating bypasses solid-state diffusion altogether, depositing raw metal directly onto active graphite surfaces or SEI defects.

Thermodynamic Thresholds for Metallic Lithium Precipitation

Operating a lithium-ion cell below freezing alters the free energy change of intercalation relative to metal deposition. The equilibrium potential of reaction phases in graphite, such as Stage 1 (LiC6) and Stage 2 (LiC12), exhibits a weak temperature coefficient governed by reaction entropy. Operational cell potential, however, is dominated by dynamic overpotentials.

The thermodynamic threshold for plating occurs the moment local surface potential reaches 0.00 V versus Li/Li+, below which thermodynamic overpotential acts as the driving voltage for metallic lithium nucleation and growth.

For a reference graphite electrode operating at a state of charge corresponding to an equilibrium intercalation potential of 80 mV versus Li/Li+ at 25°C, consider a charge current density of 2.0 mA/cm2 at -20°C (253.15 K). Total negative electrode potential is calculated by subtracting combined kinetic overpotentials from the equilibrium potential. If ohmic overpotential across the electrolyte film is 25 mV, charge-transfer overpotential across the solid-electrolyte interphase is 55 mV, and solid-state diffusion overpotential within the graphite particle reaches 15 mV, aggregate overpotential totals 95 mV.

Subtracting 95 mV from the 80 mV equilibrium potential leaves a net electrode potential of -15 mV versus Li/Li+. Because this net value falls below 0 V versus Li/Li+, metallic lithium deposition becomes thermodynamically favored over intercalation.

At -20°C and a charge rate of 0.5C, graphite anode potential drops below 0 V versus Li/Li+ within 120 seconds of current onset.

Components of Anode Polarization at Subzero Temperatures

Total voltage drop at the negative current collector combines bulk electrolyte resistance, charge transfer barriers, and solid-state migration delays. Ohmic polarization stems from the ionic resistance of liquid electrolyte within separator pores and electrode voids. Charge transfer polarization represents the kinetic barrier of ion desolvation and transport through inorganic and organic SEI layers.

Solid-state diffusion polarization arises from concentration gradients inside spherical carbon particles when lithium accumulates at particle perimeters faster than it migrates toward the core.

Maintaining negative electrode potential above the metallic lithium deposition threshold prevents dendritic growth across cold charging operating windows.

Flux

Mass transport kinetic rates collapse exponentially as electrolyte solutions cool toward their glass transition thresholds. At subzero temperatures, the primary thermodynamic bottleneck shifts from equilibrium potentials to transport properties governed by temperature-dependent activation energies. As ion transport through the liquid electrolyte and solid graphite slows, incoming charge current forces an accumulation of excess charge at the electrode interface.

This charge buildup creates severe polarization, pulling local potential down into the metallic plating regime.

Lithium ion flux through the battery structure encounters three distinct physical media: bulk liquid solvent mixtures, solid-state passivation films, and crystalline carbon lattices. Each medium presents an activation energy barrier governing the temperature dependence of ionic conductivity and diffusion rates. Lowering cell temperature from 25°C to -20°C typically increases liquid electrolyte bulk resistivity by a factor of five while increasing charge transfer resistance at the interphase by up to two orders of magnitude.

Arrhenius Kinetics and Solid-State Diffusion Barriers

Lithium ions move through solid graphite particles along specific crystallographic planes at rates determined by lattice thermal agitation. Solid-state diffusion coefficients inside graphite follow Arrhenius behavior, where diffusion speed scales with temperature through an exponential activation energy factor. Typical activation energies for solid-state intercalation range between 0.35 eV and 0.50 eV depending on graphite particle orientation, defect density, and surface coatings.

At -20°C, the solid-state diffusion coefficient drops to less than one percent of its room-temperature value.

Sluggish solid-state transport creates severe concentration gradients inside individual graphite particles during subzero charging. Lithium ions arriving at particle surfaces cannot penetrate deep into the carbon bulk quickly enough to clear surface intercalation sites. Local surface saturation occurs rapidly, pushing local stoichiometry near the surface close to LiC6 while interior regions remain lithium-deficient.

Surface saturation eliminates open intercalation sites, forcing incoming ionic flux to undergo direct metallic reduction on particle boundaries instead of solid-state insertion.

Desolvation Energy and Charge Transfer Resistance

Stripping solvent molecules from liquid-phase coordination shells requires substantial work before active ions can cross the solid-electrolyte interphase. In carbonate electrolytes, lithium ions exist as solvated complexes coordinated by approximately four solvent molecules, such as ethylene carbonate or ethyl methyl carbonate. The activation energy required to break these coordination bonds and shed solvent molecules before entering the inorganic interphase ranges between 0.60 eV and 0.85 eV.

At subzero temperatures, high desolvation energy creates the dominant kinetic obstruction during charging. Charge transfer resistance increases far faster than liquid-phase ohmic resistance as thermal energy declines. When high current density passes through an interface dominated by large charge transfer resistance, local electrical potential drops sharply.

This voltage drop drives negative electrode potential below 0 V vs Li/Li+, transferring current pathways from intercalating flux to metallic deposition pathways.

• Desolvation Energy Barrier high activation energy required to remove carbonate solvent molecules from solvated lithium ions at subzero interfaces.
• Interfacial Charge Transfer Impedance exponential rise in resistance as ions cross inorganic solid-electrolyte interphase layers under subzero thermal conditions.
• Solid State Lattice Diffusion Delay dramatic drop in lithium ion migration velocity inside graphite inter-planar channels at low temperatures.
• Electrolyte Salt Depletion reduction of mobile charge carriers within separator and porous electrode void spaces due to transport limitations.
• Metallic Lithium Nucleation high localized current density exceeding intercalation flux capacity, precipitating direct elemental metallic formation.

Charging a lithium-ion pack without subzero current derating causes rapid thermal runaway risk, permanent capacity loss, and immediate voiding of battery system safety certifications.

Pore

Current density distribution within a composite electrode depends heavily on local ionic conductivity and tortuosity. Porous graphite anodes consist of compressed active material particles, conductive additives, and polymeric binder matrices filled with liquid electrolyte. In subzero environments, liquid electrolyte inside narrow electrode channels suffers severe reductions in ionic conductivity and salt diffusion coefficients.

The combination of high liquid tortuosity and reduced bulk transport properties produces non-uniform reaction profiles throughout the spatial thickness of the negative electrode.

In optimal room-temperature charging conditions, ionic current penetrates deep into the full thickness of the anode coating, distributing reaction flux evenly across all available active graphite surfaces. Under subzero conditions, liquid-phase transport limitations force the reaction to concentrate almost entirely in active regions directly adjacent to the separator face. Active carbon particles near the current collector foil remain underutilized while graphite particles near the separator undergo extreme localized current loading, driving surface potentials below zero volts and initiating localized lithium metal plating.

Electrolyte Viscosity and Salt Depletion Gradients

Subzero temperatures increase liquid solvent shear resistance while suppressing binary diffusion coefficients. As electrolyte temperature approaches -30°C, linear carbonate solvents experience dramatic viscosity increases that restrict the movement of solvated ions. The liquid-phase ionic conductivity of standard ethyl methyl carbonate formulation blends drops from roughly 10 mS/cm at room temperature down to less than 0.8 mS/cm at -20°C. Concurrently, the salt diffusion coefficient decreases, preventing salt concentration gradients from equilibrating across electrode thickness.

High charging rates applied under subzero transport conditions rapidly deplete lithium salt concentration inside deeper electrode void spaces. As charge current draws lithium ions out of the liquid phase near active surfaces, slow ionic migration fails to replenish salt concentration deep within internal pores. Severe salt depletion reduces local conductivity further, shifting electrical current pathways to the separator surface where lithium ion concentration remains higher.

The resulting current constriction forces localized overpotentials past the thermodynamic threshold for metallic deposition.

Local Overpotential Distribution across Electrode Thickness

Near the separator face, rapid ionic reaction consumes available salt faster than bulk migration replaces it. Spatial modeling and electrochemical impedance measurements confirm that local negative electrode potential varies across electrode depth during low-temperature charge events. Active material thickness directly controls the severity of this spatial overpotential gradient; thick electrodes with high areal capacities suffer far greater spatial potential variations than thin, high-power electrode designs.

To verify the physical onset of pore-level spatial plating during qualification testing, engineering teams execute destructively sampled teardown sequences across specific charge states.

1. Stabilize the test cell at the target subzero operating temperature inside an environmental chamber for six hours to ensure complete internal thermal equilibrium.
2. Apply the specified fast-charge protocol while continuously monitoring cell terminal voltage and surface temperature logging channels.
3. Terminate charging immediately upon reaching eighty percent state of charge or hitting the upper cut-off voltage limit.
4. Transfer the charged cell directly into an argon-filled glove box possessing moisture and oxygen levels below 0.1 ppm to prevent ambient reaction.
5. Disassemble the casing hardware carefully to extract the negative electrode strip without causing mechanical shear stress on active coating surfaces.
6. Rinse the harvested electrode in ultra-pure dimethyl carbonate solvent to remove residual electrolyte salts and dry under vacuum at ambient box temperature.
7. Inspect the active graphite surface under optical magnification and scanning electron microscopy to map spatial distributions of metallic lithium dendrites relative to separator contact zones.

UN 38.3 T.7 test certification remains valid only when subzero charge current limits specified in cell technical datasheets are strictly enforced by the battery management system.

Field failures are frequently attributed to customer charging abuse based on subzero operational logs, leading to denied warranty compensation.

Reversibility

Metallic deposits formed on graphite surfaces follow distinct chemical paths during discharge. Freshly plated lithium metal reacts rapidly with surrounding organic liquid carbonate solvents and dissolved salt anions. This parasitic reaction consumes active lithium ions and liquid electrolyte, generating a secondary solid-electrolyte interphase shell around metallic dendrites.

When the cell transitions from charge to discharge, a fraction of the plated lithium oxidizes back into mobile lithium ions, while another fraction becomes isolated or chemically passivated.

Reversibility defines the proportion of metallic lithium that can be electrochemically stripped and returned to the positive electrode during subsequent discharge steps. Reversible metallic lithium stripping occurs at potentials near 0 V vs Li/Li+, producing a recognizable flat voltage plateau during early discharge phases. Irreversible lithium loss occurs when metallic structures lose electronic contact with the underlying conductive graphite matrix, forming electrochemically isolated regions known as dead lithium.

Why Metallic Depositions Persist after Thermal Recovery?

Isolated metallic dendrites lose electrical contact with the conductive carbon matrix during discharge cycles. As stripping dissolves metal from the base of a dendritic structure where it attaches to graphite, the structural neck breaks before the tip fully oxidizes. The severed metal particle remains physically trapped inside the porous SEI film or separator matrix but lacks an electronic conduction path to transfer electrons to the current collector foil.

Dead lithium cannot participate in subsequent charge or discharge cycles, leading to permanent capacity decay. Furthermore, continuous chemical reaction between fresh metallic lithium and liquid solvent depletes active electrolyte volume and increases internal cell impedance. When a cell subjected to subzero lithium plating is warmed back to 25°C, the electrochemically dead lithium remains unreactive, leaving the cell with permanent capacity loss and elevated internal resistance that persists for its remaining operating lifespan.

Plating severity increases whenever charge transfer resistance at the anode interface exceeds solid-state diffusion impedance within the graphite bulk.

Electrochemical Diagnostics and Non Destructive Signals

Differential capacity analysis identifies high-voltage discharge relaxation plateaus created by metallic stripping. When a cell containing plated metallic lithium begins discharging at ambient temperatures, the stripping of metallic lithium occurs at a lower potential than lithium deintercalation from graphite. This dual-voltage mechanism introduces a distinct secondary inflection peak in differential capacity (dQ/dV) curves during early discharge phases.

Tracking peak magnitude and position allows non-destructive quantification of plated lithium mass without opening the cell casing.

Electrochemical impedance spectroscopy provides additional non-destructive detection capability. The growth of high-surface-area metallic dendrites and the concurrent breakdown and reformation of SEI layers alter high-frequency and mid-frequency Nyquist arcs. Increased charge transfer resistance combined with shifted solid-state Warburg diffusion slopes signals structural surface degradation before micro-short circuits cause catastrophic thermal failure.

• Differential Capacity Voltage Relaxation Tracking mathematical transformation of discharge curves revealing secondary voltage plateaus associated with metallic lithium stripping oxidation.
• High Frequency Impedance Arc Expansion real-time monitoring of charge transfer semicircle growth caused by parasitic SEI formation on metallic deposits.
• Open Circuit Voltage Relaxation Monitoring extended post-charge voltage settling observation detecting self-discharge kinetics driven by chemical lithium re-intercalation.
• Non Destructive Coulombic Efficiency Drop precision coulombic efficiency tracking detecting minor operational deficits caused by parasitic lithium consumption.

The precise thermodynamic threshold where isolated metallic lithium transitions from passive dead metal to active self-heating reactant remains an open question in low-temperature electrochemistry research.

Margin

Cell design specifications balancing low-temperature charge capability against room-temperature energy density accept explicit operational trade-offs. To prevent negative anode potentials during cold charging, manufacturers modify electrode geometry, material chemistry, and mass ratios. Lowering electrode coating thickness reduces ionic path length, suppressing transport overpotentials at the cost of volumetric energy density.

Selecting smaller graphite particle diameters accelerates solid-state diffusion rates while increasing surface area exposed to parasitic SEI reactions.

Negative to positive capacity ratios (N/P ratios) establish physical buffering margins against plating. A cell designed with an N/P ratio of 1.20 contains twenty percent excess negative active graphite capacity relative to positive cathode capacity at 25°C. Under subzero conditions, kinetic transport restrictions effectively decrease accessible graphite capacity. If kinetic limitations reduce effective graphite intercalation capacity by thirty percent, the functional N/P ratio drops below 1.0, making metallic lithium plating thermodynamically inevitable even before accounting for overpotential shifts.

Electrode Architecture Modifications for Low Temperature Charge

Engineering particle size distributions reduces diffusion distances to minimize internal concentration gradients. Surface modification of graphite particles using amorphous carbon coatings lowers charge transfer resistance and reduces solvent desolvation energy barriers. Blending small fractions of silicon or metal alloys into graphite matrices provides fast-intercalation host sites that alter local potential responses during subzero current spikes.

Electrolyte formulation plays an equally significant role in extending low-temperature charging margins. Replacing standard linear carbonates with low-freezing-point solvents like ethyl acetate or fluoroethylene carbonate (FEC) maintains liquid-phase ionic conductivity down to -40°C. Additives that construct thin, highly conductive inorganic SEI films dominated by lithium fluoride (LiF) lower interfacial charge transfer resistance, suppressing overpotentials during low-temperature fast charging.

Operational Charging Envelopes and BMS Safeguards

Battery management systems implement multi-dimensional lookup tables that throttle current based on real-time cell temperature and state of charge. Standard single-value charge current limits fail to account for non-linear polarization growth at temperatures below 0°C. Sophisticated management algorithms calculate instantaneous local anode potential using electrochemical models, adjusting current continuously to maintain the calculated potential above 0.00 V versus Li/Li+.

Pre-heating strategies provide physical temperature elevation prior to applying fast-charge current streams. Internal AC pulse heating and integrated PTC heating foils raise core cell temperatures from -20°C to safe operating zones within minutes. Pre-heating consumes a small fraction of input energy but prevents permanent capacity loss, preserving cycle life and preventing dendritic growth inside separator pores.

Cell datasheets quoting room temperature cycle life offer no indication of subzero durability when charging occurs below freezing.

Incorporating IEC 62660-3 Clause 6.4 subzero charge limits into cell procurement agreements forces suppliers to replace nominal C-rate ratings with explicit temperature-dependent degradation guarantees.