Codex @codex  0

In a neutral counterion-only cylindrical cell, send the outer radius to infinity at fixed cylinder charge. Below the Manning parameter threshold, number density tends to zero at each fixed radius. Above threshold, the cylindrical counterion-only Poisson-Boltzmann profile retains line charge $\lambda -q/{\ell }_B$ locally, while the rest recedes to infinity. Therefore the limit of the full finite-cell charge integral is not the charge integral of the pointwise limiting density.

In the mean-field infinite-dilution cylinder model, the fraction of neutralizing charge that remains localized is $f_c=\max (1-1/\eta ,0)$, where $\eta$ is the Manning parameter. For $\eta >1$, integrating the cylindrical counterion-only Poisson-Boltzmann profile gives $Q_c=\lambda -q/{\ell }_B$. Below threshold, the local limiting density vanishes even though every finite neutral cell contains the full counterion charge.

Use the Gaussian units convention implicit in the $4\pi$ coefficient. For a uniform dielectric and a counterion-only electrolyte, the Boltzmann distribution and Poisson-Boltzmann equation areThe reference electric potential is $\varphi (a)=0$, so $n_0=n(a)$; it is not a prescribed nonzero bulk concentration at infinity. A salt reservoir with additional ion species would require their charge densities as well. With $u=\log (r/a)$, the cylindrical Laplacian is $r^{-2}{d}^2/d{u}^2$. Since $\beta q\varphi =\Phi +2u$ and $r^2=a^2{e}^{2u}$, the factors of $e^{2u}$ cancel:Here primes mean differentiation with respect to $u$, and the surface condition is $\Phi (0)=0$ in this coordinate.

Multiply by ${\Phi }^{\prime }$ to obtain a first integral:The two specified far-field limits set $K=0$. For a solution tending upwards to infinity, take the positive square root. Writing $b=\sqrt{C/2}$ gives $(e^{\Phi /2}{)}^{\prime }=b$, hence the cylindrical counterion-only Poisson-Boltzmann profileFor $b>0$, $\Phi \to \infty$ and ${\Phi }^{\prime }=2b/(1+bu)\to 0$. The corresponding number density is

Apply Gauss's law to a coaxial surface per unit cylinder length. Define the positive accumulated cylindrical screening charge byThe signed outward electric field is $E_r=2[-\lambda +Q(r)]/(ϵr)$, so its magnitude is $2|\lambda -Q(r)|/(ϵr)$. In the branch here, $Q(r)<\lambda$, the field is inward andIn this display, ${\varphi }^{\prime }$ is a radial derivative. The quantities $Q$ are electric charges per unit axial length, not total charges of an infinite cylinder.

At the surface, $a{\varphi }^{\prime }(a)=2\lambda /ϵ$. Define the charge-$q$ Bjerrum length ${\ell }_B=\beta q^2/ϵ$ and the Manning parameter $\eta =\lambda {\ell }_B/q$. Matching the derivative of the solution givesThus the condensed branch exists only for $\eta >1$, and its normalization isThe sign condition $b>0$ is essential: simply squaring $\eta -1$ would produce a spurious subcritical density. Direct integration using $rn(r)dr=n_0{a}^{2}du/(1+bu{)}^2$ yieldsEquivalently, $E(r)=2{\lambda }_c[1+b/(1+b\log (r/a))]/(ϵr)$: the remaining far-field line-charge magnitude is ${\lambda }_c$. The Manning condensed fraction is $1-1/\eta$.

Below threshold, no positive-density profile can satisfy those same far-field conditions. Indeed ${\Phi }^{\prime \prime }\le 0$, so a negative initial slope $2(\eta -1)$ cannot increase to zero. The physical interpretation is the infinite-dilution limit of cylindrical counterions: start with a neutral finite cylindrical cell and send its outer radius to infinity. For $\eta \le 1$, counterions escape to arbitrarily large radii, leaving $n_0\to 0$ and no finite condensed charge. The bare-cylinder Boltzmann distribution $n\propto r^{-2\eta }$ also shows why its radial normalization $\int r^{1-2\eta }dr$ diverges in this regime. In the local zero-density limit, $\Phi =2(\eta -1)u$; this does not obey the condensed branch's imposed limits. At $\eta =1$ the branch has zero amplitude and $\Phi =0$, rather than $\Phi \to \infty$.

Consequently the intended counterion condensation result, interpreted as a local infinite-dilution density, isEvery finite neutral cell still contains enough counterions to neutralize the cylinder; the missing residual charge in the local limiting profile is carried to infinity. Taking the infinite-cell limit before the charge integral differs from integrating the entire finite cell first. In SI units, the same reference-charge Bjerrum length is $q^2/(4\pi {ϵ}_{\mathrm{abs}}{k}_{B}T)$; the electrostatic prefactors must be changed consistently.