Gauss's law

q

{\displaystyle q}

) should have units of force*
*Electric field (

E

{\displaystyle E}

) should have units of force per area (the same as stress)*
*Electric constant (

ϵ
0

{\displaystyle \epsilon\_0}

) should have units of force per area (the same as bulk modulus)*
*Magnetic field (

B

{\displaystyle B}

) should have units of force per velocity per time^2*
*Magnetic constant (

μ
0

{\displaystyle \mu\_0}

) should have units of 1/(force per time^2)*
See this for more details

In physics, **Gauss's law**, also known as **Gauss's flux theorem**, is a law relating the distribution of electric charge to the resulting electric field. In its integral form, it states that the flux of the electric field out of an arbitrary closed surface is proportional to the electric charge enclosed by the surface, irrespective of how that charge is distributed. Even though the law alone is insufficient to determine the electric field across a surface enclosing any charge distribution, this may be possible in cases where symmetry mandates uniformity of the field. Where no such symmetry exists, Gauss's law can be used in its differential form, which states that the divergence of the electric field is proportional to the local density of charge.

The law was first formulated by Joseph-Louis Lagrange in 1773, followed by Carl Friedrich Gauss in 1813, both in the context of the attraction of ellipsoids. It is one of Maxwell's four equations, which form the basis of classical electrodynamics. Gauss's law can be used to derive Coulomb's law, and vice versa.

## Qualitative description

In words, Gauss's law states that

*The net electric flux through any hypothetical closed surface is equal to*

1

ε
0

{\displaystyle \frac{1}{\varepsilon\_0}}

*times the net electric charge within that closed surface*.

Gauss's law has a close mathematical similarity with a number of laws in other areas of physics, such as Gauss's law for magnetism and Gauss's law for gravity. In fact, any inverse-square law can be formulated in a way similar to Gauss's law: for example, Gauss's law itself is essentially equivalent to the inverse-square Coulomb's law, and Gauss's law for gravity is essentially equivalent to the inverse-square Newton's law of gravity.

The law can be expressed mathematically using vector calculus in integral form and differential form; both are equivalent since they are related by the divergence theorem, also called Gauss's theorem. Each of these forms in turn can also be expressed two ways: In terms of a relation between the electric field **E** and the total electric charge, or in terms of the electric displacement field **D** and the *free* electric charge.

## Equation involving the E field

Gauss's law can be stated using either the electric field **E** or the electric displacement field **D**. This section shows some of the forms with **E**; the form with **D** is below, as are other forms with **E**.

### Integral form

Gauss's law may be expressed as:

Φ
E

=

Q

ε
0

{\displaystyle \Phi\_E = \frac{Q}{\varepsilon\_0}}

where Φ*E* is the electric flux through a closed surface S enclosing any volume V, Q is the total charge enclosed within V, and *ε*0 is the electric constant. The electric flux Φ*E* is defined as a surface integral of the electric field:

Φ
E

=

{\displaystyle \Phi\_E=}

S

{\displaystyle \scriptstyle \_S}

E

⋅

d

A

{\displaystyle \mathbf{E} \cdot \mathrm{d}\mathbf{A}}

where **E** is the electric field, d**A** is a vector representing an infinitesimal element of area of the surface, and · represents the dot product of two vectors.

Since the flux is defined as an *integral* of the electric field, this expression of Gauss's law is called the *integral form*.

In problems involving conductors set at known potentials, the potential away from them is obtained by solving Laplace's equation, either analytically or numerically. The electric field is then calculated as the potential's negative gradient. Gauss's law makes it possible to find the distribution of electric charge: The charge in any given region of the conductor can be deduced by integrating the electric field to find the flux through a small box whose sides are perpendicular to the conductor's surface and by noting that the electric field is perpendicular to the surface, and zero inside the conductor.

The reverse problem, when the electric charge distribution is known and the electric field must be computed, is much more difficult. The total flux through a given surface gives little information about the electric field, and can go in and out of the surface in arbitrarily complicated patterns.

An exception is if there is some symmetry in the problem, which mandates that the electric field passes through the surface in a uniform way. Then, if the total flux is known, the field itself can be deduced at every point. Common examples of symmetries which lend themselves to Gauss's law include: cylindrical symmetry, planar symmetry, and spherical symmetry. See the article Gaussian surface for examples where these symmetries are exploited to compute electric fields.

### Differential form

By the divergence theorem, Gauss's law can alternatively be written in the *differential form*:

∇
⋅

E

=

ρ

ε
0

{\displaystyle \nabla \cdot \mathbf{E} = \frac{\rho}{\varepsilon\_0}}

where ∇ · **E** is the divergence of the electric field, *ε*0 is the electric constant, and ρ is the volume charge density (charge per unit volume).

### Equivalence of integral and differential forms

Main articles: Divergence theorem
The integral and differential forms are mathematically equivalent, by the divergence theorem. Here is the argument more specifically.

S

{\displaystyle {\scriptstyle \_S}}

E

⋅

d

A

{\displaystyle \mathbf{E} \cdot \mathrm{d}\mathbf{A}}

=

Q

ε
0

{\displaystyle = \frac{Q}{\varepsilon\_0}}

for any closed surface S containing charge Q. By the divergence theorem, this equation is equivalent to:

∭
V

∇
⋅

E

d

# V

Q

ε
0

{\displaystyle \iiint\_V \nabla \cdot \mathbf{E} \, \mathrm{d}V = \frac{Q}{\varepsilon\_0}}

for any volume V containing charge Q. By the relation between charge and charge density, this equation is equivalent to:

∭
V

∇
⋅

E

d

# V

∭
V

ρ

ε
0

d

V

{\displaystyle \iiint\_V \nabla \cdot \mathbf{E} \, \mathrm{d}V = \iiint\_V \frac{\rho}{\varepsilon\_0} \, \mathrm{d}V}

for any volume V. In order for this equation to be *simultaneously true* for *every* possible volume V, it is necessary (and sufficient) for the integrands to be equal everywhere. Therefore, this equation is equivalent to:

∇
⋅

E

=

ρ

ε
0

.

{\displaystyle \nabla \cdot \mathbf{E} = \frac{\rho}{\varepsilon\_0}.}

Thus the integral and differential forms are equivalent.

### Free, bound, and total charge

Main articles: Electric polarization
The electric charge that arises in the simplest textbook situations would be classified as "free charge"—for example, the charge which is transferred in static electricity, or the charge on a capacitor plate. In contrast, "bound charge" arises only in the context of dielectric (polarizable) materials. (All materials are polarizable to some extent.) When such materials are placed in an external electric field, the electrons remain bound to their respective atoms, but shift a microscopic distance in response to the field, so that they're more on one side of the atom than the other. All these microscopic displacements add up to give a macroscopic net charge distribution, and this constitutes the "bound charge".

Although microscopically all charge is fundamentally the same, there are often practical reasons for wanting to treat bound charge differently from free charge. The result is that the more fundamental Gauss's law, in terms of **E** (above), is sometimes put into the equivalent form below, which is in terms of **D** and the free charge only.

### Integral form

This formulation of Gauss's law states the total charge form:

Φ
D

=

Q

f
r
e
e

{\displaystyle \Phi\_D = Q\_\mathrm{free}}

where Φ*D* is the **D**-field flux through a surface S which encloses a volume V, and *Q*free is the free charge contained in V. The flux Φ*D* is defined analogously to the flux Φ*E* of the electric field **E** through S:

Φ
D

=

{\displaystyle \Phi\_D=}

S

{\displaystyle {\scriptstyle \_S}}

D

⋅

d

A

{\displaystyle \mathbf{D} \cdot \mathrm{d}\mathbf{A} }

### Differential form

The differential form of Gauss's law, involving free charge only, states:

∇
⋅

D

=

ρ

f
r
e
e

{\displaystyle \nabla\cdot\mathbf{D} = \rho\_\mathrm{free}}

where ∇ · **D** is the divergence of the electric displacement field, and *ρ*free is the free electric charge density.

## Equivalence of total and free charge statements

∇
⋅

E

=

ρ

ε
0

{\displaystyle \nabla\cdot \mathbf{E} = \dfrac{\rho}{\varepsilon\_0}}

is equivalent to the equation

∇
⋅

D

=

ρ

f
r
e
e

{\displaystyle \nabla\cdot\mathbf{D} = \rho\_\mathrm{free}}

Note that we are only dealing with the differential forms, not the integral forms, but that is sufficient since the differential and integral forms are equivalent in each case, by the divergence theorem.
We introduce the polarization density **P**, which has the following relation to **E** and **D**:

D

=

ε
0

E

P

{\displaystyle \mathbf{D}=\varepsilon\_0 \mathbf{E} + \mathbf{P}}

and the following relation to the bound charge:

ρ

b
o
u
n
d

=
−
∇
⋅

P

{\displaystyle \rho\_\mathrm{bound} = -\nabla\cdot \mathbf{P}}

Now, consider the three equations:

ρ

b
o
u
n
d

=
∇
⋅
(
−

P

)

ρ

f
r
e
e

=
∇
⋅

D

ρ

=
∇
⋅
(

ε
0

E

)

{\displaystyle \begin{align}
\rho\_\mathrm{bound} &= \nabla\cdot (-\mathbf{P}) \\
\rho\_\mathrm{free} &= \nabla\cdot \mathbf{D} \\
\rho &= \nabla \cdot(\varepsilon\_0\mathbf{E})
\end{align}}

The key insight is that the sum of the first two equations is the third equation. This completes the proof: The first equation is true by definition, and therefore the second equation is true if and only if the third equation is true. So the second and third equations are equivalent, which is what we wanted to prove.

## Equation for linear materials

In homogeneous, isotropic, nondispersive, linear materials, there is a simple relationship between **E** and **D**:

D

=
ε

E

{\displaystyle \mathbf{D} = \varepsilon \mathbf{E}}

where ε is the permittivity of the material. For the case of vacuum (aka free space), *ε* = *ε*0. Under these circumstances, Gauss's law modifies to

Φ
E

=

Q

f
r
e
e

ε

{\displaystyle \Phi\_E = \frac{Q\_\mathrm{free}}{\varepsilon}}

for the integral form, and

∇
⋅

E

=

ρ

f
r
e
e

ε

{\displaystyle \nabla \cdot \mathbf{E} = \frac{\rho\_\mathrm{free}}{\varepsilon}}

for the differential form.

## Interpretations

### In terms of fields of force

Gauss's theorem can be interpreted in terms of the lines of force of the field as follows:

The flux through a closed surface is dependent upon both the magnitude and direction of the electric field lines penetrating the surface. In general a positive flux is defined by these lines leaving the surface and negative flux by lines entering this surface. This results in positive charges causing a positive flux and negative charges creating a negative flux. These electric field lines will extend to infinite decreasing in strength by a factor of one over the distance from the source of the charge squared. The larger the number of field lines emanating from a charge the larger the magnitude of the charge is, and the closer together the field lines are the greater the magnitude of the electric field. This has the natural result of the electric field becoming weaker as one moves away from a charged particle, but the surface area also increases so that the net electric field exiting this particle will stay the same. In other words the closed integral of the electric field and the dot product of the derivative of the area will equal the net charge enclosed divided by permittivity of free space.

## Relation to Coulomb's law

### Deriving Gauss's law from Coulomb's law

Strictly speaking, Gauss's law cannot be derived from Coulomb's law alone, since Coulomb's law gives the electric field due to an individual point charge only. However, Gauss's law *can* be proven from Coulomb's law if it is assumed, in addition, that the electric field obeys the superposition principle. The superposition principle says that the resulting field is the vector sum of fields generated by each particle (or the integral, if the charges are distributed smoothly in space).

E

(

r

# )

q

4
π

ε
0

e

r

r
2

{\displaystyle \mathbf{E}(\mathbf{r}) = \frac{q}{4\pi \varepsilon\_0} \frac{\mathbf{e}\_r}{r^2}}

where

**e***r* is the radial unit vector,
r is the radius, |**r**|,
*ε*0 is the electric constant,
q is the charge of the particle, which is assumed to be located at the origin.
Using the expression from Coulomb's law, we get the total field at **r** by using an integral to sum the field at **r** due to the infinitesimal charge at each other point **s** in space, to give

E

(

r

# )

1

4
π

ε
0

∫

ρ
(

s

)
(

r

−

s

)

r

−

s

3

d

3

s

{\displaystyle \mathbf{E}(\mathbf{r}) = \frac{1}{4\pi\varepsilon\_0} \int \frac{\rho(\mathbf{s})(\mathbf{r}-\mathbf{s})}{|\mathbf{r}-\mathbf{s}|^3} \, \mathrm{d}^3 \mathbf{s}}

where ρ is the charge density. If we take the divergence of both sides of this equation with respect to **r**, and use the known theorem

∇
⋅

(

r

r

3

)

=
4
π
δ
(

r

)

{\displaystyle \nabla \cdot \left(\frac{\mathbf{r}}{|\mathbf{r}|^3}\right) = 4\pi \delta(\mathbf{r})}

where *δ*(**r**) is the Dirac delta function, the result is

∇
⋅

E

(

r

# )

1

ε
0

∫
ρ
(

s

)

δ
(

r

−

s

)

d

3

s

{\displaystyle \nabla\cdot\mathbf{E}(\mathbf{r}) = \frac{1}{\varepsilon\_0} \int \rho(\mathbf{s})\, \delta(\mathbf{r}-\mathbf{s})\, \mathrm{d}^3 \mathbf{s}}

Using the "sifting property" of the Dirac delta function, we arrive at

∇
⋅

E

(

r

# )

ρ
(

r

)

ε
0

,

{\displaystyle \nabla\cdot\mathbf{E}(\mathbf{r}) = \frac{\rho(\mathbf{r})}{\varepsilon\_0},}

which is the differential form of Gauss' law, as desired.

Since Coulomb's law only applies to stationary charges, there is no reason to expect Gauss's law to hold for moving charges based on this derivation alone. In fact, Gauss's law does hold for moving charges, and in this respect Gauss's law is more general than Coulomb's law.

Ω
⊆

R
3

{\displaystyle \Omega \subseteq R^3 }

 be a bounded open set, and

E

0

(

r

# )

1

4
π

ε
0

∫

Ω

ρ
(

r

′

)

r

−

r

′

‖

r

−

r

′

‖

3

d

r

′

≡

1

4
π

ε
0

∫

Ω

e
(

r
,

r

′

)

d

r

′

{\displaystyle \mathbf E\_0(\mathbf r) = \frac {1}{4 \pi \varepsilon\_0} \int\_{\Omega} \rho(\mathbf r')\frac {\mathbf r - \mathbf r'} {\left \| \mathbf r - \mathbf r' \right \|^3} \mathop{\mathrm{d}\mathbf{r}'}
\equiv \frac {1}{4 \pi \varepsilon\_0} \int\_{\Omega} e(\mathbf{r, \mathbf{r}'}){\mathrm{d}\mathbf{r}'}}

 be the electric field, with

ρ
(

r

′

)

{\displaystyle \rho(\mathbf r')}

 a continuous function (density of charge).
It is true for all

r

≠

r
′

{\displaystyle \mathbf{r} \neq \mathbf{r'}}

 that

∇

r

⋅

e

(

r
,

r
′

# )
0

{\displaystyle \nabla\_\mathbf{r} \cdot \mathbf{e}(\mathbf{r, r'}) = 0}

.
Consider now a compact set

V
⊆

R
3

{\displaystyle V \subseteq R^3}

 having a piecewise smooth boundary

∂
V

{\displaystyle \partial V}

 such that

Ω
∩
# V
∅

{\displaystyle \Omega \cap V = \emptyset}

. It follows that

e
(

r
,

r

′

)
∈

C
1

(
V
×
Ω
)

{\displaystyle e(\mathbf{r, \mathbf{r}'}) \in C^1(V \times \Omega)}

 and so, for the divergence theorem:

∮

∂
V

E

0

⋅
d

S

=

∫
V

∇

⋅

E

0

d
V

{\displaystyle \oint\_{\partial V} \mathbf{E}\_0 \cdot d\mathbf{S} = \int\_V \mathbf{\nabla} \cdot \mathbf{E}\_0 \, dV}

But because

e
(

r
,

r

′

)
∈

C
1

(
V
×
Ω
)

{\displaystyle e(\mathbf{r, \mathbf{r}'}) \in C^1(V \times \Omega)}

,

∇

⋅

E

0

(

r

# )

{\displaystyle \mathbf{\nabla} \cdot \mathbf{E}\_0(\mathbf{r}) = }

1

4
π

ε
0

∫

Ω

∇

r

⋅
e
(

r
,

r

′

)

d

r

′

{\displaystyle \frac {1}{4 \pi \varepsilon\_0} \int\_{\Omega} \nabla\_\mathbf{r} \cdot e(\mathbf{r, \mathbf{r}'}){\mathrm{d}\mathbf{r}'}}

 = 0 for the argument above (

Ω
∩
# V
∅

⟹

∀

r

∈
V

∀

r
′

∈
Ω

r

≠

r
′

{\displaystyle \Omega \cap V = \emptyset \implies \forall \mathbf{r} \in V \ \ \forall \mathbf{r'} \in \Omega \ \ \ \mathbf{r} \neq \mathbf{r'} }

 and then

∇

r

⋅

e

(

r
,

r
′

# )
0

{\displaystyle \nabla\_\mathbf{r} \cdot \mathbf{e}(\mathbf{r, r'}) = 0}

)
Therefore the flux through a closed surface generated by some charge density outside (the surface) is null.
Now consider

r

0

∈
Ω

{\displaystyle \mathbf{r}\_0 \in \Omega}

, and

B
R

(

r

0

)
⊆
Ω

{\displaystyle B\_R(\mathbf{r}\_0)\subseteq \Omega}

 as the sphere centered in

r

0

{\displaystyle \mathbf{r}\_{0}}

 having

R

{\displaystyle R}

 as radius (it exists because

Ω

{\displaystyle \Omega}

 is an open set).
Let

E

B
R

{\displaystyle \mathbf{E}\_{B\_R}}

 and

E

C

{\displaystyle \mathbf{E}\_C}

 be the electric field created inside and outside the sphere respectively. Then,

E

B
R

{\displaystyle \mathbf{E}\_{B\_R}}

 =

1

4
π

ε
0

∫

B
R

(

r

0

)

e
(

r
,

r

′

)

d

r

′

{\displaystyle \frac {1}{4 \pi \varepsilon\_0} \int\_{B\_R(\mathbf{r}\_0)} e(\mathbf{r, \mathbf{r}'}){\mathrm{d}\mathbf{r}'}}

,

E

C

{\displaystyle \mathbf{E}\_C}

 =

1

4
π

ε
0

∫

Ω
∖

B
R

(

r

0

)

e
(

r
,

r

′

)

d

r

′

{\displaystyle \frac {1}{4 \pi \varepsilon\_0} \int\_{\Omega \setminus B\_R(\mathbf{r}\_0)} e(\mathbf{r, \mathbf{r}'}){\mathrm{d}\mathbf{r}'}}

 and

E

B
R

{\displaystyle \mathbf{E}\_{B\_R}}

E

C

{\displaystyle \mathbf{E}\_C}

 =

E

0

{\displaystyle \mathbf{E}\_0}

Φ
(
R
# )

∮

∂

B
R

(

r

0

)

E

0

⋅
d

S

=

∮

∂

B
R

(

r

0

)

E

B
R

⋅
d

S

∮

∂

B
R

(

r

0

)

E

C

⋅
d

S

=

∮

∂

B
R

(

r

0

)

E

B
R

⋅
d

S

{\displaystyle \Phi(R) =
\oint\_{\partial B\_R(\mathbf{r}\_0)} \mathbf{E}\_0 \cdot d\mathbf{S} =
\oint\_{\partial B\_R(\mathbf{r}\_0)} \mathbf{E}\_{B\_R} \cdot d\mathbf{S} +
\oint\_{\partial B\_R(\mathbf{r}\_0)} \mathbf{E}\_C \cdot d\mathbf{S} =
\oint\_{\partial B\_R(\mathbf{r}\_0)} \mathbf{E}\_{B\_R} \cdot d\mathbf{S} }

The last equality follows by observing that

(
Ω
∖

B
R

(

r

0

)
)
∩

B
R

(

r

0

# )
∅

{\displaystyle (\Omega \setminus B\_R(\mathbf{r}\_0)) \cap B\_R(\mathbf{r}\_0) = \emptyset}

, and the argument above.
The RHS is the electric flux generated by a charged sphere, and so:

Φ
(
R
# )

Q
(
R
)

ε
0

=

1

ε
0

{\displaystyle \Phi(R) =\frac {Q(R)}{\varepsilon\_0} = \frac {1}{\varepsilon\_0} }

∫

B
R

(

r

0

)

ρ
(

r

′

)

d

r

′

=

1

ε
0

ρ
(

r

c
′

)

B
R

(

r

0

)

w
i
t
h

r
c
′

∈

B
R

(

r

0

)

{\displaystyle \int\_{B\_R(\mathbf{r}\_0)} \rho(\mathbf r'){\mathrm{d}\mathbf{r}'} =
\ \ with \ r'\_c \in \ B\_R(\mathbf{r}\_0)}

Where the last equality follows by the mean value theorem for integrals. Using the squeeze theorem and the continuity of

ρ

{\displaystyle \rho}

, one arrives at:

∇

⋅

E

0

(

r

0

# )

lim

R
→
0

1

B
R

(

r

0

)

Φ
(
R
# )

1

ε
0

ρ
(

r

0

)

{\displaystyle \mathbf{\nabla} \cdot \mathbf{E}\_0(\mathbf{r}\_0) =
\lim\_{R \to 0} \frac{1}{|B\_R(\mathbf{r}\_0)|}\Phi(R) =
\frac {1}{\varepsilon\_0} \rho(\mathbf r\_0) }

### Deriving Coulomb's law from Gauss's law

Strictly speaking, Coulomb's law cannot be derived from Gauss's law alone, since Gauss's law does not give any information regarding the curl of **E** (see Helmholtz decomposition and Faraday's law). However, Coulomb's law *can* be proven from Gauss's law if it is assumed, in addition, that the electric field from a point charge is spherically symmetric (this assumption, like Coulomb's law itself, is exactly true if the charge is stationary, and approximately true if the charge is in motion).

∮
S

E

⋅
d

A

=

Q

ε
0

{\displaystyle \oint\_S\mathbf{E}\cdot d\mathbf{A} = \frac{Q}{\varepsilon\_0} }

By the assumption of spherical symmetry, the integrand is a constant which can be taken out of the integral. The result is

4
π

r
2

r

^

⋅

E

(

r

# )

Q

ε
0

{\displaystyle 4\pi r^2\hat{\mathbf{r}}\cdot\mathbf{E}(\mathbf{r}) = \frac{Q}{\varepsilon\_0}}

where **r̂** is a unit vector pointing radially away from the charge. Again by spherical symmetry, **E** points in the radial direction, and so we get

E

(

r

# )

Q

4
π

ε
0

r

^

r
2

{\displaystyle \mathbf{E}(\mathbf{r}) = \frac{Q}{4\pi \varepsilon\_0} \frac{\hat{\mathbf{r}}}{r^2}}

which is essentially equivalent to Coulomb's law. Thus the inverse-square law dependence of the electric field in Coulomb's law follows from Gauss' law.

## See also

* Method of image charges
* Uniqueness theorem for Poisson's equation

## Notes

1. ↑ More specifically, the infinitesimal area is thought of as planar and with area d*A*. The vector d**A** is normal to this area element and has magnitude d*A*.