Dielectric screening is the rearrangement of a medium's charges in response to an applied charge. In a metal, mobile electrons accumulate near a positive perturbation and reduce its long-range potential. Thomas-Fermi screening treats the electron gas locally as an equilibrium gas whose density responds to the local electrostatic potential, ignoring gradients of the density on the microscopic wavelength scale.
Take , background permittivity , and electrostatic potential . An electron's potential energy is , so the local kinetic Fermi energy becomes . For a three-dimensional spin-degenerate parabolic band at zero temperature,
where is the spin-summed density of states per volume per energy. The induced charge is . Combining it with Poisson equation gives
For a point external charge , spherical symmetry and decay at infinity select
Thus the bare potential is cut off exponentially. Fourier transformation makes the same result transparent: , so the static longitudinal dielectric factor is relative to the background. For a free electron band, and .
Figure 1.
Local Thomas-Fermi screening and optical/acoustic modes of an illustrative diatomic chain
. Details.
The local linear approximation requires a weak perturbation and slowly varying fields on the relevant microscopic scale. It omits nonlocal Friedel oscillations and atom-scale structure; a nonlinear local-density treatment is needed near a sufficiently strong charge. Finite temperature replaces by the appropriate compressibility derivative . The nondegenerate limit gives classical screening rather than the degenerate Thomas-Fermi length. A larger density of states means that a smaller potential can move more carriers, hence stronger screening.
Use Thomas-Fermi screening, treating a homogeneous three-dimensional conductor in static linear response. For a positive test charge the screened electric potential and bare Coulomb potential are
Here is the background permittivity before the mobile carriers screen; it is in the simplest free-electron model. Both curves have the same leading divergence near the origin, while the screened curve lies below the bare one and decays much faster.
A positive potential attracts electrons. If is their total spin-summed density of states per volume and energy, and . Inserting this charge into Poisson equation gives the Thomas-Fermi screening length
A larger density of states lets the conductor rearrange more charge for the same potential change and shortens the screening length. For a parabolic free electron model, , hence .
Figure 1.
Bare and Thomas-Fermi-screened potentials of a positive test charge
. Original dimensionless model curves on the same axes, normalized by . The near-origin divergence is clipped for display. The local screening approximation omits microscopic nonlocal effects. Details.
The exponential form is the intended local screening model. A microscopic zero-temperature metal also has Friedel oscillations, so the exponential alone is not an exact description at every length scale.
Substitute harmonic amplitudes with the specified convention:
One electron contributes a dipole moment , so the electric polarization of a density is
Using gives
Thus the paper's dimensionless is the relative permittivity, rather than the absolute permittivity . The electron's negative charge disappears from the electric susceptibility because both the force and dipole moment contain .
For monovalent mobile ions, the linearized Debye–Hückel screening length is
Here is a number density in , not a molarity. Use the question's idealized neutral-water convention, , and . Both ion species contribute:
For added salt concentration in , the ionic strength in this model is . With the numerical molarity ,
The water contribution is negligible at that salt level. The factor converts molarity into ; omitting it or counting only one ion species gives a wrong screening length. Exact water autoionization and activities vary with temperature, so the first estimate uses the simplified neutral concentrations implied by the question. At the higher salt level, ion-size effects, nonideal activities and changes in permittivity limit the numerical accuracy of this continuum estimate.
At fixed temperature and solvent permittivity, the Debye–Hückel screening length satisfies , where is ionic strength. Fully dissociated sodium chloride has two monovalent ion species, giving M. Fully dissociated magnesium sulfate has two divalent species, giving
The multivalent-electrolyte screening-length ratio is therefore
The squared valences are essential: counting only the total ion concentration misses a factor of four in the screening strength per formula unit. This is the comparable linear-screening estimate, assuming unchanged dielectric response and full dissociation. ion pairs and concentrated-electrolyte nonideality can make a real divalent-salt solution depart from it; no such correction is specified.
Stabilization and screening. Charge stabilization of colloids means suppressing aggregation through repulsion between similarly charged particles and their electric double layers. In DLVO theory, this repulsion competes with attractive Van der Waals forces. A substantial repulsive barrier compared with thermal energy can make aggregation slow, even when the dispersed state is only metastable.
The Debye–Hückel screening length is the characteristic range over which a weak electric potential is attenuated by rearrangement of mobile ions. It measures the thickness of a diffuse screening cloud, not a molecular diameter. For a plane, a weak potential decays as into the bulk.
Derivation and the source convention. Put a planar interface at , with electrolyte in , bulk potential zero, and absolute permittivity . For a symmetric electrolyte, let be the bulk number density of each ionic species. The Boltzmann distribution gives
so the mobile charge density is
Substitution into the SI Poisson equation, , gives the Poisson-Boltzmann equation
For , linearize the hyperbolic sine. The Debye–Hückel approximation is then
Here is the ionic valence magnitude, the elementary electric charge, the Boltzmann constant, absolute temperature, the vacuum permittivity, and the solvent's relative permittivity. For several ionic species, the numerator is replaced by .
The original PDF includes an additional in the denominator. That factor does not follow from SI Poisson equation with number densities and the stated permittivities. The Bjerrum length provides a useful check: and the consistent result is . A Gaussian-unit convention moves factors of together with the charge and permittivity conventions; it does not justify inserting this extra factor while retaining all the given SI quantities. Both numerical interpretations are reported below rather than silently changing the source.
Conditions favoring stability. Increasing the surface-potential or surface-charge magnitude usually strengthens the electrostatic barrier. Lower ionic strength increases and generally makes repulsion longer ranged; added salt and especially multivalent counterions shorten it and can promote aggregation. Suitable pH can keep an ionizable surface charged, while conditions near an isoelectric state weaken charge repulsion. The solvent-dependent Van der Waals forces, charge regulation and particle geometry also affect the barrier, so a long screening length is not by itself sufficient. At fixed density and valence, the linear formula gives , but actual temperature changes can also alter the permittivity and surface chemistry.
Surface charge at small potential. Take the usual one-sided planar-interface model, with negligible electric displacement on the nonaqueous side. The solution decaying into the solution is
The electric field is . The electrostatic interface boundary conditions then give the surface charge density
One can also derive it by integrating the diffuse charge density:
Thus the diffuse layer neutralizes the surface charge. This is the planar Debye-Hückel double layer. If instead the object is an isolated sheet wetted by identical electrolyte on both sides, each side screens half its total charge and . The one-sided assumption is the conventional local model for the exterior of a large charged colloid; it must be stated because the question does not specify the material on the other side.
Numerical ion density. For , , and , the consistent SI result is
This is for each ion species, and hence also the salt concentration for a monovalent salt. The combined number density of both signs is .
Literal substitution into the formula printed in the PDF instead gives
These values are distinct consequences of the source's extra factor; they cannot both describe the same physical screening length in SI with the same definition of .
Surface potential and validity. The requested linear estimate in the one-sided model is
However, the thermal voltage is , so . Thus the small-potential approximation is not self-consistent for the supplied surface charge. The boxed value is its formal extrapolation, not a reliable prediction of the actual surface potential.
For a valid nonlinear comparison within the ideal one-sided planar Poisson-Boltzmann equation, multiply the differential equation by and use zero potential and field in the bulk. This gives
Using the surface boundary condition yields the Grahame equation and its inverse:
With the SI screening length, the latter gives . This is a nonlinear mean-field estimate under explicit assumptions, not a unique experimentally determined potential: compact interfacial layers, finite ion size and surface chemistry were not supplied. The source's normalization problem and its use of a linear formula far outside the small-potential regime are separate issues.
Let be the permittivity, , the elementary charge, and the bulk salt formula-unit number density. The Boltzmann distribution gives
Consequently the exterior Poisson-Boltzmann equation is
For small , meaning , linearize the exponentials. The Debye–Hückel approximation gives
This is also : the divalent cation contributes and the two monovalent anions contribute . In terms of molar salt concentration in , replace by ; if is in , use . The molar ionic strength is in the same concentration units.
The boundary conditions are and, by Gauss's law,
Here the surface charge is uniform and the spherically symmetric interior field is zero. Set ; the radial equation reduces to . Decay excludes the growing exponential, so write . Then , giving
It solves the linearized Poisson-Boltzmann equation, not the full nonlinear equation at finite . The radial ion profiles, consistently to first order, are
The exponential Boltzmann distribution formulas above express the full mean-field relation, but substituting a linearized potential into them is accurate only to the retained order. A useful electroneutrality check is
The result is the Debye-Hückel potential of a charged sphere.
Use SI units, uniform absolute permittivity and a neutral bulk reservoir. Local equilibrium of ideal mobile ions gives . Combining their charge with the Poisson equation yields
The reference is in the bulk. Boundary data specify fixed charge, fixed potential or a surface reaction/charge regulation law. These are physically different problems. The mean-field theory assumes point ions, thermal equilibrium and a continuum solvent; ion-size packing, dielectric structure and strong electrostatic correlations can invalidate it.
For a symmetric monovalent electrolyte let count each sign and define , . Outside fixed sources,
For , the Debye-Hückel approximation is . A planar wall bounding electrolyte in has ; assuming zero field on the other side, its charge is . A sphere of radius and total charge instead gives
The factors arise from geometry and boundary conditions; planar screening cannot be substituted unchanged for a small sphere. The weak-potential condition must hold where the potential is largest, particularly at the surface.
The nonlinear planar Poisson-Boltzmann profile is also exactly integrable. Multiply by and use at infinity:
The sign selects a profile decaying to the bulk for either sign of surface potential. Integration gives
For weak this reduces to the linear solution. Strong surface potentials have a bounded far-field amplitude , while the near-wall ion density still grows strongly. There is no need to permit the negative ion densities produced by an invalid extrapolation of a linearized density formula.
The same first integral yields the Grahame equation and its diffuse-layer differential capacitance per unit area:
These describe the diffuse layer. A compact interfacial layer adds a voltage drop and can place its capacitance in series with ; a measured electrode capacitance need not equal the ideal diffuse contribution.[ref]
For comparison, a planar counterion-only double layer next to a negative plane has positive counterions and no salt reservoir. Its profile is algebraic:
Integrating over gives , checking charge neutrality. This problem cannot be obtained by putting a finite reservoir density into the symmetric-electrolyte answer without changing its boundary conditions.
Applications include screening around synthetic colloids, diffuse-layer electrode capacitance and forces between charged plates. In planar geometry, the bulk-subtracted normal stress is constant:
At the midpoint of identical plates, symmetry sets , so the remaining excess ionic osmotic pressure determines their repulsion. A comparison with force or capacitance measurements should vary salt and surface potential independently, test charge regulation and adsorption, and distinguish compact-layer effects from failures of diffuse screening. Nonlinear Poisson-Boltzmann theory retains finite-potential Boltzmann populations, but still omits correlations and finite-ion-volume saturation.
Figure 1.
Original linear and nonlinear planar screening profiles and diffuse-layer capacitance
. Details.
The constants specify isotropic, nondispersive linear response: controls electric polarization, and controls magnetization. Their numerical relation to conventional electric susceptibility and magnetic susceptibility depends on the signs and field normalization in the action. In the rationalized units used here, the printed action gives the constitutive relations
Thus its relative permittivity and magnetic permeability are
Here conventional is defined by , not by . With the alternative repaired action that yields the PDF's target equation, is unchanged but and . In particular, a positive conventional electric susceptibility corresponds to negative in the printed conventions; the problem does not specify signs for these constants.
Permittivity Created 2026-10-05 Updated 2026-10-18
Permittivity describes how an electric field and a material's electric polarization combine into the electric displacement field. In a static, local response with linearity and isotropy, and no electric polarization at zero electric field, the constitutive equation is
Here is vacuum permittivity, while is the relative permittivity, a dimensionless quantity. The SI unit of permittivity is farad per metre, : multiplying it by an electric field in gives an electric displacement field in . A larger positive permittivity lets the same parallel-plate capacitor hold more plate electric charge at the same electric potential difference.
The material contribution is electric polarization , the electric dipole moment per unit volume. The electric displacement field and electric susceptibility satisfy
Substitution gives the first constitutive equation. The electric field here is the averaged, total field inside the material, not just the field of the external plate electric charges. The continuum approximation averages over many chemical atoms or molecules; locality means that the electric polarization at a position depends on the electric field there. Static equilibrium, fixed material conditions and sufficiently small electric fields are needed for this scalar proportionality. [ref].
For an ideal parallel-plate capacitor completely filled by a uniform dielectric, hold the free plate electric charge magnitude fixed and neglect the fringing field of a capacitor. Writing for the plate area, Gauss's law gives . Thus
The electric dipoles point along the electric field. Their surface electric charges, described by polarization density, oppose the field from the free plate electric charges. For a hypothetical dielectric with , the internal electric field is one quarter of its vacuum value, while . At fixed electric potential difference instead, the internal electric field stays approximately and the plate electric charge grows fourfold. These are different experimental constraints. [ref].
Figure 1. . An original schematic and normalized model comparison for a uniform dielectric with relative permittivity four. The arrows for electric dipole moments run from negative to positive electric charge; the two field contributions oppose each other. The bar values compare electric displacement field , , and electric polarization , which all have units . The schematic is not to microscopic scale. Details.
Permittivity is a useful material summary, not a universal constant of the material under every condition. Temperature, frequency and electric field strength can change the response. In a linear anisotropic dielectric, permittivity is a tensor, , so the electric displacement field need not be parallel to the electric field. At nonzero frequency, delayed electric polarization requires complex permittivity and may produce dielectric loss. [ref].
Piezoelectric polymer 2026-10-20
A polymer exhibiting piezoelectricity. A compliant, low-density film can accept larger bending or strain than a brittle ceramic, while low permittivity can yield a large voltage coefficient despite a modest charge coefficient. Extractable energy still depends on mechanical input, capacitance and electrical loading.
For a coaxial cable with inner and outer radii , a radial interface , and permittivities , define . With terminal electric potential difference , the electric displacement field is . The electric field is in each layer . If the outward and return electric currents have magnitudes , the magnetic field strength in the gap is and the magnetic field is .
The proof uses Gauss's law without free interface electric charge, the electric potential difference integral, and the magnetostatic Ampère law. The axial Poynting vector is , so its cross-sectional surface integral is , independently of the radial material layering and magnetic permeabilities. The empty inner and exterior regions have zero electric field and magnetic field for the regular solution with bounded exterior electric potential and no applied magnetic field.
Rayleigh scattering by a small sphere Created 2026-10-14 Updated 2026-10-20
Consider a fixed, homogeneous, nonabsorbing, nonmagnetic sphere of radius in vacuum. Let be its relative refractive index, a fixed real number; in vacuum this equals the sphere's refractive index. Treat as independent of wavelength over the range being compared, neglecting optical dispersion. Air can be approximated by vacuum here because its refractive index is close to one. These assumptions describe a toy particle, rather than a literal model of an individual air molecule.
The Rayleigh scattering cross-section is an effective area: the total scattered power is , where is the incident light intensity. For its leading approximation is
The units are area, since has units . This Rayleigh scattering cross-section measures scattering, not conversion of the incident electromagnetic radiation into heat. The definition follows [ref].
To see the mechanism, the nearly uniform applied electric field produces an electric dipole moment with peak magnitude
where is the incident electric field peak wave amplitude and is vacuum permittivity. This is the sphere's response in electrostatics, with its permittivity divided by vacuum permittivity equal to under the nonmagnetic assumption. [ref]. The averaged electric dipole radiation power is , with angular frequency and speed of light . Dividing by gives the displayed Rayleigh scattering cross-section. Thus doubling radius multiplies the electric dipole moment by eight and the scattered power by , provided the larger sphere still satisfies . [ref].
For a uniform linear medium with isotropy and local constitutive response, Maxwell equations give . Here is permittivity and is magnetic permeability. Since , the relative permittivity and the relative magnetic permeability obey
If both responses are real and positive with negligible loss, and the ordinary positive-index branch is intended, then . If additionally , then . All these responses must be evaluated at the electromagnetic wave's frequency; a static relative permittivity cannot simply be substituted at an optical frequency. [ref]
With loss, use a complex index and explicitly fix the phase convention for complex amplitudes. For with , its wave amplitude has factor . A forward attenuating branch has , while its wave phase is controlled by . The complex relation alone leaves a choice of sign: the physical solution also needs causality and the appropriate radiation or decay condition. In negative-index media the physically selected real part can be negative. The positive square root is therefore not a universal rule for complex or negative-index media. Direction and polarization of an electromagnetic wave also matter in a plane wave in an anisotropic dielectric, for which one scalar index need not describe every propagation branch.
Relative permittivity Created 2026-10-14 Updated 2026-10-20
The dimensionless ratio of a medium's permittivity to vacuum permittivity. The constant scalar value used in a Wannier-Mott exciton is an approximation to the relevant material response.
Reststrahlen effect 2026-10-25
In a polar insulating crystal, a band with negative real permittivity can occur between the transverse and longitudinal optical resonances. A lossless bulk transverse wave is evanescent there, producing high normal-incidence reflectance. Damping rounds the edges and permits absorption. The reflection minimum beyond the band can occur where the refractive index approaches one; it is not exactly the zero of permittivity.
Surface plasmon polariton Created 2026-10-18 Updated 2026-10-20
A transverse magnetic polarization electromagnetic wave bound to an interface through coupling to the material's collective charge response. Its fields propagate along the interface and decay away from it on both sides. In an ideal nonmagnetic lossless planar vacuum-material interface, a real negative relative permittivity below supports the bound branch described by the lossless planar surface-plasmon dispersion. Material losses generally make the propagation constant complex.
Thomas-Fermi screening 2026-10-18
In the long-wavelength static response of a degenerate conductor, a small electric potential changes the electron density by , where is the spin-summed density of states per volume per energy. The induced charge is therefore . With background permittivity , Poisson equation becomes
A point charge then has . Larger electronic density of states gives stronger screening. This local response approximation does not capture the microscopic oscillatory tail of Friedel oscillations.[ref]
Transverse electrostatic field inside a dielectric cylinder Created 2026-10-20 Updated 2026-10-20
For an infinitely long circular dielectric cylinder in a uniform transverse field in vacuum, the inside electrostatic potential is linear in the transverse coordinate. Matching the tangential electric field and normal electric displacement field to the outside dipolar solution gives the displayed uniform internal field. The result assumes no free interface charge and neglects end effects; an exterior dielectric changes the permittivity ratio.
Replace two well-separated cylindrical conductors by opposite axial line charges. Their surface-to-surface potential difference is , giving the displayed capacitance per unit length. The approximation needs and a homogeneous surrounding permittivity.