Solution (source code)

= Solution

In a conductor, mobile <electrons> redistribute in response to an <electric potential>, producing an induced <charge density> that opposes its source. Finite electronic compressibility means a finite density change requires a change of local kinetic <Fermi energy>. Screening is therefore not perfect at every microscopic distance. On scales comparable to atomic structure or the <Fermi wavelength>, a continuum local-density model itself fails.

<Thomas-Fermi screening> treats electrons as a locally uniform degenerate <free electron gas>, replaces the ionic background by a smooth neutralizing charge, and neglects nonlocal response, exchange and correlation corrections. The electrochemical <chemical potential> is uniform in equilibrium. Because an electron has charge $-e$, its electrostatic energy is $-eV$, giving
$$
\mu=E_F(n)-eV,\qquad E_F(n)=\frac{\hbar^2}{2m}(3\pi^2n)^{2/3}.
$$
At zero potential $\mu=E_{F0}$. Linearizing at fixed $\mu$, with $|eV|\ll E_{F0}$, gives
$$
\frac{2E_{F0}}{3n_0}\delta n=eV,\qquad
\boxed{\rho_{\rm ind}=-e\delta n=-\frac{3e^2n_0}{2E_{F0}}V.}
$$
Insert this in <Poisson's equation> to obtain $(\nabla^2-\xi^{-2})V=-\rho_{\rm free}/\epsilon_0$, where the <Thomas-Fermi screening length> is
$$
\boxed{\xi=\sqrt{\frac{2\epsilon_0E_{F0}}{3e^2n_0}}.}
$$
Away from a point source, direct radial differentiation shows $\nabla^2(e^{-r/\xi}/r)=\xi^{-2}e^{-r/\xi}/r$. The decaying solution is therefore $V=Ce^{-r/\xi}/r$. Integrating <Poisson's equation> over a small sphere containing charge $Q$ fixes $C=Q/(4\pi\epsilon_0)$. For a sufficiently large source, the linearized formula is only a far-field approximation, not a justified description of its singular core.

For monovalent <sodium>, the electron number density is its atomic number density,
$$
n_0=\frac{968}{0.023}N_A=2.53\times10^{28}\ \mathrm{m^{-3}}.
$$
This gives $E_{F0}=3.15\ \mathrm{eV}$ and
$$
\boxed{\xi=6.76\times10^{-11}\ \mathrm m=0.0676\ \mathrm{nm}.}
$$
Room temperature is far below the corresponding <Fermi temperature>, so degeneracy is well justified, and sodium's simple conduction band makes a free-electron estimate useful. However $k_F\xi=0.615$ and the mean atomic spacing is $n_0^{-1/3}=0.340\ \mathrm{nm}$. The potential varies on a scale shorter than $k_F^{-1}$ and the atomic spacing, so the slowly varying local-density and continuum assumptions are not quantitatively controlled at the screening length itself. A nonlocal treatment includes <Friedel oscillations>, absent from the pure exponential.

Nearly unchanged interlayer force constants in two-layer and twenty-layer films are consistent with very local screening and bonding forces: reducing the distant bulk need not strongly change the force between adjacent layers. They also indicate that metallic electronic response can survive in the thin film. The observation does not by itself validate the numerical Thomas-Fermi approximation, establish an identical full phonon spectrum, or rule out surface-induced electronic changes; it directly constrains the measured local force constant.