# Electrostatics and Magnetostatics¶

## Electric Field¶

The electric potential energy of a point charge in an electric field is given by

In the case of fixed point charge \(Q\) field, the amount electric energy of test charge \(q\) is

Suppose we have N point charges places randomly in the space, the total electric energy of the system can be calculated through the a superposition of the electric potential.

in which the half is due to the double counting of the interactions.

### Energy Density¶

The next question we are going to consider is the electric potential energy of a bulk charged object.

where

Derivation of Electric Energy of Charged Object

Construct the system by adding small amount of charge to it from zero charge distribution.

Charge distribution \(\rho(\vec x)\) can be changed linearly \(\lambda \rho(\vec x')\) where \(\lambda\) is a scalar number.

The energy change due to adding charge to the system to change the charge distribution linearly \(\delta \lambda\). Energy change is

where

Integrate from 0 to 1 will give us the total energy for an object with charge distribution \(\rho(\vec x)\).

It is obvious that the \(1/2\) comes from the linear dependence of \(\lambda \rho(\vec x)\).

Now we have the expression for energy of electric field, it is straightforward to find the energy density.

- To calculate that, we need to put together two parts of energy, i.e., E-field outside of the conductor and the surface electric potential energy of the conductor.
- Then apple the divergence theorem, counting the infinite surface and also all the surfaces of the conductor.
- Notice the surface integral of the potential energy on the conductor surface is exactly the negative of the surface energy density we put in to the expression in the first step.
- The surface integral at infinity surface is zero.
- The only term left is \(-\frac{1}{8\pi}\int d^3x E_i \partial_i \phi=\frac{1}{8\pi}\int d^3x E_i E_i\).

### Electric Forces¶

Alright, move on to the concept of force in E&M. Coloumb force is

Maxwell’s equations tells that

So force can be rewritten as

Now notice that \(\partial_i(E_j E_j)=2E_j\partial_iE_j\) and \(\partial_i E_j = \partial_j E_i\).

Proof

So \(\nabla\times \vec E = 0\) means \(\partial_i E_j - \partial_J E_i = 0\).

Force becomes

Recall that force is momentum per unit time. What is inside the integral means some force density or momentum flow density per unit time, for some reason we use the definition

which is symmetric under the exchange of i and j.

Stress Tensor in General Relativity or Fluid Dynamics

Electromagnetic energy momentum tensor using Gaussian units and +2 signature is

where \(\mathbf F = \mathbf d \mathbf A \Rightarrow F_{\mu\nu} = \partial_\mu A_\nu - \partial_\nu A_\mu\).

To change the index to \(F^\mu_{\phantom{\mu}\nu}\) we just use the Minkowski metric which just put plus and minus signs on the components,

Why is this the case? This is because of the Lorentz group requires it. Generators of \(L_+^\uparrow\) which means the orthochronous patch combine together with the infinitesimal change to form the tiny change in Lie algebra, that is to form \(\omega\) in \(L = I +\omega\)

where the six matrices are basically the generators to construct the Lorentz group.

Energy momentum tensor can be decomposed using the generators,

where we are using the Einstein summation rule.

For more information about this please read arXiv:physics/0005084 .

Force on Capacitor

Suppose we have a capacitor with two parallel round plates of radius \(R\) and separation \(d\). The top plate has charge \(Q\) and bottom plate has charge \(-Q\). Now ask the question, what is the force on the top plate?

Is its magnitude \(F=QE\) because of the fact that the charge \(Q\) is in a electric field \(\vec E\)? NO.

We can calculate the force using several methods like

- principle of virtual work,
- stress tensor.

However the result shows that the magnitude of force is only \(F=\frac{1}{2} Q E\), half of the expectation we had at first.

The intuition is that electric field only exists on one side of the plate while \(F=qE\) describes the force of a charged particle emerged in the electric field. Once one knows this, the half comes from the fact that stress only applies on one side of the top plate.

### Multipole Expansion¶

For a charge object with charge distribution \(\rho(\vec x)\), the electic potential is

For any general distribution, we can always have \(\frac{1}{\lvert \vec x - \vec x' \rvert}\) expanded and define differential multipoles.

Monopole and Dipole

Before we go into this expansion, a review of the idea of monopole and dipole is very useful.

Monopole is the case of spherical symmetric distribution of electric potential.

Dipole, using the simplest model, is a system of two charged particles with \(\pm Q\) and directed from negative charge to positive charge.

In the language of math, \(\vec q=Q\vec d\). The electric potential is the superposition of the two charges.

I need some time to make clear of different methods.

Dipole comes from the expansion at \(\vec x \gg \vec x'\),

thus the dipole part is actually

Multipole

A nth multipole has the dimnsion of \(\frac{[Q]}{[r]} \left(\frac{d}{r}\right)^2\).

In most of the cases, we have a very small dipole compared to the distance at where the field is measured. So we could **take the limit that the dipole is point like**. The actually meaning of a point dipole is that the dipole is a constant viewed from very far away.

The next question to be answered is the electric field generated by this point dipole.

which is **only true for** \(\vec x \neq 0\).

Force feels by a dipole is

Torque can also be calculated

Force and Torque of Dipole

Force can be calculated using principle of virtual work. \(-\vec p \cdot \vec E\) is the energy which can be explained using simple two charge model.

Torque has two parts, one with a precession-like term which is called spin torque while another one is a orbital torque since it tries to align the dipole with the field lines.

## Dielectric Material¶

In the presence of dielectric material, a new set of quantities will be introduced.

### Polarization¶

Why do we introduce this polarization \(\vec P\)? To make the picture and the math easier to understand, in some sense.

Imagine dielectric material in electric field. Since the material is dielectric, the external electric field induces additional electric field.

We use subscript P for quantities that is induced.

For a neutral object, one conservation law is

where \(\rho_P\) is the charge volume density distribution inside the object while \(\sigma_P\) is the surface charge density. We need the surface charge because charge will be induced on the surface.

Maxwell’s equations are very beautiful, so we definitely want to preserve Gauss’ law. Thus we define

Polarization and Electric Field

So what is the relation of \(\vec P\) and \(\vec E\)? This is not trivial. The physics here is the electric field (vector), space distribution, polarization (vector).

How many ranks do we need to describe a spatially distribution in Euler space? Two maximum. So

However, the tensor \(\chi_{ij}\) can be a function of \(E_j\) and even \(P_i\).

For **isotropic** material, we would have no direction dependence of the spatial distribution

**This is because no shear stress is possible thus :math:`P_i` doesn’t depend on E field on other directions.**

\(\chi\) can still be a function of \(E_i\) though. Now we require that this material is **linear**, which means \(\chi\) is independent of \(E_i\).

With the requirement of **isotropic and linear**, we have a simple relation,

where, for some reason, \(\chi_e\) is called **electric suscptibility**.

Why is this isotropic and linear approximation useful?

Think about the microscopic structure of material.

### Electric Potential¶

The potential should be related dipole since this is dielectric material. (<- This is a joke.) Anyway, it is related to dipole.

The electric potential is composed of two parts, one produced by the volume charge density and another produced by the surface charge density.

Recall that the electric potential by a dipole is

The electrical potential generated by induced charge distribution in dielectric material is a integral of a lot dipole electric potential fields. This can be seen by the following approximation in \(x\) direction,

then

**At this point we only considered the induced field.**

### Total Field¶

Gauss’ law tells us the field about the total charge, both the original and the induced, and the total field, whatever generated it.

where \(\rho\) means the free charge.

As we have introduced

Combine terms

The point is we don’t want to mess with \(\vec P\), so we define a **displancement vector**

so that the Gauss’ law only involve with free charge \(\rho\) and the displacement vector,

For **isotropic and linear** material, we already know \(\vec P = \chi_e \vec E\),

then we define

which is the **permittivity** by definition.

### Energy and Stress¶

Just as a reference,

which indicates that the energy density is

Stress tensor \(T_{ij}\) which appears in the formula of force \(\vec F = \int \vec{dS} \cdot \overleftrightarrow{T}\) is

## Magnetic Field¶

Lorentz force is

### Ampere‘s Law¶

Force between two charged wire is given by

Force on A Charged in Electric Field

Force on a charge in a electric field generated by \(\rho(\vec x)\) is

Charge Conservation in Magnetostatistics

Conservation of charge means

we require the source of current density 0 which means

The question is how can the electric current feels each other? By intermediate field we say.

### Biot-Savart Law¶

The magnetic field generated by electric current is given by Biot-Savart law.

Starting from magnetic field we could find the force,

Force of Charge Distribution in Electric Field

vec F = int d^3x’ frac{rho(vec x’)} vec E(vec x’) .

### Divergence and Curl of Magnetic Field¶

Using the formula of magnetic field generated by current,

The last line is true since

The curl on the other hand is not zero.

Apply a loop integral,

### Vector Potential¶

For a electric field, we have a scalar potential,

since electrostatics is curl free,

Magnetostatics, on the other hand, is divergence free. Thus we expect

Not surprisingly, we have

is also a valid vector potential for \(\vec B\).

Gauges

Coulomb gauge is

Example of Gauge Freedom

To have a magnetic field \(\vec B = B_0 \hat {e_z}\), we could use different vector potentials,

This is the gauge freedom of the magnetic field.