Period 3 power,
built for the grid.

01 / Periodic table origins

Why moving from Period 2 (Li, Z=3) to Period 3 (Na, Z=11) changes everything.

■ Sodium-Ion (SIB: Prussian White, Layered Oxide, Polyanionic NFPP)■ Lithium-Ion (LIB: LFP, LMFP, NMC-811)● White marker = Cell manufacturing cost floor ($/kWh)

Every engineering trade-off between sodium-ion (SIB) and lithium-ion (LIB) batteries follows directly from adding a closed \(n=2\) neon core shell \((1s^2 2s^2 2p^6)\) when moving down Group 1 from Lithium \(([\text{He}]2s^1,\; M_{\text{Li}}=6.94\text{ g/mol})\) to Sodium \(([\text{Ne}]3s^1,\; M_{\text{Na}}=22.99\text{ g/mol})\):

• Shannon Octahedral Radius & Hard Carbon Intercalation: The 6-coordinate ionic radius expands from \(r_{\text{Li}^+} = 0.76\text{ \AA}\) to \(r_{\text{Na}^+} = 1.02\text{ \AA}\) (a \((1.02/0.76)^3 \approx 2.42\times\) ionic volume expansion). Because \(\text{Na}^+\)-graphite binary staging compounds are thermodynamically unstable above ~35 mAh/g due to π-d orbital repulsion, SIBs replace graphite with disordered hard carbon (turbostratic nanodomains + closed nanopores delivering 300–350 mAh/g).
• Aluminum Foil on Both Electrodes & 0 V Transport Safety: Unlike \(\text{Li}\), which alloys with Aluminum at low potentials (~0.3 V vs. \(\text{Li/Li}^+\)) and forces LIBs to use expensive, heavy Copper foil on the anode, \(\text{Na}\) forms no binary alloy with \(\text{Al}\). SIBs use inexpensive Aluminum current collectors on both cathode and anode — cutting current-collector cost by ~70% and allowing cells to be shipped and stored completely shorted at 0.0 V with zero copper dissolution or thermal runaway hazard.
• Standard Potential & Energy Density Ceiling: By Faraday's law \(Q_{\text{th}} = nF / (3.6 M)\), metal specific capacity drops from 3,861 mAh/g (Li) to 1,166 mAh/g (Na) and \(E^\circ\) shifts from \(-3.04\text{ V}\) to \(-2.71\text{ V}\) vs. SHE, placing practical SIB cells at 135–165 Wh/kg versus 190 Wh/kg for LFP and 275 Wh/kg for NMC-811.

02 / Preventing massive grid blackouts

Why lower Stokes desolvation energy makes Na-ion the ideal polar-vortex and Dunkelflaute BESS.

Why does a larger ion move faster across the electrolyte interface in extreme cold? Because \(\text{Na}^+\) has a larger crystallographic radius (1.02 Å vs. 0.76 Å), its surface electrostatic charge density \(q / (4\pi r^2)\) is 44% lower than \(\text{Li}^+\). That weaker electric field binds carbonate/ether solvent molecules less tightly, reducing the interfacial desolvation activation barrier \(E_a^{\text{desolv}}\) by 25%–30% and yielding a smaller Stokes solvated radius in solution.

During winter polar-vortex blackout emergencies (such as the 2021 Texas ERCOT freeze or Canadian prairie cold snaps at −20 °C to −35 °C), LFP grid batteries suffer severe charge-transfer impedance spikes and lose 36%–72% of their usable capacity unless they burn auxiliary HVAC heating power. Polyanionic (NFPP) and Prussian White sodium-ion grid banks retain 88%–92% capacity at −20 °C, achieve 8,000–10,000 cycles, and eliminate lithium supply bottlenecks for 4-to-12 hour utility peak shaving and synthetic inertia.

"""Reproducible Pandas analysis for ST-05: Sodium Battery Technology vs. Lithium-Ion."""
import numpy as np
import pandas as pd

# 1. Periodic Table & Faraday Law Fundamental Ratios
F = 96485.3321  # C/mol
M_Li, M_Na = 6.94, 22.99  # g/mol
r_Li, r_Na = 0.76, 1.02   # Shannon 6-coord octahedral ionic radius (Angstroms)
E0_Li, E0_Na = -3.04, -2.71  # V vs SHE

q_Li = F / (3.6 * M_Li)  # 3861 mAh/g
q_Na = F / (3.6 * M_Na)  # 1166 mAh/g
vol_ratio = (r_Na / r_Li) ** 3

print(f"Theoretical metal specific capacity: Li = {q_Li:.0f} mAh/g vs. Na = {q_Na:.0f} mAh/g")
print(f"Shannon ionic volume expansion (Na+/Li+)^3: {vol_ratio:.2f}x | Potential shift: +{E0_Na - E0_Li:.2f} V")

# 2. Load 6-chemistry benchmark table & compute Levelized Cost of Storage (LCOS $/MWh-cycle)
df = pd.read_csv("studies-sodium-chemistries.csv")
df["cost_per_1000_cycles_usd_kwh"] = (df["cellCostUsdKwh"] / (df["cycleLife"] / 1000.0)).round(2)
df["winter_effective_wh_kg_minus20C"] = (df["gravWhKg"] * df["tempRetMinus20Pct"] / 100.0).round(1)

print("\n=== Sodium-Ion vs. Lithium-Ion Chemistry Matrix (Pandas) ===")
print(df[["id", "gravWhKg", "cellCostUsdKwh", "cycleLife", "tempRetMinus20Pct", "cost_per_1000_cycles_usd_kwh", "winter_effective_wh_kg_minus20C"]])