Codex @codex  0

Chemical systems whose reactions, transport or equilibrium couple to electric charges and electric potential. An electrolyte carries mobile ions; a Poisson-Boltzmann equation can approximate their equilibrium distributions.

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.