Following are the definitions for continuous charge distributions.67
The linear charge density is the ratio of an infinitesimal electric charge dQ (SI unit: C) to an infinitesimal line element, λ q = d Q d ℓ , {\displaystyle \lambda _{q}={\frac {dQ}{d\ell }}\,,} similarly the surface charge density uses a surface area element dS σ q = d Q d S , {\displaystyle \sigma _{q}={\frac {dQ}{dS}}\,,} and the volume charge density uses a volume element dV ρ q = d Q d V , {\displaystyle \rho _{q}={\frac {dQ}{dV}}\,,}
Integrating the definitions gives the total charge Q of a region according to line integral of the linear charge density λq(r) over a line or 1d curve C, Q = ∫ L λ q ( r ) d ℓ {\displaystyle Q=\int _{L}\lambda _{q}(\mathbf {r} )\,d\ell } similarly a surface integral of the surface charge density σq(r) over a surface S, Q = ∫ S σ q ( r ) d S {\displaystyle Q=\int _{S}\sigma _{q}(\mathbf {r} )\,dS} and a volume integral of the volume charge density ρq(r) over a volume V, Q = ∫ V ρ q ( r ) d V {\displaystyle Q=\int _{V}\rho _{q}(\mathbf {r} )\,dV} where the subscript q is to clarify that the density is for electric charge, not other densities like mass density, number density, probability density, and prevent conflict with the many other uses of λ, σ, ρ in electromagnetism for wavelength, electrical resistivity and conductivity.
Within the context of electromagnetism, the subscripts are usually dropped for simplicity: λ, σ, ρ. Other notations may include: ρℓ, ρs, ρv, ρL, ρS, ρV etc.
The total charge divided by the length, surface area, or volume will be the average charge densities: ⟨ λ q ⟩ = Q ℓ , ⟨ σ q ⟩ = Q S , ⟨ ρ q ⟩ = Q V . {\displaystyle \langle \lambda _{q}\rangle ={\frac {Q}{\ell }}\,,\quad \langle \sigma _{q}\rangle ={\frac {Q}{S}}\,,\quad \langle \rho _{q}\rangle ={\frac {Q}{V}}\,.}
In dielectric materials, the total charge of an object can be separated into "free" and "bound" charges.
Bound charges set up electric dipoles in response to an applied electric field E, and polarize other nearby dipoles tending to line them up, the net accumulation of charge from the orientation of the dipoles is the bound charge. They are called bound because they cannot be removed: in the dielectric material the charges are the electrons bound to the nuclei.8
Free charges are the excess charges which can move into electrostatic equilibrium, i.e. when the charges are not moving and the resultant electric field is independent of time, or constitute electric currents.9
In terms of volume charge densities, the total charge density is: ρ = ρ f + ρ b . {\displaystyle \rho =\rho _{\text{f}}+\rho _{\text{b}}\,.} as for surface charge densities: σ = σ f + σ b . {\displaystyle \sigma =\sigma _{\text{f}}+\sigma _{\text{b}}\,.} where subscripts "f" and "b" denote "free" and "bound" respectively.
The bound surface charge is the charge piled up at the surface of the dielectric, given by the dipole moment perpendicular to the surface:10 q b = d ⋅ n ^ | s | {\displaystyle q_{b}={\frac {\mathbf {d} \cdot \mathbf {\hat {n}} }{|\mathbf {s} |}}} where s is the separation between the point charges constituting the dipole, d {\displaystyle \mathbf {d} } is the electric dipole moment, n ^ {\displaystyle \mathbf {\hat {n}} } is the unit normal vector to the surface.
Taking infinitesimals: d q b = d d | s | ⋅ n ^ {\displaystyle dq_{b}={\frac {d\mathbf {d} }{|\mathbf {s} |}}\cdot \mathbf {\hat {n}} } and dividing by the differential surface element dS gives the bound surface charge density: σ b = d q b d S = d d | s | d S ⋅ n ^ = d d d V ⋅ n ^ = P ⋅ n ^ . {\displaystyle \sigma _{b}={\frac {dq_{b}}{dS}}={\frac {d\mathbf {d} }{|\mathbf {s} |dS}}\cdot \mathbf {\hat {n}} ={\frac {d\mathbf {d} }{dV}}\cdot \mathbf {\hat {n}} =\mathbf {P} \cdot \mathbf {\hat {n}} \,.} where P is the polarization density, i.e. density of electric dipole moments within the material, and dV is the differential volume element.
Using the divergence theorem, the bound volume charge density within the material is q b = ∫ ρ b d V = − ∮ S P ⋅ n ^ d S = − ∫ ∇ ⋅ P d V {\displaystyle q_{b}=\int \rho _{b}\,dV=-\oint _{S}\mathbf {P} \cdot {\hat {\mathbf {n} }}\,dS=-\int \nabla \cdot \mathbf {P} \,dV} hence: ρ b = − ∇ ⋅ P . {\displaystyle \rho _{b}=-\nabla \cdot \mathbf {P} \,.}
The negative sign arises due to the opposite signs on the charges in the dipoles, one end is within the volume of the object, the other at the surface.
A more rigorous derivation is given below.11
The electric potential due to a dipole moment d is: φ = 1 4 π ε 0 ( r − r ′ ) ⋅ d | r − r ′ | 3 {\displaystyle \varphi ={\frac {1}{4\pi \varepsilon _{0}}}{\frac {(\mathbf {r} -\mathbf {r} ')\cdot \mathbf {d} }{|\mathbf {r} -\mathbf {r} '|^{3}}}}
For a continuous distribution, the material can be divided up into infinitely many infinitesimal dipoles d d = P d V = P d 3 r {\displaystyle d\mathbf {d} =\mathbf {P} dV=\mathbf {P} d^{3}\mathbf {r} } where dV = d3r′ is the volume element, so the potential is the volume integral over the object: φ = 1 4 π ε 0 ∭ ( r − r ′ ) ⋅ P | r − r ′ | 3 d 3 r ′ {\displaystyle \varphi ={\frac {1}{4\pi \varepsilon _{0}}}\iiint {\frac {(\mathbf {r} -\mathbf {r} ')\cdot \mathbf {P} }{|\mathbf {r} -\mathbf {r} '|^{3}}}d^{3}\mathbf {r'} }
Since ∇ ′ ( 1 | r − r ′ | ) ≡ ( e x ∂ ∂ x ′ + e y ∂ ∂ y ′ + e z ∂ ∂ z ′ ) ( 1 | r − r ′ | ) = r − r ′ | r − r ′ | 3 {\displaystyle \nabla '\left({\frac {1}{|\mathbf {r} -\mathbf {r} '|}}\right)\equiv \left(\mathbf {e} _{x}{\frac {\partial }{\partial x'}}+\mathbf {e} _{y}{\frac {\partial }{\partial y'}}+\mathbf {e} _{z}{\frac {\partial }{\partial z'}}\right)\left({\frac {1}{|\mathbf {r} -\mathbf {r} '|}}\right)={\frac {\mathbf {r} -\mathbf {r} '}{|\mathbf {r} -\mathbf {r} '|^{3}}}} where ∇′ is the gradient in the r′ coordinates, φ = 1 4 π ε 0 ∭ P ⋅ ∇ ′ ( 1 | r − r ′ | ) d 3 r ′ {\displaystyle \varphi ={\frac {1}{4\pi \varepsilon _{0}}}\iiint \mathbf {P} \cdot \nabla '\left({\frac {1}{|\mathbf {r} -\mathbf {r} '|}}\right)d^{3}\mathbf {r'} }
Integrating by parts φ = 1 4 π ε 0 ∭ [ ∇ ′ ⋅ ( P | r − r ′ | ) − 1 r − r ′ ( ∇ ′ ⋅ P ) ] d 3 r ′ {\displaystyle \varphi ={\frac {1}{4\pi \varepsilon _{0}}}\iiint \left[\nabla '\cdot \left({\frac {\mathbf {P} }{|\mathbf {r} -\mathbf {r} '|}}\right)-{\frac {1}{\mathbf {r} -\mathbf {r} '}}(\nabla '\cdot \mathbf {P} )\right]d^{3}\mathbf {r'} } using the divergence theorem:
which separates into the potential of the surface charge (surface integral) and the potential due to the volume charge (volume integral):
that is σ b = P ⋅ n ^ , ρ b = − ∇ ⋅ P {\displaystyle \sigma _{b}=\mathbf {P} \cdot \mathbf {\hat {n}} \,,\quad \rho _{b}=-\nabla \cdot \mathbf {P} }
The free charge density serves as a useful simplification in Gauss's law for electricity; the volume integral of it is the free charge enclosed in a charged object - equal to the net flux of the electric displacement field D emerging from the object:
See Maxwell's equations and constitutive relation for more details.
For the special case of a homogeneous charge density ρ0, independent of position i.e. constant throughout the region of the material, the equation simplifies to: Q = V ρ 0 . {\displaystyle Q=V\rho _{0}.}
Start with the definition of a continuous volume charge density: Q = ∫ V ρ q ( r ) d V . {\displaystyle Q=\int _{V}\rho _{q}(\mathbf {r} )\,dV.}
Then, by definition of homogeneity, ρq(r) is a constant denoted by ρq, 0 (to differ between the constant and non-constant densities), and so by the properties of an integral can be pulled outside of the integral resulting in: Q = ρ q , 0 ∫ V d V = ρ 0 V {\displaystyle Q=\rho _{q,0}\int _{V}\,dV=\rho _{0}V} so, Q = V ρ q , 0 . {\displaystyle Q=V\rho _{q,0}.}
The equivalent proofs for linear charge density and surface charge density follow the same arguments as above.
For a single point charge q at position r0 inside a region of 3d space R, like an electron, the volume charge density can be expressed by the Dirac delta function: ρ q ( r ) = q δ ( r − r 0 ) {\displaystyle \rho _{q}(\mathbf {r} )=q\delta (\mathbf {r} -\mathbf {r} _{0})} where r is the position to calculate the charge.
As always, the integral of the charge density over a region of space is the charge contained in that region. The delta function has the shifting property for any function f: ∫ R d 3 r f ( r ) δ ( r − r 0 ) = f ( r 0 ) {\displaystyle \int _{R}d^{3}\mathbf {r} f(\mathbf {r} )\delta (\mathbf {r} -\mathbf {r} _{0})=f(\mathbf {r} _{0})} so the delta function ensures that when the charge density is integrated over R, the total charge in R is q: Q = ∫ R d 3 r ρ q = ∫ R d 3 r q δ ( r − r 0 ) = q ∫ R d 3 r δ ( r − r 0 ) = q {\displaystyle Q=\int _{R}d^{3}\mathbf {r} \,\rho _{q}=\int _{R}d^{3}\mathbf {r} \,q\delta (\mathbf {r} -\mathbf {r} _{0})=q\int _{R}d^{3}\mathbf {r} \,\delta (\mathbf {r} -\mathbf {r} _{0})=q}
This can be extended to N discrete point-like charge carriers. The charge density of the system at a point r is a sum of the charge densities for each charge qi at position ri, where i = 1, 2, ..., N: ρ q ( r ) = ∑ i = 1 N q i δ ( r − r i ) {\displaystyle \rho _{q}(\mathbf {r} )=\sum _{i=1}^{N}\ q_{i}\delta (\mathbf {r} -\mathbf {r} _{i})}
The delta function for each charge qi in the sum, δ(r − ri), ensures the integral of charge density over R returns the total charge in R: Q = ∫ R d 3 r ∑ i = 1 N q i δ ( r − r i ) = ∑ i = 1 N q i ∫ R d 3 r δ ( r − r i ) = ∑ i = 1 N q i {\displaystyle Q=\int _{R}d^{3}\mathbf {r} \sum _{i=1}^{N}\ q_{i}\delta (\mathbf {r} -\mathbf {r} _{i})=\sum _{i=1}^{N}\ q_{i}\int _{R}d^{3}\mathbf {r} \delta (\mathbf {r} -\mathbf {r} _{i})=\sum _{i=1}^{N}\ q_{i}}
If all charge carriers have the same charge q (for electrons q = −e, the electron charge) the charge density can be expressed through the number of charge carriers per unit volume, n(r), by ρ q ( r ) = q n ( r ) . {\displaystyle \rho _{q}(\mathbf {r} )=qn(\mathbf {r} )\,.}
Similar equations are used for the linear and surface charge densities.
Further information: classical electromagnetism and special relativity and relativistic electromagnetism
In special relativity, the length of a segment of wire depends on velocity of observer because of length contraction, so charge density will also depend on velocity. Anthony French12 has described how the magnetic field force of a current-bearing wire arises from this relative charge density. He used (p 260) a Minkowski diagram to show "how a neutral current-bearing wire appears to carry a net charge density as observed in a moving frame." When a charge density is measured in a moving frame of reference it is called proper charge density.131415
It turns out the charge density ρ and current density J transform together as a four-current vector under Lorentz transformations.
See also: Density functional theory and Hartree–Fock method
In quantum mechanics, charge density ρq is related to wavefunction ψ(r) by the equation ρ q ( r ) = q | ψ ( r ) | 2 {\displaystyle \rho _{q}(\mathbf {r} )=q|\psi (\mathbf {r} )|^{2}} where q is the charge of the particle and |ψ(r)|2 = ψ*(r)ψ(r) is the probability density function i.e. probability per unit volume of a particle located at r. When the wavefunction is normalized - the average charge in the region r ∈ R is Q = ∫ R q | ψ ( r ) | 2 d 3 r {\displaystyle Q=\int _{R}q|\psi (\mathbf {r} )|^{2}\,d^{3}\mathbf {r} } where d3r is the integration measure over 3d position space.
For system of identical fermions, the number density is given as sum of probability density of each particle in :
n ( r ) = ∑ i ⟨ ψ | δ 3 ( r − r i ′ ) | ψ ⟩ {\displaystyle n(\mathbf {r} )=\sum _{i}\langle \psi |\delta ^{3}(\mathbf {r} -\mathbf {r} _{i}')|\psi \rangle }
n ( r ) = ∑ i ∫ d 3 r 2 ⋯ ∫ d 3 r N Ψ ∗ ( r , r 2 , … , r i = r ′ , … , r N ) Ψ ( r , r 2 , … , r i = r ′ , … , r N ) . {\displaystyle n(\mathbf {r} )=\sum _{i}\int {\mathrm {d} }^{3}\mathbf {r} _{2}\cdots \int {\mathrm {d} }^{3}\mathbf {r} _{N}\,\Psi ^{*}(\mathbf {r} ,\mathbf {r} _{2},\dots ,\mathbf {r} _{i}=\mathbf {r} ',\dots ,\mathbf {r} _{N})\Psi (\mathbf {r} ,\mathbf {r} _{2},\dots ,\mathbf {r} _{i}=\mathbf {r} ',\dots ,\mathbf {r} _{N}).}
Using symmetrization condition: n ( r ) = N ∫ d 3 r 2 ⋯ ∫ d 3 r N Ψ ∗ ( r , r 2 , … , r N ) Ψ ( r , r 2 , … , r N ) . {\displaystyle n(\mathbf {r} )=N\int {\mathrm {d} }^{3}\mathbf {r} _{2}\cdots \int {\mathrm {d} }^{3}\mathbf {r} _{N}\,\Psi ^{*}(\mathbf {r} ,\mathbf {r} _{2},\dots ,\mathbf {r} _{N})\Psi (\mathbf {r} ,\mathbf {r} _{2},\dots ,\mathbf {r} _{N}).} where ρ q ( r ) = q ⋅ n ( r ) {\textstyle \rho _{q}(\mathbf {r} )=q\cdot n(\mathbf {r} )} is considered as the charge density.
The potential energy of a system is written as: ⟨ ψ | U | ψ ⟩ = ∫ V ( r ) n ( r ) δ 3 r {\displaystyle \langle \psi |U|\psi \rangle =\int V(\mathbf {r} )n(\mathbf {r} )\delta ^{3}\mathbf {r} } The electron-electron repulsion energy is thus derived under these conditions to be: U e e [ n ] = J [ n ] = 1 2 ∫ δ 3 r ′ ∫ δ 3 r ( ( e n ( r ) ) ( e n ( r ′ ) ) | r − r ′ | ) = 1 2 ∫ δ 3 r ′ ∫ δ 3 r ( ρ ( r ) ρ ( r ′ ) | r − r ′ | ) {\displaystyle U_{ee}[n]=J[n]={\frac {1}{2}}\int \delta ^{3}\mathbf {r} '\int \delta ^{3}\mathbf {r} \left({\frac {(en(\mathbf {r} ))(en(\mathbf {r} '))}{|\mathbf {r} -\mathbf {r} '|}}\right)={\frac {1}{2}}\int \delta ^{3}\mathbf {r} '\int \delta ^{3}\mathbf {r} \left({\frac {\rho (\mathbf {r} )\rho (\mathbf {r} ')}{|\mathbf {r} -\mathbf {r} '|}}\right)} Note that this is excluding the exchange energy of the system, which is a purely quantum mechanical phenomenon, has to be calculated separately.
Then, the energy is given using Hartree-Fock method as:
E [ n ] = I + J − K {\displaystyle E[n]=I+J-K}
Where I is the kinetic and potential energy of electrons due to positive charges, J is the electron electron interaction energy and K is the exchange energy of electrons.1617
The charge density appears in the continuity equation for electric current, and also in Maxwell's Equations. It is the principal source term of the electromagnetic field; when the charge distribution moves, this corresponds to a current density. The charge density of molecules impacts chemical and separation processes. For example, charge density influences metal-metal bonding and hydrogen bonding.18 For separation processes such as nanofiltration, the charge density of ions influences their rejection by the membrane.19
P.M. Whelan, M.J. Hodgson (1978). Essential Principles of Physics (2nd ed.). John Murray. ISBN 0-7195-3382-1. 0-7195-3382-1 ↩
"Physics 2: Electricity and Magnetism, Course Notes, Ch. 2, p. 15-16" (PDF). MIT OpenCourseware. Massachusetts Institute of Technology. 2007. Retrieved December 3, 2017. https://ocw.mit.edu/courses/physics/8-02-physics-ii-electricity-and-magnetism-spring-2007/readings/summary_w01d2.pdf ↩
Serway, Raymond A.; Jewett, John W. (2013). Physics for Scientists and Engineers, Vol. 2, 9th Ed. Cengage Learning. p. 704. ISBN 9781133954149. 9781133954149 ↩
Purcell, Edward (2011-09-22). Electricity and Magnetism. Cambridge University Press. ISBN 9781107013605. 9781107013605 ↩
I.S. Grant; W.R. Phillips (2008). Electromagnetism (2nd ed.). Manchester Physics, John Wiley & Sons. ISBN 978-0-471-92712-9. 978-0-471-92712-9 ↩
D.J. Griffiths (2007). Introduction to Electrodynamics (3rd ed.). Pearson Education, Dorling Kindersley. ISBN 978-81-7758-293-2. 978-81-7758-293-2 ↩
French, A. (1968). "8:Relativity and electricity". Special Relativity. W. W. Norton. pp. 229–265. /wiki/W._W._Norton ↩
Mould, Richard A. (2001). "Lorentz force". Basic Relativity. Springer Science & Business Media. ISBN 0-387-95210-1. 0-387-95210-1 ↩
Lawden, Derek F. (2012). An Introduction to Tensor Calculus: Relativity and Cosmology. Courier Corporation. p. 74. ISBN 978-0-486-13214-3. 978-0-486-13214-3 ↩
Vanderlinde, Jack (2006). "11.1:The Four-potential and Coulomb's Law". Classical Electromagnetic Theory. Springer Science & Business Media. p. 314. ISBN 1-4020-2700-1. 1-4020-2700-1 ↩
Sakurai, Jun John; Napolitano, Jim (2021). Modern quantum mechanics (3rd ed.). Cambridge: Cambridge University Press. pp. 443–453. ISBN 978-1-108-47322-4. 978-1-108-47322-4 ↩
Littlejohn, Robert G. "The Hartree-Fock Method in Atoms" (PDF). https://bohr.physics.berkeley.edu/classes/221/1112/notes/hartfock.pdf ↩
R. J. Gillespie & P. L. A. Popelier (2001). "Chemical Bonding and Molecular Geometry". Environmental Science & Technology. 52 (7). Oxford University Press: 4108–4116. Bibcode:2018EnST...52.4108E. doi:10.1021/acs.est.7b06400. PMID 29510032. /wiki/Bibcode_(identifier) ↩
Razi Epsztein, Evyatar Shaulsky, Nadir Dizge, David M Warsinger, Menachem Elimelech (2018). "Ionic Charge Density-Dependent Donnan Exclusion in Nanofiltration of Monovalent Anions". Environmental Science & Technology. 52 (7): 4108–4116. Bibcode:2018EnST...52.4108E. doi:10.1021/acs.est.7b06400. PMID 29510032.{{cite journal}}: CS1 maint: multiple names: authors list (link) /wiki/Bibcode_(identifier) ↩