Past exam of the biological physics course of the University of Oxford / 2021 / 4 / d / Solution

Adding $40000$ monovalent counterions makes the combined genome-plus-ion electric charge nominally zero. It does not make the local electric field or free energy vanish: opposite charges have spatial distributions, and confining mobile ions costs entropy.

A definite screened continuation of b replaces the bare Coulomb potential by $e^{-\kappa r}/(4\pi ϵr)$, with $\kappa =1/{\lambda }_D$ and ${\lambda }_D=1\mathrm{nm}$. In the Debye–Hückel approximation, a uniformly distributed fixed charge density $\rho =Q/V_c$ satisfies $({\nabla }^2-{\kappa }^2)\phi =-\rho /ϵ$ inside, with zero fixed charge outside and matching potential and radial field at $R$. Its interior solution isThe screened charging free energy $U_{\kappa }=\frac{1}{2}\int \rho \phi dV$ is thereforeFor $R\gg {\lambda }_D$, the bracket approaches one: only charges within roughly a screening length contribute appreciably, yielding the local estimate $U_{\kappa }\simeq Q^2{\lambda }_D^2/(2ϵV_c)$. With the same ${ϵ}_r=80$ as before, $x=20$ gives bracket $0.92519$ andabout $173$ times smaller than the unscreened packed-sphere self-energy. This is a specified continuum estimate of packed-state charging free energy; a full packaging free energy also subtracts the screened unpacked reference and includes ion redistribution. It should not be presented as a unique exact packaging work fixed by the problem's few inputs.

The supplied counterion count also permits an independent estimate of the entropic contribution. For $N$ ideal confined ions, $F_{\mathrm{ion}}=-N{k}_{B}T\ln V+\text{constant}$, so compression from an accessible reference volume $V_0$ to $V_f$ costsFor example, approximate the original screening atmosphere as an annulus of thickness ${\lambda }_D$ around a rod. With $a=0.5\mathrm{nm}$, $V_0=\pi [(a+{\lambda }_D{)}^2-a^2]L=4.27\times 10^4\mathrm{n}{\mathrm{m}}^3$. Taking $V_f=V_c-V_D=1.35\times 10^4\mathrm{n}{\mathrm{m}}^3$ gives $\Delta F_{\mathrm{ion}}\simeq 1.91\times 10^5\mathrm{pN}\mathrm{nm}$. This useful order-of-magnitude ion-entropy model is distinct from the screened charging model; the two estimates must not simply be added as independent terms without a consistent ionic free energy functional. The rod cutoff, fixed DNA volume and screening-atmosphere geometry are rough geometric assumptions, not a microscopic reconstruction of the source sketch.

There is a further consistency limit. The nominal ion number density is $n=N/V_c=1.19\mathrm{n}{\mathrm{m}}^{-3}$. For a uniform ideal monovalent screening gas, ${\lambda }_D^2=ϵk_{B}T/(n{e}^2)$ would give ${\lambda }_D\simeq 0.309\mathrm{nm}$ at $300\mathrm{K}$, rather than $1\mathrm{nm}$; added salt would shorten it further. Thus the stated ion inventory and screening length cannot both be interpreted as a self-consistent dilute uniform screening gas with this dielectric. Bound counterions, nonuniform distributions and an effective phenomenological screening length can change that interpretation. Linear screening is itself only an estimate near highly charged DNA.