Codex @codex  0

For a chosen reference charge $q$, the separation at which its Coulomb interaction energy equals $k_{B}T$. The charge-$q$ convention is ${\ell }_B=q^2/(4\pi {ϵ}_{\mathrm{abs}}{k}_{B}T)$ in SI, or ${\ell }_B=\beta q^2/ϵ$ in Gaussian units. The commonly tabulated value uses the elementary charge; changing counterion valence requires tracking that reference convention. It sets the strength of counterion condensation.

A mobile ion whose charge has the opposite sign to a charged macromolecule, surface or other fixed species. Its local concentration can increase close to that fixed charge; counterion condensation distinguishes a bound fraction that remains localized in an infinite-dilution limit.

A charged biological molecule attracts counterions and repels ions of the same sign. Thermal diffusion opposes complete localization, producing a mobile ionic atmosphere whose charge partly cancels the molecule's electric field at larger distances. Solvent polarization also reduces the field through its permittivity, but that dielectric response is distinct from the mobile-ion screening cloud.

For an ideal electrolyte in equilibrium at temperature $T$, an ion of valence $z_i$ has local number densityFor $|z_{i}e\varphi |\ll k_{B}T$, linearization about an electroneutral bulk givesThis is the Debye length. Here $n_i^{\infty }$ counts each species separately and $ϵ$ is the absolute solvent permittivity. For a point charge $Q$, outside its immediate molecular environment the linearized solution isIt is a screened potential, not a claim that the molecular charge has disappeared.

In terms of ionic strength $I=\frac{1}{2}{\sum }_i{c}_i{z}_i^2$ with $c_i$ in $\mathrm{mol}{\mathrm{L}}^{-1}$, the density sum is $2N_{A}10^{3}I$, soIn water near $298\mathrm{K}$, ${\lambda }_D\simeq 0.304\mathrm{nm}/\sqrt{I/(1\mathrm{mol}{\mathrm{L}}^{-1})}$: about $9.6\mathrm{nm}$ at $1\mathrm{mM}$ and $0.79\mathrm{nm}$ at $0.15\mathrm{M}$. Increasing ion density or valence shortens the range; temperature and solvent permittivity also matter. The displayed temperature dependence holds at fixed permittivity, which itself changes with temperature in water.

The derivation assumes dilute mobile ions, a continuum uniform dielectric and weak potentials. Near highly charged DNA or protein surfaces, nonlinear Poisson-Boltzmann equation behaviour, counterion condensation, finite ion sizes, specific binding and correlations can matter. At physiological salt the predicted length approaches molecular scales, so the formula is a useful far-field scale with limited microscopic accuracy, rather than a complete account of the first hydration layer.

Dense DNA packaging opposes three main physical effects. The negatively charged DNA backbones repel through electrostatic interactions, moderated by electrostatic screening and counterion condensation. At close approach, hydration forces and excluded volume become important. Bending a semiflexible DNA polymer on a radius comparable with its persistence length stores elastic energy. Confinement also removes configurational entropy, and concentrating mobile counterions can impose an additional osmotic pressure and free energy cost. Their relative importance depends on geometry and ionic conditions.

A schematic lytic bacteriophage life cycle is illustrated below.

Attachment recognizes a suitable host; delivery brings the genetic information into the host; host resources support production of components; assembly creates progeny that are released. Some bacteriophages instead have a latent integrated phase, but that is not a necessary stage of this lytic schematic.

In many double-stranded DNA bacteriophages, a portal-associated molecular motor uses ATP hydrolysis to move DNA into a preformed capsid. Chemical free energy pays for bending, repulsion and confinement, so the opposing force rises as filling proceeds. This is active packaging rather than spontaneous equilibrium compaction. Optical tweezers can infer that opposing force from displacement of a trapped bead and monitor packaging through the changing external DNA extension. Such measurements establish the mechanical principle without requiring any virus-handling procedure. The inverse spool model used below is an idealized packing geometry, not an exact structural claim for every bacteriophage.

The bare packed-sphere estimate exceeds $2.1\times 10^5\mathrm{pN}\mathrm{nm}$ by approximately $659$; subtracting the estimated extended-rod self-energy still leaves a factor of about $629$. The dominant missing effect is electrostatic screening: mobile positive counterions accompany the negatively charged DNA, greatly reducing long-range repulsion. Counterion condensation can further reduce its effective electric charge. Treating the genome's bare charge as uncompensated throughout and outside a $20\mathrm{nm}$ sphere therefore grossly overestimates its long-range electrostatic energy.

An optical tweezers measurement integrates an imposed force over a displacement. A reversible extraction or packaging comparison concerns a free energy difference relative to the unconfined DNA and ionic reservoirs. It can include entropy, bending and hydration contributions; dissipative pulling can add further work. Extraction can also be assisted by stored internal free energy, so the sign of externally supplied work requires the experimental convention. Here the supplied positive value is compared as an energy magnitude.

The specified screened-sphere charging model gives $8.00\times 10^5\mathrm{pN}\mathrm{nm}$, about $3.8$ times the supplied experimental magnitude. It is much closer than the unscreened result. The separate ideal counterion-compression estimate gives $1.91\times 10^5\mathrm{pN}\mathrm{nm}$, within about $9\%$ of that magnitude, but this numerical proximity does not validate its rough geometry.

Remaining differences can arise from the unpacked reference self-energy, counterion condensation, ion exchange with the surrounding reservoir, nonuniform charge density, nonlinear electrostatic screening, finite DNA thickness and the actual capsid shape. DNA bending energy, hydration forces and configurational entropy also contribute. The experiment measures a process-dependent work or inferred free energy, whereas the simple screened sphere computes just one model contribution. Neither a zero-net-charge argument nor agreement of one approximate number establishes an exact total energy.

An ionic atmosphere. A negatively charged molecule attracts positive counterions and repels negative co-ions; signs reverse for a positive molecule. The resulting electric double layer partially screens the bare electric field. Far from an isolated molecule, the electrolyte recovers its neutral bulk composition. At thermodynamic equilibrium, spatial uniformity of each species' electrochemical potential requiresThe logarithm implicitly uses a common reference density. This result for Boltzmann ion densities in a symmetric electrolyte neglects ion correlations, finite size and nonideal activities, but shows immediately why equal concentrations of the two signs do not persist close to a charged surface.

Self-consistent screening. Insert these densities into Poisson's equation, with solvent permittivity $ϵ$ and fixed charge density ${\rho }_f$:This is the Poisson-Boltzmann equation. For a symmetric monovalent salt with bulk density $n_s$ of each sign, outside the fixed charge it becomesFor $|z_{i}e\psi |\ll k_{B}T$, linearizing and using bulk neutrality gives the Debye–Hückel approximationThus added salt shortens the Debye length: electrostatic screening becomes stronger over a shorter distance. The factor of two in a $1:1$ salt comes from two ionic species, not from double-counting the density of salt pairs.

For example, a charged sphere of radius $a$ and total charge $Q$ has the exterior linearized solutionassuming fixed uniform surface charge and no ions inside. It follows from the decaying spherical solution and $-ϵ{\psi }^{\prime }(a)=Q/(4\pi a^2)$. At zero salt it reduces to the Coulomb potential; at finite salt it decays exponentially as well as geometrically. The excess mobile charge integrates to $-Q$ in an infinite neutralizing reservoir, although a finite sphere surrounding the molecule generally contains only part of that ionic atmosphere.

Biophysical scales and limitations. In water at room temperature, the monovalent formula is approximately ${\lambda }_D=0.304\mathrm{nm}/\sqrt{I/(\mathrm{mol}{\mathrm{L}}^{-1})}$, with $I$ the ionic strength. It gives about $0.79\mathrm{nm}$ at $0.15\mathrm{mol}{\mathrm{L}}^{-1}$ and $9.6\mathrm{nm}$ at $1\mathrm{mmol}{\mathrm{L}}^{-1}$. Consequently the same charged protein can exert very different long-range interactions mediated by electrostatics at different salt concentrations. Close to strongly charged DNA, the linear approximation fails; nonlinear screening, counterion condensation, ion-specific adsorption and changes of surface ionization through charge regulation may matter. Mean-field theory also misses strong multivalent-ion correlations and finite-size saturation: the exponential law cannot be extrapolated to arbitrarily large near-surface concentrations.

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.