We are given a task to characterize the solution to maxwell's equations that enables a non radiative source term.

We will base our analysis of subset of solution to maxwells equations in free space of the following form and for that we will use the linear mapping:

Worth noting is that we have the plane wave expanded in spherical harmoics as:

This and the orthogonality of the spherical hormonics means that

Also because every multiplicatoin of spherical harmonics can be represented as

we can integrate all those combinations, more specifilly if \(l'' = l_0\) then we have

More specifically if \(F(\hat n)\) can be decomposed in a set of linear combinations of spherical harmoics \(Y_{l',m'}\) so that for all those we have the same \(l_0\) then

With some analysis one can show that \(\hat \varphi_{\hat n}\) and \(\hat \theta_{\hat n}sin(\theta)\) has coordinates that can be decomposed in spherical harmonics with \(l' = 1\) and that if we take \(l=0\) we obviously get

and

Set

and

and

e.g. a force free field.

we now have

and

Now

So,

A final rule to note is that

if \(F\) is independent of \(\hat r\),

with this we see that

and

and

and

Let's define the scalar potential \(\phi\), with \(\rhoi\) constant as

and the vector potential, with \(a,b\) constant as

Now assume the Lorenz gauge and the we must fullfill the continuity equation:

Which translates to

\(f_1\) and \(f_2\) is orthogonal to \(\hat n\) and \(\hat n \cdot \hat n = 1\) so we have:

Set \(d -> \rho d\) and we have

Furthermore we have that \(\phi\) shall satisfy the wave equation with the speed \(c\) so that means:

and hence

With this we have a solution to maxwells equations in free space. The next step is to introduce source terms at a spherical shell at a certain distance. This is a bit tricky. First we must patch solutions together so that the postential is continuous as we move radially. How to do this? well we note from the discusion about \(f_1,f_2\) that

and

and hence let assume at a certain distance \(r_0\),

This means that according to theory, not only do we know that Maxwell's equations will not radiate, we can take the scalar potentials equal to zero and patch with a new solution for \(r > r_0\) where wi replace \(j_0(kr)\) with \(y_0(kr)\) and assume \(j_1'(k r_0)=h y_1'(k r_0)\) which imply

Because \(j_i,y_i\) both follows the same recurrence equation we coclude that the resulting solution indeed solves Maxwells equations e.g we use

and

We have \(E = - \nabla \phi - \frac{\partial}{\partial t}A\) that is

Now for a change in the interface for \(E\) over a phase we know that A is continuous and the difference comes from \(\phi\). Therefore we calculate

The derivative with respect to the spherical harmonics can be neglected because that the bessel function is untouched and hence at \(r_0\) zero. So:

So over the interface 1.

Turning the focus to the magnetic part we have \(B = \nabla \times A\) which we can write

Turning over to the boundary condition we see that if we take

where \([F]_{i,j}'\) is \([F]_{i,j}\) with \(j_i\) replaced by \(y_i\). e.g. we assume that the maxwell solution

and

for \(r > r_0\). $$

because the radial part has \(j_1(kr_0) = 0\) as a factor.

Furthermore:

As said the radial part of \(B1\) dissapears and left is a form of type

But

The divergence of this is

Furthermore for the source terms we have

for a flow on the sphere. Now we can match \(b\) so that the continuity equation for the source terms are satisfied and we are done. note that a is arbritary in this evaluation and can be taken 0 or we could have design it with the help of complex numbers so that both \(B\) satisfy \(\nabla \times B = k B\) and \(\nabla \times B' = k B'\) e.g. force free.

## Comments

comments powered by Disqus