Quantifying Sodium Quasi-Metallic Cluster Nucleation Potentials inside Sub-Nanometer Closed Carbon Pores

Nucleus

Sodium storage within non-graphitizable hard carbon anodes occurs through two distinctly separated thermodynamic regimes during electrochemical discharge. The high-potential sloping region between 1.20 V and 0.10 V versus Na/Na+ originates from ionic adsorption at defective graphene edge sites and surface heteroatom functional groups alongside interstitial insertion between expanded turbostratic layers. The low-potential quasi-flat plateau sitting between 0.10 V and 0.00 V versus Na/Na+ delivers the primary share of reversible capacity through the condensation of metallic-like sodium clusters inside closed sub-nanometer pores.

Accurately decoupling these mechanisms requires evaluating the change in Gibbs free energy associated with pore confinement, where nanometer-scale spatial constraints suppress bulk phase growth while altering the chemical potential of trapped sodium species.

When sodium species enter closed carbon pores with internal diameters below 1.0 nm, spatial confinement forces a transition from solvated ion transport to dense atomic aggregation. The thermodynamic driving force for this process deviates from classical bulk metal deposition described by the Nernst equation. Inside closed pores, capillary condensation dynamics and strong van der Waals interactions between sodium atoms and the surrounding carbon wall stabilize small sodium clusters at potentials slightly above zero volts relative to reference bulk sodium metal.

The shift in nucleation potential is governed by the modified Gibbs-Thomson relationship, expressed as:

ΔV = (2 γ V_m) / (F r_p)

where γ represents the interfacial energy between sodium metal and the carbon surface, V_m represents the molar volume of metallic sodium, F represents the Faraday constant, and r_p represents the effective pore radius. For pore radii between 0.25 nm and 0.40 nm, the calculated nucleation potential shift ranges between 8 mV and 38 mV positive relative to bulk sodium deposition. This positive potential shift establishes a narrow thermodynamic window where quasi-metallic cluster formation proceeds without triggering outer-surface metallic plating.

Thermodynamic Driving Force for Sodium Cluster Formation

Reversible sodium storage at low voltage hinges on maintaining the localized chemical potential of the sodium cluster below the threshold for phase separation into crystalline sodium. Within sub-nanometer closed pores, sodium atoms arrange in non-crystalline geometries containing between 4 and 13 atoms. Electronic structure calculations indicate that these quasi-metallic clusters exhibit electron density distributions intermediate between isolated atomic species and delocalized metallic conduction bands.

The lack of long-range metallic bonding suppresses the formation of continuous Fermi surfaces, reducing localized reactivity with surrounding carbon atoms.

Kinetic energy barriers during the transition from adsorbed ion to quasi-metallic cluster depend heavily on charge transfer dynamics at the pore entrance. As Na+ ions approach an open defect or pore throat, the energetic penalty of shedding the primary solvation sheath, such as ethylene carbonate or propylene carbonate molecules, dictates the local activation energy. If the desolvation energy exceeds the energetic benefit gained through pore confinement, the charge transfer overpotential pushes the local anode potential below 0.00 V versus Na/Na+.

This condition converts harmless interior pore filling into aggressive dendrite nucleation on the exterior carbon surface.

A 0.35 nm closed pore stabilized at 25 °C demonstrates a 22 mV positive potential shift relative to bulk sodium plating during C/10 discharge.

Quantifying Critical Nucleation Radius in Pore Confinement

The transition from isolated atomic filling to collective cluster nucleation requires a critical density of sodium atoms within the closed pore volume. Classical nucleation theory dictates that the critical radius for stable phase formation decreases as cathodic overpotential increases. Inside sub-nanometer carbon pores, physical geometric boundaries replace the thermodynamic critical radius when the pore radius falls below the classical critical nucleus size.

The pore wall acts as a physical boundary, restricting cluster growth to the internal volume of the cavity and preventing classic spherical growth geometries.

Determining the optimal pore dimension involves balancing gravimetric storage density against kinetic impedance. Closed pores smaller than 0.30 nm exhibit strong steric repulsion that prevents multi-atom cluster organization, leaving interior pore capacity unused. Pores larger than 1.20 nm lack sufficient van der Waals interaction from opposing pore walls, causing the nucleation potential to collapse down to 0.00 V versus Na/Na+ and matching the bulk plating potential.

Hard carbon materials synthesized to maximize closed pores in the 0.50 nm to 0.85 nm range achieve maximum plateau capacity while providing a safety buffer against exterior dendrite formation.

Selecting hard carbon specifications with pore distribution profiles skewed outside the 0.45 nm to 0.99 nm window drives plateau capacity loss, elevates initial coulombic inefficiency, and causes cell thermal runaway during low-temperature fast charging.

Slit

Carbon microstructures in non-graphitizable hard carbons consist of random, highly disordered turbostratic domains where short-range graphitic crystallites intertwine to form an interconnected network of closed voids and narrow slit-like channels. The spatial arrangement of these graphene sheets determines both the volumetric density of the active material and the mass transport pathways available to mobile sodium ions. While parallel graphene layers separated by d002 interlayer spacings above 0.37 nm facilitate fast interstitial sodium diffusion during the sloping voltage region, the entrance throat dimensions of the closed pores dictate whether sodium clusters can nucleate within the interior volume during the low-voltage plateau.

Slit-pore geometry differs fundamentally from spherical void models in its directional transport constraints and physical aspect ratios. Slit openings are bounded by the basal planes and defect-rich edges of adjacent turbostratic carbon stacks. The effective width of these channels determines the steric hindrance experienced by desolvated sodium species entering the internal cavity.

High-resolution transmission electron microscopy and automated image processing demonstrate that optimal hard carbon architectures exhibit curved, onion-like graphitic domain walls that seal off internal void space while maintaining ultra-narrow access necks.

Steric Constraints and Desolvation Energetics at Pore Throats

The solvation structure of the sodium ion in conventional organic ester electrolytes comprises four to six solvent molecules forming a primary coordination shell with a hydrodynamic diameter between 0.70 nm and 0.90 nm. Entering a sub-nanometer carbon pore neck measuring less than 0.60 nm across requires complete shedding of this coordination sphere at the carbon-electrolyte interface. The activation energy for this stripping process contributes directly to the total charge transfer resistance of the anode.

When carbon pore necks fall below 0.33 nm, the bare sodium ion, possessing an ionic radius of 0.102 nm, encounters electrostatic Pauli repulsion from the dense pi-electron clouds projecting off the graphitic pore walls. This energetic barrier restricts ion entry, effectively rendering the internal closed volume electrochemically inactive. Tuning the defect concentration at pore entrances allows engineered hard carbon materials to accelerate desolvation kinetics without enlarging the interior pore throat to the point where solvent molecules infiltrate and decompose inside the cavity.

Structural Metrics Governing Interlayer Spacing and Closed Void Volume

X-ray diffraction analysis provides quantitative metrics regarding the average crystallite dimensions within hard carbon matrices. The interlayer spacing, calculated from the (002) diffraction peak positioning using Bragg’s law, ranges between 0.365 nm and 0.400 nm in high-performance sodium-ion battery carbons. Values below 0.360 nm limit interstitial insertion capacity, whereas values exceeding 0.410 nm indicate excessively disordered, low-density structures susceptible to high initial capacity losses.

Beyond interlayer spacing, the average stack height (Lc) and lateral crystallite size (La) derived from Scherrer formulations quantify the degree of structural order. Closed void volume correlates inversely with Lc and La values; smaller, highly contorted graphitic domains yield a higher spatial frequency of enclosed nanovoids upon thermal consolidation. Precision density measurements combining helium pycnometry with skeletal density extraction reveal that optimal hard carbon materials contain a closed pore volume between 0.12 cm³/g and 0.22 cm³/g, supporting theoretical plateau capacities exceeding 250 mAh/g.

• Pore Throat Constriction limits sodium ion entrance when localized structural defects collapse during high-temperature thermal treatment steps, blocking access to internal closed cavities.
• Solvent Molecule Infiltration occurs when pore throat diameters exceed 0.85 nm, leading to continuous electrolyte reduction inside closed voids and permanent capacity loss.
• Graphitic Crystallite Coalescence reduces closed pore volume through excessive parallel alignment of graphene layers during prolonged carbonization profiles.
• Surface Oxygen Trapping at slit openings generates irreversible sodium-oxygen complexes that increase interfacial impedance and decrease first-cycle coulombic efficiency.

Closed pores sealed by contorted graphitic sheets resist electrolyte penetration while maintaining high mechanical rigidity against structural swelling during sodium cluster insertion.

Resonance

Differentiating between intercalated sodium, surface-adsorbed ions, metallic sodium dendrites, and quasi-metallic clusters inside sub-nanometer pores requires diagnostic techniques sensitive to localized electronic environments. Nuclear magnetic resonance spectroscopy provides non-destructive, isotope-specific analysis of sodium environments through 23Na chemical shift measurements. Because 23Na possesses a nuclear spin of 3/2, its resonance spectrum reflects both localized magnetic shielding and quadrupolar interactions driven by electric field gradients in the immediate coordination environment.

In hard carbon anodes subjected to electrochemical sodium loading, distinct chemical shift regimes correspond to specific storage states. Ionic sodium adsorbed on surface functional groups or residing within liquid electrolyte appears in a narrow chemical shift range between -20 ppm and +10 ppm relative to standard aqueous sodium chloride solutions. Intercalated sodium positioned between expanded turbostratic carbon layers yields resonance signals shifted downfield to between +15 ppm and +50 ppm.

Metallic sodium dendrites deposited on outer particle surfaces display a pronounced Knight shift between +1100 ppm and +1150 ppm, caused by conduction electron spin susceptibility associated with bulk metallic states.

Nuclear Magnetic Resonance Knight Shifts as Electronic Identifiers

Quasi-metallic sodium clusters confined inside sub-nanometer closed carbon pores produce an intermediate 23Na NMR shift profile located between +100 ppm and +750 ppm. The magnitude of this Knight shift correlates directly with cluster size and localized electron density delocalization. As sodium atoms aggregate inside a closed void during low-potential charging, the 23Na resonance peak continuously shifts downfield from +120 ppm toward +650 ppm, reflecting the emergence of quasi-metallic electronic character within the confined atomic grouping.

Variable-temperature 23Na solid-state NMR experiments allow calculation of the electronic density of states at the Fermi level for these confined clusters. Unlike bulk sodium metal, which maintains a temperature-independent Knight shift down to cryogenic temperatures (Pauli paramagnetism), quasi-metallic clusters show temperature-dependent shifts below 150 K. This temperature dependence confirms quantum confinement effects, proving that the sodium species exist as discrete nanometer-scale aggregates rather than continuous metallic phases.

Solid-state 23Na NMR spectra recorded at 20 °C resolve a discrete quasi-metallic peak at +280 ppm, separated from the bulk metallic signal by 850 ppm.

Scattering Techniques for Closed Pore Volume Extraction

Small-angle X-ray scattering and small-angle neutron scattering offer non-invasive characterization of enclosed pore geometries inaccessible to conventional gas adsorption techniques like nitrogen or argon physisorption. Gas molecules cannot penetrate true closed pores, leading gas physisorption to drastically underestimate total porosity. Scattering techniques exploit the electron density contrast between the solid carbon matrix and the empty or sodium-filled internal cavities.

Analyzing the scattering intensity vector I(q) across the low-angle range (q = 0.01 Å⁻¹ to 0.5 Å⁻¹) using Guinier and Porod approximations yields precise pore size distribution curves. Contrast matching experiments, wherein fluid media with matched scattering length densities fill only open surface defects, isolate scattering signals originating purely from closed internal voids. Operando SAXS measurement performed during cell discharge reveals changes in scattering contrast as sodium clusters nucleate inside closed pores, enabling real-time tracking of internal pore filling kinetics as a function of applied state of charge.

Calculations derived solely from skeletal density can suggest high closed pore volumes while obscuring broad pore size distributions that include macro-voids incapable of stabilizing quasi-metallic clusters.

Hysteresis

Voltage polarization and kinetic hysteresis between charge and discharge cycles present engineering challenges in commercial hard carbon sodium-ion cells. While the high-potential sloping region exhibits minimal potential separation between insertion and extraction pathways, the low-potential plateau displays kinetic voltage hysteresis that broadens under elevated current rates. This hysteresis arises from structural restructuring of sodium clusters during extraction, charge-transfer resistance across desolvated pore throats, and localized concentration polarization inside narrow carbon channels.

During low-potential charging at high C-rates, ohmic drop and mass transport limitations sum with the thermodynamic equilibrium potential. This combined overpotential can drive the instantaneous surface potential of the hard carbon particle below 0.00 V versus Na/Na+ even when the interior closed pores remain partially unfilled. When local surface potential drops below 0.00 V, the thermodynamic barrier for bulk metallic plating collapses, initiating dendrite growth on outer grain surfaces and compromising cell operating safety.

Overpotential Drift during High C-Rate Low-Potential Charging

Quantifying total overpotential during the plateau region requires isolating individual resistance contributions through electrochemical impedance spectroscopy and galvanostatic intermittent titration techniques. The total cell overpotential (η_total) represents the sum of ohmic drop (η_Ω), solid electrolyte interphase resistance (η_SEI), charge transfer overpotential (η_ct), and solid-state pore diffusion overpotential (η_diff):

η_total = I R_Ω + η_SEI + η_ct + η_diff

As current density increases from C/20 to 2C, the solid-state diffusion overpotential inside sub-nanometer carbon channels increases exponentially due to spatial congestion along single-file transport pathways. This spatial congestion increases charge-transfer resistance at pore throats, forcing a voltage drop that terminates low-potential plateau charging prematurely and reduces usable energy density.

Can Sub-Nanometer Pores Suppress Dendrite Formation at Sub-Zero Voltages?

Operating a hard carbon anode at negative potentials relative to bulk sodium reference electrodes risks immediate plating if surface defects act as unpassivated nucleation sites. Sub-nanometer closed pores mitigate this risk by providing localized energy wells that draw sodium inward via strong capillary forces. However, this protective capillary action functions only when the current density remains below the critical diffusion limit of the carbon pore neck network.

Exceeding the critical diffusion current density causes sodium ions to accumulate at the carbon particle boundary, building an electric double-layer field that favors surface plating over internal pore diffusion. Fast-charging protocols must monitor differential voltage curves (dQ/dV) to detect the onset of sharp, low-voltage stripping peaks during subsequent discharge steps, which indicate unwanted surface metal plating.

1. Mount the target hard carbon electrode inside a temperature-controlled three-electrode test cell utilizing a sodium metal reference electrode and an ultra-pure non-aqueous ester electrolyte.
2. Perform three preliminary conditioning cycles between 2.00 V and 0.01 V versus Na/Na+ at C/10 to establish a stable solid electrolyte interphase layer.
3. Apply a multi-step galvanostatic discharge profile, stepping current density progressively from C/20 up to 3C while monitoring the precise potential transition into the low-voltage plateau.
4. Extract solid-state chemical diffusion coefficients using galvanostatic intermittent titration by recording potential relaxation curves following 10-minute current pulses across the plateau region.
5. Map differential capacity curves across temperature steps ranging from -20 °C to +55 °C to identify the temperature dependence of nucleation overpotentials.

Standard supply agreements specify that hard carbon anodes operated at 1C rates down to 5 mV vs Na/Na+ must retain 92% plateau capacity retention after 1000 cycles at 25 °C.

Standard master procurement contracts include strict limits on maximum allowable voltage hysteresis during C/2 discharge, where exceeding a 35 mV hysteresis threshold grants the buyer full rights to reject the delivered material lot.

Pyrolysis

Synthesis of non-graphitizable hard carbon with tailored sub-nanometer closed porosity requires careful selection of organic precursors and tight control over high-temperature thermal treatment profiles. Natural biopolymers such as sucrose, starch, cellulose, and coconut shells, along with synthetic polymers including phenolic resins and coal-tar pitches, serve as standard carbon feedstocks. The chemical structure of the precursor, specifically its oxygen-to-carbon ratio and crosslinking density, governs whether the material undergoes graphitization or consolidates into a disordered turbostratic framework during thermal decomposition.

Biomass precursors require extensive pre-treatment and washing protocols to eliminate inorganic impurity elements, such as potassium, sodium, silicon, and iron, which catalyze localized graphitization and form unwanted open mesopores. Synthetic precursors, while offering higher chemical purity and structural batch-to-batch consistency, carry higher raw material costs that impact the final landed cost per kilowatt-hour of cell storage capacity.

Thermal Processing Windows for Optimal Closed Pore Generation

The carbonization thermal profile splits into two critical temperature regimes: low-temperature stabilization (200 °C to 500 °C) and high-temperature final annealing (1000 °C to 1600 °C). During stabilization, oxidative crosslinking introduces oxygen-containing functional groups that anchor aromatic domains, preventing the alignment of carbon sheets into parallel graphitic stacks. Inadequate crosslinking during this phase yields high graphitic order, reducing internal closed pore volume.

High-temperature annealing dictates the ultimate pore closure mechanism. As carbonization temperatures ascend from 1000 °C to 1300 °C, volatile heteroatoms like hydrogen and oxygen evolve, creating an abundance of temporary open micropores. Raising the temperature further into the 1300 °C to 1500 °C window induces local structural reorganization where flexible graphene sheets curve, crosslink, and collapse at their edges, sealing open micropores into sub-nanometer closed internal cavities.

Thermal treatment exceeding 1550 °C leads to excessive pore collapse, shrinking internal void volume and drastically reducing plateau storage capacity.

Quality Control Parameters for Incoming Hard Carbon Shipments

Industrial qualification of incoming hard carbon shipments involves rigorous verification of physical, structural, and electrochemical properties to ensure compliance with cell manufacturing tolerances. Surface area determination via nitrogen gas physisorption using Brunauer-Emmett-Teller analysis serves as a primary screening metric. Standard battery-grade hard carbon materials must maintain a BET surface area below 5.0 m²/g, with premium grades achieving values below 2.0 m²/g.

Low surface area minimizes irreversible electrolyte reduction and solid electrolyte interphase formation during the initial formation cycle.

Particle size distribution, measured via laser diffraction, directly affects slurry rheology and electrode coating consistency. Parameters D10, D50, and D90 must be controlled within tight windows to prevent particle segregation during slurry mixing and slot-die coating. Particle morphology, spherical versus irregular morphology, influences tap density, which dictates electrode packing density and volumetric cell energy density.

• BET Surface Area Screening rejects carbon lots displaying surface areas above 5.0 m²/g to prevent excessive first-cycle capacity loss during solid electrolyte interphase formation.
• Tap Density Verification ensures incoming carbon lots exceed a minimum threshold of 0.80 g/cm³ for maintaining packed volumetric energy standards.
• X-Ray Diffraction Profiling confirms that the (002) interlayer spacing remains strictly within the target range of 0.370 nm to 0.390 nm.
• Ash Content Analysis limits inorganic mineral impurities to below 200 ppm via thermogravimetric burning under ambient air conditions.
• Slurry Viscosity Consistency Checks track particle surface chemistry variations that disrupt solvent-binder dispersion dynamics during electrode mixing.

A hard carbon lot carbonized at 1400 °C achieving a BET surface area of 1.4 m²/g yields an initial coulombic efficiency of 88.5% in standard ester electrolytes.

How do chemical surface etching treatments using mild carbon dioxide gas exposures at 800 °C alter entrance throat diameters without damaging internal closed void volumes?

Parity

Sodium-ion battery cell technology achieves market penetration by presenting a lower landed-cost alternative to lithium iron phosphate (LFP) systems in stationary energy storage and low-speed electric transportation segments. Because sodium cathode materials like sodium transition metal layered oxides and polyanionic compounds exhibit lower gravimetric energy densities than lithium equivalents, hard carbon anode performance governs cell-level parity. The volumetric capacity ceiling of hard carbon, tied directly to its skeletal density and closed pore volume, establishes the ultimate physical boundary for sodium-ion energy density.

Commercial sodium-ion prismatic cells using hard carbon anodes currently achieve gravimetric energy densities between 150 Wh/kg and 170 Wh/kg, compared to 180 Wh/kg to 210 Wh/kg for modern LFP cells. Volumetric energy density trails further behind, with sodium-ion cell designs reaching 280 Wh/L to 350 Wh/L against LFP’s 380 Wh/L to 450 Wh/L. This volumetric penalty stems directly from hard carbon’s lower skeletal density (1.50 g/cm³ to 1.80 g/cm³) compared to synthetic graphite (2.20 g/cm³ to 2.26 g/cm³).

Volumetric Density Penalties and Cathode Active Material Matching

To optimize pack-level spatial utilization, cell designers must carefully match hard carbon anode loading with cathode capacities. Initial coulombic efficiency (ICE) mismatches represent a primary source of capacity loss. Hard carbon materials typically demonstrate ICE values between 83% and 90%, whereas layered oxide cathodes yield ICE values between 88% and 94%.

This disparity requires over-balancing the cathode loading or introducing sacrificial sodium additives during electrode slurry preparation to offset first-cycle sodium loss.

Pre-metallation techniques, such as applying sodium metal powder or chemical pre-sodium agents to the hard carbon electrode, elevate first-cycle coulombic efficiency above 96%. However, these processing steps add capital expenditure and manufacturing complexity, partially offsetting raw material cost savings. Cell manufacturers must evaluate whether the volumetric capacity gain achieved through pre-sodium processes balances the added processing cost per kilowatt-hour produced.

Landed Cost Models and Levelized Cost of Storage Calculations

Levelized cost of storage (LCOS) models capture the combined operational economic impacts of round-trip efficiency, cycle life, calendar degradation, capital cost, and thermal management burdens. Despite a lower volumetric density, sodium-ion batteries employing hard carbon anodes achieve favorable LCOS metrics in stationary energy storage applications due to abundance of raw sodium materials, elimination of copper current collectors on the anode side (aluminum foil functions reliably at sodium potentials without alloy formation), and low low-temperature heating requirements.

Replacing copper current collectors with aluminum foils on the anode reduces bill-of-materials costs by approximately $4.00 to $7.00 per kilowatt-hour. Combined with precursor cost advantages, sodium-ion cells present a path toward cell manufacturing costs under $50.00/kWh at gigawatt-scale production volumes. Commercial execution depends on securing high-yield hard carbon supply agreements that guarantee tight structural specifications across multi-ton production runs.

Landed cost calculations must include dangerous goods freight classifications under UN 3480 and UN 38.3 test requirements, where sodium-ion packs transported at zero percent state of charge qualify for relaxed thermal hazard provisions during air and maritime transit.