Effect of Nanopore Wall Geometry on Electrical Double-Layer Charging Dynamics

Confinement strongly influences electrochemical systems, where structural control has enabled advances in nanofluidics, sensing, and energy storage. In electric double-layer capacitors (EDLCs), or supercapacitors, energy density is governed by the accessible surface area of porous electrodes. Continuum models, built on first-principles transport equations, have provided critical insight into electrolyte dynamics under confinement but have largely focused on pores with straight walls. In such geometries, a fundamental trade-off emerges: wider pores charge faster but store less energy, while narrower pores store more charge but charge slowly. Here, we apply perturbation analysis to the Poisson-Nernst-Planck (PNP) equations for a single pore of gradually varying radius, focusing on the small potential and slender aspect ratio regime. Our analysis reveals that sloped pore walls induce an additional ionic flux, enabling simultaneous acceleration of charging and enhancement of charge storage. The theoretical predictions closely agree with direct numerical simulations while reducing computational cost by 5-6 orders of magnitude. We further propose a modified effective circuit representation that captures geometric variation along the pore and demonstrate how the framework can be integrated into pore-network models. This work establishes a scalable approach to link pore geometry with double-layer dynamics and offers new design principles for optimizing supercapacitor performance.

💡 Research Summary

This paper investigates how the geometry of nanopores influences the charging dynamics of electric double‑layer capacitors (EDLCs), a key factor in the energy density of supercapacitors. While previous continuum studies have focused on pores with straight walls—cylindrical or slit‑shaped—showing a trade‑off where wider pores charge quickly but store less charge and narrower pores store more charge but charge slowly, the authors ask whether a gradual change in pore radius can break this limitation.

To answer this, they consider a single axisymmetric pore whose radius α(z) varies smoothly along its axial coordinate z. The electrolyte is modeled as a symmetric 1:1 binary electrolyte with equal diffusivities for cations and anions. The pore walls are assumed ideally blocking (no Faradaic reactions), and the pore length ℓp is much larger than its characteristic radius ap, giving a slender‑aspect‑ratio parameter δ = ap²/ℓp² ≪ 1. The system is also subjected to a small wall potential Φw (dimensionless voltage) such that Φw ≪ 1, allowing linearization of the governing equations.

The governing equations are the Poisson‑Nernst‑Planck (PNP) set: continuity equations for ion concentrations, Nernst‑Planck fluxes, and Poisson’s equation for the electrostatic potential. The authors nondimensionalize the equations using the bulk concentration c0, the thermal voltage kBT/e, and the pore length ℓp. They then perform a two‑parameter regular perturbation expansion in Φw and δ, retaining terms up to O(Φw δ⁰). At leading order (Φw⁰) the solution is trivial: no charge density, uniform concentration, and zero potential.

At first order in Φw, the radial dependence collapses: the combination ρ₁₀ + 2Φ₁₀ is independent of the radial coordinate R, implying that the axial electric field is uniform across the pore cross‑section. The axial dynamics are governed by a diffusion‑type equation that includes an additional term proportional to α dα/dZ · ∂Z(ρ₁₀ + 2Φ₁₀). This term originates from the no‑flux boundary condition on the sloping wall and represents a “geometric flux” induced by the gradual change in pore radius. Physically, when the pore widens (dα/dZ > 0) the geometry drives ions forward, accelerating charging; when the pore narrows (dα/dZ < 0) the same term slows axial ion transport, allowing more charge to accumulate locally. Importantly, this effect appears for any Debye length (any κ = ap/λ), so it is robust across electrolyte concentrations and solvent properties.

The authors validate the analytical predictions against full two‑dimensional numerical solutions of the PNP equations. Across a range of Debye lengths, wall potentials up to Φw ≈ 4 (≈ 100 mV), and several wall‑shape functions α(Z), the analytical charge density and potential profiles match the numerical results to within five to six significant figures.

Beyond the single‑pore theory, the paper proposes a modified transmission‑line circuit model that incorporates spatially varying resistance and capacitance derived from the analytical solution. This “effective circuit” can be embedded in pore‑network models, enabling rapid simulation of thousands of interconnected pores while preserving the essential physics of geometry‑induced fluxes. Computationally, the perturbation‑based approach reduces the cost by five to six orders of magnitude compared with direct PNP simulations, making it feasible for large‑scale electrode design studies.

The key insight is that by engineering the slope of pore walls, one can simultaneously accelerate charging and increase stored charge—effectively breaking the classic trade‑off between pore size, charging speed, and energy density. This design principle opens new avenues for high‑power, high‑energy supercapacitors required in electric vehicles, grid‑scale storage, and fast‑response sensors. The authors suggest future extensions to include nonlinear potentials, ion‑size effects, dielectric decrement, and convective flows, as well as experimental validation of the predicted performance gains.