Electrostatic Screening in Nanotubes: A Tubular Response Function Framework

TL;DR

This work develops tubular response functions as a cylindrical analogue to planar surface response for evaluating electrostatic interactions in nanotubes. By deriving these functions from linear response and applying them to dielectric, 1D electron gas, and metallic armchair CNTs (via Luttinger theory), the authors show that metallic walls enforce exponential screening of long-range Coulomb interactions, a property that persists even in realistic CNTs due to circumferential quantum confinement. They also quantify how strong dielectric anisotropy of nanoconfined water enhances short- and intermediate-range interactions and compute the ion self-energy contributions arising from both dielectric background changes and wall polarization. The framework provides a versatile, quantitative tool for predicting ionic correlations, charge storage, and dynamics in nanotube-based nanofluidic devices and can be extended to other confined electrolytes and rolled-up 2D materials.

Abstract

The structure and transport of electrolytes in nanoscale channels are known to be affected by the electronic properties of the confining walls. This influence is particularly pronounced in quasi-one-dimensional nanotubes, where the high surface-to-volume ratio makes the wall the dominant source of electrostatic screening. For instance, ideal metallic tubes suppress long-range Coulomb interactions between ions exponentially. Yet, there exists no generic framework for evaluating electrostatic interactions in tubular confinement. Here, we introduce tubular response functions - a generalisation of surface response functions that captures how nanotubes with arbitrary electronic properties screen Coulomb interactions. Using this framework, we evaluate the interaction potential of ions confined in a metallic carbon nanotube, treating its electronic properties exactly within a Luttinger liquid model. We demonstrate that the long-range exponential screening characteristic of ideal metals persists in realistic metallic nanotubes, regardless of their electron density. We trace the origin of this perfect screening property to the quantum confinement of electrons along the tube circumference. Our framework opens the way for quantitative descriptions of ionic correlations and charge storage in nanotube-based electrodes, and can be further extended to address confined ion dynamics.

Figures (4)

• Figure 1: Water fluctuations in quasi-1D confinement. a Schematic illustration of an ion in a (6,6) armchair carbon nanotube. b Schematic of a nanopore of radius $R$ embedded into an infinite dielectric medium. The radial, azimuthal and axial dielectric constants, $\epsilon_r$, $\epsilon_{\varphi}$ and $\epsilon_z$, are indicated with arrows. c Illustration of the tubular response function $g_s$ as a reflection coefficient for the external potential $\phi^{\text{ext}}$ applied to the solid from the liquid side. The external potential polarises the solid, leading to an induced potential $\phi^{\text{ind}}$ in the liquid. At the solid-liquid interface, the induced potential is $\phi^{\text{ind}}(R) = -g^s \phi^{\text{ext}}(R)$.
• Figure 2: Band structures of a quasi-1D electron gas and a carbon nanotube. a Band structure of a quasi-1D electron gas on an infinite cylinder with radius $R=4\angstrom$. b Band structure of a (6,6) armchair CNT with radius $R \approx 4\angstrom$. a,b For the chosen electron density $n=5e13cm^{-2}$, only the lowest bands are occupied in both cases.
• Figure 3: Tubular response function and screened ion potential. a Tubular response function $g_{0,q}$ for various nanotube systems, defined in the main text, with radius $R=4\angstrom$. The dielectric medium is isotropic with dielectric constants $\epsilon^s_{r,\varphi,z} = 4$. The perfect metal has, by definition, $g_{0,q}=1$. All other models are based on the same electron density $n=5e13cm^{-2}$. The armchair CNT has the chiral indices $(6,6)$ and the radius $R \approx 4.07\angstrom$. b Bare and screened Coulomb potential of an ion at the centre of various nanotubes. The potential is computed along the tube axis with the same parameters as for panel a. The integral of Eq. \ref{['eq:total-ion-potential']} is computed on a logarithmically spaced set of $q$-values between $1e-8\angstrom^{-1}$ and $1e3\angstrom^{-1}$. The inset magnifies the data for small distances $z$. The screened Coulomb potentials of the ideal metal, the quasi-1D electron gas (1DEG) and the metallic armchair CNT fall nearly on top of each other. c Same data as panel b, but in double-logarithmic scale. The slight curvature of the bare potential $\phi^{\text{ext}}$, and consequently all screened potentials, at small $\abs{z}$ is an artifact of the finite radius $\rho=0.03\angstrom$, which regularises the computation of the bare Coulomb potential.
• Figure 4: Self-energy of a unit-charge ion with radius $R_0=1\angstrom$ at the centre of a polarisable nanotube. a Self-energy $\Sigma_{\epsilon}$ due to the change in dielectric background between bulk and nanoconfined water, as a function of the anisotropy factor $\alpha = \sqrt{\epsilon_r/\epsilon_z}$. b Nanotube polarisation contribution $\Sigma_{\text{NT}}$ to the self-energy of an ion at the centre of various nanotubes with radius $R=4\angstrom$ and identical parameters as in Fig. \ref{['fig3']}. c Same as b, here depending on the tube radius $R$. The data for a perfect metal, a 1DEG, a 2DEG and an armchair CNT are almost identical.