We shall assume that the magnetic field $\overrightarrow{B}$ and the angular velocity vector $\overrightarrow{\omega}$ point in the same direction ($\overrightarrow{\omega}\cdot \overrightarrow{B}>0$). The Lorentz force causes the positive charges to move outward in the sphere, so a negative space charge will build up in the interior, while the surface charge will be positive.

The charges distribute themselves in such a way that the electric field balances the Lorentz force. The sphere is not a perfect conductor, so in a steady state the charges must move with exactly the same velocity $\overrightarrow{v}=\overrightarrow{\omega}\times \overrightarrow{r}$ as the sphere at every point. (To keep the charges in motion relative to the sphere an electromotive force would be necessary.)

$$\overrightarrow{E}=-\overrightarrow{v}\times \overrightarrow{B}$$ | (1) |

The space charge density $\rho $ can be found by taking the divergence of this equation.

Here we used the fact that a mixed product does not change when its terms are cyclically permuted: $\overrightarrow{a}\cdot (\overrightarrow{b}\times \overrightarrow{c})=\overrightarrow{b}\cdot (\overrightarrow{c}\times \overrightarrow{a})$. According to Maxwell’s equations $\nabla \times \overrightarrow{B}={\mu}_{0}\overrightarrow{j}={\mu}_{0}\rho \overrightarrow{v}$.

$$\nabla \times \overrightarrow{v}=\nabla \times (\overrightarrow{\omega}\times \overrightarrow{r})=\overrightarrow{\omega}\nabla \cdot \overrightarrow{r}-(\overrightarrow{\omega}\cdot \nabla )\overrightarrow{r}=2\overrightarrow{\omega}$$

Substituting these results back into equation (2), $\rho ={\epsilon}_{0}{\mu}_{0}{v}^{2}\rho -2{\epsilon}_{0}\overrightarrow{\omega}\cdot \overrightarrow{B}$. Using the formula ${\epsilon}_{0}{\mu}_{0}=1\u2215{c}^{2}$ we get

$$\rho (1-{v}^{2}\u2215{c}^{2})=-2{\epsilon}_{0}\overrightarrow{\omega}\cdot \overrightarrow{B}.$$ | (3) |

Since $v\ll c$, $1-{v}^{2}\u2215{c}^{2}\approx 1$. To use equation (3) to get the space charge density, the magnetic field $\overrightarrow{B}$ needs to be known. But the homogeneous magnetic field will be distorted by the field originating from the circulating charges in the sphere. If the sphere is not spinning very fast, this correction is relatively small and the following method can be used to approximate $\overrightarrow{B}$:

First we assume that the applied homogeneous magnetic field is unchanged and denote it by ${\overrightarrow{B}}_{0}$, then calculate the space charge $\rho $ and current densities $\overrightarrow{j}=\rho \overrightarrow{v}$. Next, we use these currents to find the correction ${\overrightarrow{B}}_{1}$ to the homogeneous field, and find the corrected charge densities and currents using ${\overrightarrow{B}}_{0}+{\overrightarrow{B}}_{1}$. With these new current densities the next correction ${\overrightarrow{B}}_{2}$ can be calculated. The result can be made more accurate by doing more iterations.

Without carrying out the actual calculations, we know that ${B}_{1}~{\mu}_{0}I\u2215R$, where $I$ is a current-like quantity: $I~\omega R\rho {R}^{2}$. But $\rho ~{\epsilon}_{0}\omega {B}_{0}$, so ${B}_{1}~{\mu}_{0}{\epsilon}_{0}{\omega}^{2}{R}^{2}{B}_{0}={v}^{2}\u2215{c}^{2}{B}_{0}$. Thus, when using the $v\ll c$ approximation, the corrections to the homogeneous field due to the induced rotating space charge can be neglected. Further on, $\overrightarrow{B}$ shall be considered homogeneous.

Now we proceed with calculating the electric field and potential inside the sphere. Substituting $\overrightarrow{v}=\overrightarrow{\omega}\times \overrightarrow{r}$ into equation (1) we get $\overrightarrow{E}=-\overrightarrow{r}(\overrightarrow{\omega}\cdot \overrightarrow{B})+\overrightarrow{\omega}(\overrightarrow{B}\cdot \overrightarrow{r})$. $\overrightarrow{\omega}$ and $\overrightarrow{B}$ are parallel, so $\overrightarrow{E}=-\overrightarrow{a}(\overrightarrow{\omega}\cdot \overrightarrow{B})$, where $a$ is the distance from the axis of rotation: if we choose $\widehat{z}\parallel \overrightarrow{\omega}$, then $\overrightarrow{a}=x\widehat{x}+y\u0177$.

It must be noted here that the electric field can compensate the Lorentz force only if $\nabla \times (\overrightarrow{v}\times \overrightarrow{B})=0$. This is true for a homogeneous $\overrightarrow{B}$-field approximation when the axis of rotation is parallel to the magnetic field, but not otherwise. If the Lorentz force can not be completely cancelled, there will always be some charge motion relative to the sphere until the energy dissipation causes the spinning sphere to stop.

The electric potential inside the sphere is

$$\varphi =\frac{{a}^{2}}{2}(\overrightarrow{\omega}\cdot \overrightarrow{B})+{\varphi}_{0}$$

In spherical coordinates $\varphi ={\varphi}_{0}-{r}^{2}{sin}^{2}\theta \phantom{\rule{3.26212pt}{0ex}}(\overrightarrow{\omega}\cdot \overrightarrow{B})\u22152$, where $\theta $ is the angle between $\overrightarrow{r}$ and the axis of rotation. Unlike $\overrightarrow{E}$, the potential is continuous over the surface of the sphere, so the boundary condition $\varphi ={\varphi}_{0}-{R}^{2}{sin}^{2}\theta \phantom{\rule{3.26212pt}{0ex}}(\overrightarrow{\omega}\cdot \overrightarrow{B})\u22152$ may be used to also find $\varphi $ outside the sphere of radius $R$.

Let us approximate a potential that satisfies this boundary condition with multipole terms! The monopole term must be zero because the total charge of the sphere is zero. The dipole term is also zero because the dipole field does not match the symmetry of this system. Therefore the dominant part of the sought potential must be a quadrupole potential. It is easy to see that the potential ${\varphi}^{\left(2\right)}$ of a quadrupole that was created by superposing two dipoles along their axes (to satisfy the requirement of azimuthal symmetry) satisfies the boundary conditions exactly:

$${\varphi}^{\left(2\right)}=\left(\text{someconstant}\right)\frac{3{sin}^{2}\theta -2}{{r}^{3}}.$$

We need to seek no further: outside the sphere the potential will be

$$\varphi =\frac{{R}^{5}(\overrightarrow{\omega}\cdot \overrightarrow{B})}{2{r}^{3}}\left({sin}^{2}\theta -\frac{2}{3}\right).$$

The surface charge density is $\sigma ={\left.{\epsilon}_{0}({\partial}_{r}{\varphi}_{outside}-{\partial}_{r}{\varphi}_{inside})\right|}_{r=R}$.

$$\sigma ={\epsilon}_{0}R(\overrightarrow{\omega}\cdot \overrightarrow{B})\left(\frac{5}{2}{sin}^{2}\theta -1\right)$$

By integrating this over the surface of the sphere we see that indeed the total surface charge equals the total volume charge in magnitude, so the sphere has no net charge.

The potential inside the sphere is

$$\varphi =\frac{{a}^{2}}{2}(\overrightarrow{\omega}\cdot \overrightarrow{B})+{\varphi}_{0}=\frac{{r}^{2}{sin}^{2}\theta}{2}(\overrightarrow{\omega}\cdot \overrightarrow{B})+{\varphi}_{0}.$$

The electric field inside the sphere is

$$\begin{array}{llll}\hfill \overrightarrow{E}& =-\overrightarrow{a}(\overrightarrow{\omega}\cdot \overrightarrow{B})\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill {E}_{r}& =-r{sin}^{2}\theta \phantom{\rule{3.26212pt}{0ex}}(\overrightarrow{\omega}\cdot \overrightarrow{B})\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\\ \hfill {E}_{\theta}& =-rsin\theta cos\theta \phantom{\rule{3.26212pt}{0ex}}(\overrightarrow{\omega}\cdot \overrightarrow{B}).\phantom{\rule{2em}{0ex}}& \hfill & \phantom{\rule{2em}{0ex}}\end{array}$$The potential outside the sphere is

$$\varphi =\frac{{R}^{5}(\overrightarrow{\omega}\cdot \overrightarrow{B})}{2{r}^{3}}\left({sin}^{2}\theta -\frac{2}{3}\right).$$

The electric field outside the sphere is

$${E}_{r}=-\frac{3{R}^{5}(\overrightarrow{\omega}\cdot \overrightarrow{B})}{2{r}^{4}}\left({sin}^{2}\theta -\frac{2}{3}\right)=\frac{{R}^{5}(\overrightarrow{\omega}\cdot \overrightarrow{B})}{4{r}^{4}}(1+3cos2\theta )$$

$${E}_{\theta}=\frac{{R}^{5}(\overrightarrow{\omega}\cdot \overrightarrow{B})}{{r}^{4}}cos\theta sin\theta .$$

The surface charge density is

$$\sigma ={\epsilon}_{0}R(\overrightarrow{\omega}\cdot \overrightarrow{B})\left(\frac{5}{2}{sin}^{2}\theta -1\right).$$