docs/physics/electromagnetism/electrostatics.md: update

docs/physics/electromagnetism/electrostatics.md

@@ -2,7 +2,7 @@
 
 A Notational remark, let the vector from origin to source point be given by $\mathbf{r}'$ and the vector from origin to field point by $\mathbf{r}$. The vector from source to field point is given by $\bm{\mathfrak{r}} = \mathbf{r} - \mathbf{r}'$.
 
-Let us start with electrostatics in vacuum, where no current, no magnetic field and no time dependence are taken into account. We call this the **electrostatic regime**.
+In the **electrostatic regime** there is no current, no magnetic fields and time dependence is not taken into account
 
 ## Electrostatics in vacuum
 
@@ -12,16 +12,16 @@ Electrostatics builds entirely on the following axiom:
 > 
 > $$
 >   \begin{align*}
->       \nabla \cdot \mathbf{E} &= \frac{\rho}{\varepsilon_0},\\
+>       \nabla \cdot \mathbf{E} &= \frac{\rho}{\epsilon_0},\\
 >       \nabla \times \mathbf{E} &= \mathbf{0},
 >   \end{align*}
 > $$
 >
-> with $\rho: \mathbf{r} \mapsto \rho(\mathbf{r})$ the **space charge density** and $\varepsilon_0$ the **permittivity of space**. 
+> with $\rho: \mathbf{r} \mapsto \rho(\mathbf{r})$ the **volume charge density** and $\epsilon_0$ the **vacuum permittivity**. 
 
 The following definition connects the electric field $\mathbf{E}: \mathbf{r} \mapsto \mathbf{E}(\mathbf{r})$ to the Newtonian formalism of mechanics:
 
-> *Definition 1*: The force field $\mathbf{F}$ on a point charge $q$ at a field point $\mathbf{r}$ is defined as
+> *Definition 1*: The electric force $\mathbf{F}$ on a point charge $q$ at a field point $\mathbf{r}$ is defined as
 >
 > $$
 >   \mathbf{F}(\mathbf{r}) = q \mathbf{E}(\mathbf{r}),
@@ -34,7 +34,7 @@ From *Axiom 1* we may proof that the following theorem holds (Coulomb law):
 > *Theorem 1*: The electric field $\mathbf{E}: \mathbf{r} \mapsto \mathbf{E}(\mathbf{r})$ may be described as
 >
 > $$
->   \mathbf{E}(\mathbf{r}) = \begin{cases} \frac{1}{4\pi\varepsilon_0} \cdot \frac{q}{\mathfrak{r}^2} \mathbf{e}_\mathfrak{r} \qquad &\mathrm{0D},\\ \frac{1}{4\pi\varepsilon_0} \int \frac{\lambda}{\mathfrak{r}^2} \mathbf{e}_\mathfrak{r} dl \qquad &\mathrm{1D},\\ \frac{1}{4\pi\varepsilon_0} \iint \frac{\sigma}{\mathfrak{r}^2} \mathbf{e}_\mathfrak{r} da \qquad &\mathrm{2D},\\ \frac{1}{4\pi\varepsilon_0} \iiint \frac{\rho}{\mathfrak{r}^2} \mathbf{e}_\mathfrak{r} d\tau \qquad &\mathrm{3D}, \end{cases}
+>   \mathbf{E}(\mathbf{r}) = \begin{cases} \frac{1}{4\pi\epsilon_0} \cdot \frac{q}{\mathfrak{r}^2} \mathbf{e}_\mathfrak{r} \qquad &\mathrm{0D},\\ \frac{1}{4\pi\epsilon_0} \int \frac{\lambda}{\mathfrak{r}^2} \mathbf{e}_\mathfrak{r} dl' \qquad &\mathrm{1D},\\ \frac{1}{4\pi\epsilon_0} \iint \frac{\sigma}{\mathfrak{r}^2} \mathbf{e}_\mathfrak{r} da' \qquad &\mathrm{2D},\\ \frac{1}{4\pi\epsilon_0} \iiint \frac{\rho}{\mathfrak{r}^2} \mathbf{e}_\mathfrak{r} d\tau' \qquad &\mathrm{3D}, \end{cases}
 > $$
 >
 > with $\lambda$ and $\sigma$ the line and surface charge density.
@@ -46,7 +46,7 @@ From *Axiom 1* we may proof that the following theorem holds (Coulomb law):
 From the divergence theorem it follows as well that 
 
 $$
-    \oint_S \mathbf{E} \cdot d\mathbf{a} = \int_V \nabla \cdot \mathbf{E} d\tau = \int_V \frac{\rho}{\varepsilon_0} d\tau = \frac{Q}{\varepsilon_0},
+    \oint_S \mathbf{E} \cdot d\mathbf{a} = \int_V \nabla \cdot \mathbf{E} d\tau = \int_V \frac{\rho}{\epsilon_0} d\tau = \frac{Q}{\epsilon_0},
 $$
 
 with $Q$ the enclosed charge. This result is called Gauß' law. 
@@ -75,14 +75,14 @@ One may simply proof that from *Definition 2* it follows that:
 
     Follows from gradient theorem.
 
-Now from (*Axiom 1*) $\nabla \cdot \mathbf{E} = \frac{\rho}{\varepsilon_0}$ it follows that $\nabla^2 V = - \frac{\rho}{\varepsilon_0}$.
+Now from (*Axiom 1*) $\nabla \cdot \mathbf{E} = \frac{\rho}{\epsilon_0}$ it follows that $\nabla^2 V = - \frac{\rho}{\epsilon_0}$.
 
 Similarly to *Theorem 1* we may state the following about the electric potential $V$:
 
 > *Theorem 3*: The electric potential $V: \mathbf{r} \mapsto V(\mathbf{r})$ may be described as
 >
 > $$
->   V(\mathbf{r}) = \begin{cases} \frac{1}{4\pi\varepsilon_0} \cdot \frac{q}{\mathfrak{r}} \qquad &\mathrm{0D},\\ \frac{1}{4\pi\varepsilon_0} \int \frac{\lambda}{\mathfrak{r}} dl \qquad &\mathrm{1D},\\ \frac{1}{4\pi\varepsilon_0} \iint \frac{\sigma}{\mathfrak{r}} da \qquad &\mathrm{2D},\\ \frac{1}{4\pi\varepsilon_0} \iiint \frac{\rho}{\mathfrak{r}} d\tau \qquad &\mathrm{3D}.\end{cases}
+>   V(\mathbf{r}) = \begin{cases} \frac{1}{4\pi\epsilon_0} \cdot \frac{q}{\mathfrak{r}} \qquad &\mathrm{0D},\\ \frac{1}{4\pi\epsilon_0} \int \frac{\lambda}{\mathfrak{r}} dl' \qquad &\mathrm{1D},\\ \frac{1}{4\pi\epsilon_0} \iint \frac{\sigma}{\mathfrak{r}} da' \qquad &\mathrm{2D},\\ \frac{1}{4\pi\epsilon_0} \iiint \frac{\rho}{\mathfrak{r}} d\tau' \qquad &\mathrm{3D}.\end{cases}
 > $$
 
 ??? note "Proof:"
@@ -120,7 +120,7 @@ As we set the reference point at infinity, then $qV$ could be interpreted as the
     $$
         \begin{align*}
             W_i &= q_i \sum_j V_j(\mathbf{r}_i),\\ 
-                &= q_i \sum_j \frac{1}{4\pi\varepsilon_0} \frac{q_j}{\mathfrak{r}_{ij}},
+                &= q_i \sum_j \frac{1}{4\pi\epsilon_0} \frac{q_j}{\mathfrak{r}_{ij}},
         \end{align*}
     $$
 
@@ -129,14 +129,14 @@ As we set the reference point at infinity, then $qV$ could be interpreted as the
     $$
         \begin{align*}
             W &= \sum_{i=1}^N W_i,\\
-              &= \frac{1}{2} \sum_{i=1}^N \sum_{j \neq i} \frac{1}{4\pi\varepsilon_0} \frac{q_j}{\mathfrak{r}_{ij}},\\
+              &= \frac{1}{2} \sum_{i=1}^N \sum_{j \neq i} \frac{1}{4\pi\epsilon_0} \frac{q_j}{\mathfrak{r}_{ij}},\\
               &=\frac{1}{2} \sum_{i=1}^N q_i V(\mathbf{r}_i).
         \end{align*}
     $$
 
 Is the energy stored in the system.
 
-It follows now that for a space charge density $\rho$ the work $W$ needed to construct the system is given by
+It follows now that for a volume charge density $\rho$ the work $W$ needed to construct the system is given by
 
 $$
     W = \frac{1}{2} \int_\mathscr{V} \rho V d\tau,
@@ -147,7 +147,7 @@ and the integrals for line and surface charge densities ($\lambda,\sigma$) are o
 > *Theorem 5*: The work $W$ required to construct the system can be expressed in terms of the electric field $\mathbf{E}$ as
 >
 > $$
->   W = \frac{\varepsilon_0}{2} \int_{\mathbb{R}^3} \|\mathbf{E}\|^2 d\tau.
+>   W = \frac{\epsilon_0}{2} \int_{\mathbb{R}^3} \|\mathbf{E}\|^2 d\tau.
 > $$
 
 ??? note "Proof:"
@@ -156,16 +156,16 @@ and the integrals for line and surface charge densities ($\lambda,\sigma$) are o
 
     $$
         \begin{align*}
-            W &= \frac{\varepsilon_0}{2} \int_\mathscr{V} (\nabla \cdot \mathbf{E}) V d\tau,\\
+            W &= \frac{\epsilon_0}{2} \int_\mathscr{V} (\nabla \cdot \mathbf{E}) V d\tau,\\
-              &= \frac{\varepsilon_0}{2} \Big(-\int_\mathscr{V} \mathbf{E} \cdot (\nabla V) d\tau + \oint_{\partial \mathscr{V}} V \mathbf{E} \cdot d\mathbf{a}\Big),\\
+              &= \frac{\epsilon_0}{2} \Big(-\int_\mathscr{V} \mathbf{E} \cdot (\nabla V) d\tau + \oint_{\partial \mathscr{V}} V \mathbf{E} \cdot d\mathbf{a}\Big),\\
-              &= \frac{\varepsilon_0}{2} \Big(\int_\mathscr{V} \|\mathbf{E}\|^2 d\tau + \oint_\mathscr{\partial \mathscr{V}} V \mathbf{E} \cdot d\mathbf{a}\Big).
+              &= \frac{\epsilon_0}{2} \Big(\int_\mathscr{V} \|\mathbf{E}\|^2 d\tau + \oint_\mathscr{\partial \mathscr{V}} V \mathbf{E} \cdot d\mathbf{a}\Big).
         \end{align*}
     $$
 
     If we now set $\mathscr{V} = \mathbb{R}^3$ then the integral over $\partial \mathscr{V}$ goes to zero and we are left with:
 
     $$
-        W = \frac{\varepsilon_0}{2} \int_{\mathbb{R}^3} \|\mathbf{E}\|^2 d\tau.
+        W = \frac{\epsilon_0}{2} \int_{\mathbb{R}^3} \|\mathbf{E}\|^2 d\tau.
     $$
 
 That is performing the integral over all of space. Which is mathematically rather nice, but imposes that the construction of point charges requires an infinite amount of energy. A result of the introduction of infinities. 
@@ -177,11 +177,11 @@ As may be observed, this result does not obey the superposition principle.
 In the electrostatic regime the following properties of conductors are valid:
 
 1. $\mathbf{0} = \mathbf{J} = \sigma \mathbf{E} \implies \mathbf{E} = \mathbf{0}$ inside a conductor.
-2. $\mathbf{0} = \nabla \cdot \mathbf{E} = \frac{\rho}{\varepsilon_0} \implies \rho = 0$ inside a conductor.
+2. $\mathbf{0} = \nabla \cdot \mathbf{E} = \frac{\rho}{\epsilon_0} \implies \rho = 0$ inside a conductor.
 
 Which implies that any net charge resides on the surface, that a conductor is an equipotential and that the electric field is normal to the surface of the conductor.
 
-Since the electric field $\mathbf{E}$ is proportional to the space charge density $\rho$ or charge $Q$, so also is $V$. Which enables us to define a constant of proportionality $C = \frac{Q}{V}$, the capacitance. 
+Since the electric field $\mathbf{E}$ is proportional to the volume charge density $\rho$ or charge $Q$, so also is $V$. Which enables us to define a constant of proportionality $C = \frac{Q}{V}$, the capacitance. 
 
 The work required to charge a capacitor can then be expressed as
 
@@ -194,7 +194,7 @@ $$
 The primary task of electrostatics is to find the electric field $\mathbf{E}$ of a given stationary charge distribution $\rho$. In principle, this purpose is accomplished by Coulombs's law as by the result of *Theorem 1*. But these integrals are generally difficult to solve for any but the simplest charge distributions. In some case symmetries in the charge distribution allow Gauß' law to be used. But more often than not we require to determine the potential $V$ with the result in *Theorem 3* or with Poisson's equation:
 
 $$
-    \nabla^2 V = -\frac{\rho}{\varepsilon_0},
+    \nabla^2 V = -\frac{\rho}{\epsilon_0},
 $$
 
 where in the region where $\rho = 0$ this reduces to Laplace's equation:
@@ -223,13 +223,13 @@ This result brings us to the consideration of the boundary condition for which t
 
 and consequentially:
 
-> *Corollary 2*: An electric potential $V: \mathbf{r} \mapsto V(\mathbf{r})$ that satisfies $\nabla^2 V=-\frac{\rho}{\varepsilon_0}$ in some volume $\mathscr{V}$ is unique if $\rho: \mathbf{r} \mapsto \rho(\mathbf{r})$ is known $\forall \mathbf{r} \in \mathscr{V}$ and $V$ is specified on the boundary surface $\partial \mathscr{V}$.
+> *Corollary 2*: An electric potential $V: \mathbf{r} \mapsto V(\mathbf{r})$ that satisfies $\nabla^2 V=-\frac{\rho}{\epsilon_0}$ in some volume $\mathscr{V}$ is unique if $\rho: \mathbf{r} \mapsto \rho(\mathbf{r})$ is known $\forall \mathbf{r} \in \mathscr{V}$ and $V$ is specified on the boundary surface $\partial \mathscr{V}$.
 
 The combination of *Theorem 7* and *Corollary 2* is known as the **first uniqueness theorem**.
 
 The **second uniqueness theorem** is relevant when the boundary is defined by conductors:
 
-> *Theorem 8*: An electric potential $V: \mathbf{r} \mapsto V(\mathbf{r})$ that satisfies $\nabla^2 V = -\frac{\rho}{\varepsilon_0}$ in some volume $\mathscr{V}$ that is surrounded by conductors is unique if $\rho: \mathbf{r} \mapsto \rho(\mathbf{r})$ is known $\forall \mathbf{r} \in \mathscr{V}$ and the total charge on the conductors is known.
+> *Theorem 8*: An electric potential $V: \mathbf{r} \mapsto V(\mathbf{r})$ that satisfies $\nabla^2 V = -\frac{\rho}{\epsilon_0}$ in some volume $\mathscr{V}$ that is surrounded by conductors is unique if $\rho: \mathbf{r} \mapsto \rho(\mathbf{r})$ is known $\forall \mathbf{r} \in \mathscr{V}$ and the total charge on the conductors is known.
 
 The method of images to solve for the electric potential $V$ makes use of these uniqueness theorems. In essence this method uses a solveable image system that obeys the same boundary conditions, then by unicity the same solution should be obtained. Though, the electric work required to construct the system is not always the same.
 
@@ -242,7 +242,7 @@ Approximation is the de facto option if no exact solution of the electric potent
 > *Theorem 9*: The electric potential $V: \mathbf{r} \mapsto V(\mathbf{r})$ of a charge distribution $\rho: \mathbf{r} \mapsto \rho(\mathbf{r})$ may be described as
 >
 > $$
->   V(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0} \sum_{n=0}^\infty \frac{1}{r^{n+1}} \iiint (r')^n P_n(\cos \alpha) \rho(\mathbf{r}') d\tau',
+>   V(\mathbf{r}) = \frac{1}{4\pi\epsilon_0} \sum_{n=0}^\infty \frac{1}{r^{n+1}} \iiint (r')^n P_n(\cos \alpha) \rho(\mathbf{r}') d\tau',
 > $$
 >
 > with $\cos \alpha = \frac{\langle \mathbf{r}, \mathbf{r}'\rangle}{r r'}$.
@@ -252,7 +252,7 @@ Approximation is the de facto option if no exact solution of the electric potent
     We start with
 
     $$
-        V(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0} \iiint \frac{\rho(\mathbf{r}')}{\mathfrak{r}} d\tau'.
+        V(\mathbf{r}) = \frac{1}{4\pi\epsilon_0} \iiint \frac{\rho(\mathbf{r}')}{\mathfrak{r}} d\tau'.
     $$
 
     We may write $\mathfrak{r}$ as
@@ -282,13 +282,13 @@ Approximation is the de facto option if no exact solution of the electric potent
     which obtains
 
     $$
-        V(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0} \sum_{n=0}^\infty \frac{1}{r^{n+1}} \iiint (r')^n P_n(\cos \alpha) \rho(\mathbf{r}') d\tau'.
+        V(\mathbf{r}) = \frac{1}{4\pi\epsilon_0} \sum_{n=0}^\infty \frac{1}{r^{n+1}} \iiint (r')^n P_n(\cos \alpha) \rho(\mathbf{r}') d\tau'.
     $$
 
 The monopole $(n=0)$ term of $V$ is then given by
 
 $$
-    V_0(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0} \frac{1}{r} \int \rho(\mathbf{r}') d\tau',
+    V_0(\mathbf{r}) = \frac{1}{4\pi\epsilon_0} \frac{1}{r} \iiint \rho(\mathbf{r}') d\tau',
 $$
 
 is zero if the net charge is zero.
@@ -297,24 +297,24 @@ The dipole $(n=1)$ term of $V$ is given by
 
 $$
     \begin{align*}
-        V_1(\mathbf{r}) &= \frac{1}{4\pi\varepsilon_0} \cdot \frac{1}{r^2} \langle \frac{1}{r} \mathbf{r}, \int \mathbf{r}' \rho(\mathbf{r}') d\tau' \rangle,\\
+        V_1(\mathbf{r}) &= \frac{1}{4\pi\epsilon_0} \cdot \frac{1}{r^2} \Bigg\langle \frac{1}{r} \mathbf{r}, \iiint \mathbf{r}' \rho(\mathbf{r}') d\tau' \Bigg\rangle,\\
-        &= \frac{1}{4\pi\varepsilon_0} \cdot \frac{1}{r^2} \langle \mathbf{e}_r, \mathbf{p} \rangle,
+        &= \frac{1}{4\pi\epsilon_0} \cdot \frac{1}{r^2} \langle \mathbf{e}_r, \mathbf{p} \rangle,
     \end{align*}
 $$
 
-with $\mathbf{p} = \int \mathbf{r}' \rho(\mathbf{r}') d\tau'$ the **dipole moment**.
+with $\mathbf{p} = \iiint \mathbf{r}' \rho(\mathbf{r}') d\tau'$ the **electric dipole moment**.
 
-For a collection of $N \in \mathbb{N}$ point charges we have that the dipole moment $\mathbf{p}$ is given by
+For a collection of $N \in \mathbb{N}$ point charges we have that the electric dipole moment $\mathbf{p}$ is given by
 
 $$
     \mathbf{p} = \sum_{i=1}^N q_i \mathbf{r}_i'.
 $$
 
-Note that if the net charge of the charge distribution is zero then the dipole moment is independent of the choice of origin. 
+Note that if the net charge of the charge distribution is zero then the electric dipole moment is independent of the choice of origin. 
 
 ## Electrostatics in matter
 
-The dipole moment $\mathbf{p}$ of an induced dipole due to an external field $\mathbf{E}$ is given by
+The electric dipole moment $\mathbf{p}$ of an induced dipole due to an external field $\mathbf{E}$ is given by
 
 $$
     \mathbf{p} = \bm{\alpha} \lrcorner \mathbf{E},
@@ -322,7 +322,7 @@ $$
 
 with $\bm{\alpha} \in \mathscr{T}^2_0$ the polarisability.
 
-> *Definition 3*: The induced dipole moment $\mathbf{p}$ in a medium may be expressed by the **polarisation** $\mathbf{P}$ of the medium defined in terms of
+> *Definition 3*: The induced electric dipole moment $\mathbf{p}$ in a medium may be expressed by the **polarisation** $\mathbf{P}$ of the medium defined in terms of
 >
 > $$
 >   \mathbf{p} = \int_\mathscr{V} \mathbf{P} d\tau,
@@ -330,7 +330,7 @@ with $\bm{\alpha} \in \mathscr{T}^2_0$ the polarisability.
 >
 > with $\mathscr{V}$ the volume of the medium.
 
-The polarisation $\mathbf{P}$ is in essence a sort of dipole moment per unit volume of the medium.
+The polarisation $\mathbf{P}$ is in essence a sort of electric dipole moment per unit volume of the medium.
 
 > *Theorem 10*: The polarisation of the medium $\mathbf{P}$ adheres to
 >
@@ -341,42 +341,48 @@ The polarisation $\mathbf{P}$ is in essence a sort of dipole moment per unit vol
 >   \end{align*}
 > $$
 >
-> with $\rho_b: \mathbf{r} \mapsto \rho_b(\mathbf{r})$ the **bound space charge density** and $\sigma_b: \mathbf{r} \mapsto \sigma_b(\mathbf{r})$ the **bound surface charge density**.
+> with $\rho_b: \mathbf{r} \mapsto \rho_b(\mathbf{r})$ the **bound volume charge density** and $\sigma_b: \mathbf{r} \mapsto \sigma_b(\mathbf{r})$ the **bound surface charge density**.
 
 ??? note "Proof:"
 
     The dipole term of the potential in terms of the polarisation $\mathbf{P}$ is given by
 
     $$
-        V(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0} \int_\mathscr{V} \frac{\mathbf{P}(\mathbf{r}') \cdot \mathbf{e}_\mathfrak{r}}{\mathfrak{r}^2} d\tau'.
+        V(\mathbf{r}) = \frac{1}{4\pi\epsilon_0} \int_\mathscr{V} \frac{\mathbf{P}(\mathbf{r}') \cdot \mathbf{e}_\mathfrak{r}}{\mathfrak{r}^2} d\tau'.
     $$
     
-    Now observe that $\nabla' \frac{1}{\mathfrak{r}} = \frac{1}{\mathfrak{r}^2} \mathbf{e}_\mathfrak{r}, obtains
+    Now observe that $\nabla' \frac{1}{\mathfrak{r}} = \frac{1}{\mathfrak{r}^2} \mathbf{e}_\mathfrak{r}$, obtains
 
     $$
-        V(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0} \int_\mathscr{V} \mathbf{P}(\mathbf{r}') \cdot \nabla' \frac{1}{\mathfrak{r}} d\tau',
+        V(\mathbf{r}) = \frac{1}{4\pi\epsilon_0} \int_\mathscr{V} \mathbf{P}(\mathbf{r}') \cdot \nabla' \frac{1}{\mathfrak{r}} d\tau',
+    $$
+
+    and by integration by parts
+
+    $$
+        V(\mathbf{r}) = \frac{1}{4\pi\epsilon_0} \Bigg(\int_\mathscr{V} \nabla' \cdot \frac{1}{\mathfrak{r}} \mathbf{P}(\mathbf{r}') d\tau' - \int_\mathscr{V} \frac{1}{\mathfrak{r}} \nabla' \cdot \mathbf{P}(\mathbf{r}') d\tau'\Bigg),
     $$
 
     such that by the divergence theorem we obtain
 
     $$
-        V(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0} \Bigg(\oint_{\partial \mathscr{V}} \frac{1}{\mathfrak{r}} \mathbf{P}(\mathbf{r}') \cdot d\mathbf{a}' - \int_\mathscr{V} \frac{1}{\mathfrak{r}} \nabla' \cdot \mathbf{P}(\mathbf{r}') d\tau' \Bigg).
+        V(\mathbf{r}) = \frac{1}{4\pi\epsilon_0} \Bigg(\oint_{\partial \mathscr{V}} \frac{1}{\mathfrak{r}} \mathbf{P}(\mathbf{r}') \cdot d\mathbf{a}' - \int_\mathscr{V} \frac{1}{\mathfrak{r}} \nabla' \cdot \mathbf{P}(\mathbf{r}') d\tau' \Bigg).
     $$
 
     Setting $\rho_b = - \nabla \cdot \mathbf{P}$ and $\sigma_b = \mathbf{P} \cdot \mathbf{e}_n$ we obtain
 
     $$
-        V(\mathbf{r}) = \frac{1}{4\pi\varepsilon_0} \Bigg(\oint_{\partial \mathscr{V}} \frac{\sigma_b(\mathbf{r}')}{\mathfrak{r}} da' + \int_\mathscr{V} \frac{\rho_b(\mathbf{r}')}{\mathfrak{r}} d\tau'\Bigg),
+        V(\mathbf{r}) = \frac{1}{4\pi\epsilon_0} \Bigg(\oint_{\partial \mathscr{V}} \frac{\sigma_b(\mathbf{r}')}{\mathfrak{r}} da' + \int_\mathscr{V} \frac{\rho_b(\mathbf{r}')}{\mathfrak{r}} d\tau'\Bigg),
     $$
 
-    implying that the potential of a polarised object is the same as that produced by a space charge density $\rho_b$ plus a a surface charge density $\sigma_b$. 
+    implying that the potential of a polarised object is the same as that produced by a volume charge density $\rho_b$ plus a a surface charge density $\sigma_b$. 
 
-We may as well define the **free space charge density** $\rho_f$ in terms of the space charge density $\rho$ and the bound space charge density $\rho_b$ by notion of $\rho_f = \rho - \rho_b$. 
+We may as well define the **free volume charge density** $\rho_f$ in terms of the volume charge density $\rho$ and the bound volume charge density $\rho_b$ by notion of $\rho_f = \rho - \rho_b$. 
 
 > *Definition 4*: Let the electric displacement $\mathbf{D}$ of the medium be defined as
 >
 > $$
->   \mathbf{D} = \varepsilon_0 \mathbf{E} + \mathbf{P},
+>   \mathbf{D} = \epsilon_0 \mathbf{E} + \mathbf{P},
 > $$
 >
 > with $\mathbf{E}$ the electric field and $\mathbf{P}$ the polarisation of the medium.
@@ -392,30 +398,38 @@ The usefullness of *Definition 4* may become apparent with the following result:
 >   \end{align*}
 > $$
 >
-> with $\rho_f$ the free space charge density of the medium and $\mathbf{P}$ the polarisation of the medium.
+> with $\rho_f$ the free volume charge density of the medium and $\mathbf{P}$ the polarisation of the medium.
 
 ??? note "Proof:"
     
     From *Definition 4* we may write the electric field in terms of the electric displacement $\mathbf{D}$ and the polarisation $\mathbf{P}$ of the medium
 
     $$
-        \mathbf{E} = \frac{1}{\varepsilon_0} \bigg(\mathbf{D} - \mathbf{P}\bigg),
+        \mathbf{E} = \frac{1}{\epsilon_0} \bigg(\mathbf{D} - \mathbf{P}\bigg),
     $$
 
     such that with *Axiom 1* we obtain
 
     $$
-        \varepsilon_0 \nabla \cdot \mathbf{E} = \nabla \cdot \mathbf{D} - \rho_b \implies \nabla \cdot \mathbf{D} = \rho_f,
+        \epsilon_0 \nabla \cdot \mathbf{E} = \nabla \cdot \mathbf{D} + \rho_b \implies \nabla \cdot \mathbf{D} = \rho_f,
     $$
 
     and
 
     $$
-        \varepsilon_0 \nabla \times \mathbf{E} = \nabla \times \mathbf{D} - \nabla \times \mathbf{P} \implies \nabla \times \mathbf{D} = \nabla \times \mathbf{P}.
+        \epsilon_0 \nabla \times \mathbf{E} = \nabla \times \mathbf{D} - \nabla \times \mathbf{P} \implies \nabla \times \mathbf{D} = \nabla \times \mathbf{P}.
     $$
 
+Recall that from the divergence theorem we have
+
+$$
+    \oint_{\partial \mathscr{V}} \mathbf{D} \cdot d\mathbf{a} = \int_\mathscr{V} \rho_f d\tau = Q_f,
+$$
+
+with $Q_f$ the enclosed free charge, if $\nabla \times \mathbf{D} = \mathbf{0}$.
+
 ### Linear media
 
-In linear media we have $\mathbf{P} = \varepsilon_0 \chi_e \mathbf{E}$ with $\chi_e$ the electric susceptibility of the medium. 
+In linear media we have $\mathbf{P} = \epsilon_0 \chi_e \mathbf{E}$ with $\chi_e$ the electric susceptibility of the medium. 
 
-Furthermore $\mathbf{D} = \varepsilon_0 (1 + \chi_e) \mathbf{E} = \varepsilon \mathbf{E}$ with $\varepsilon$ the permittivity and $\varepsilon_r = 1 + \chi_e = \frac{\varepsilon}{\varepsilon_0}$ the dielectric constant.
+Furthermore $\mathbf{D} = \epsilon_0 (1 + \chi_e) \mathbf{E} = \epsilon \mathbf{E}$ with $\epsilon$ the permittivity and $\epsilon_r = 1 + \chi_e = \frac{\epsilon}{\epsilon_0}$ the dielectric constant.