Find the field inside and outside a spherical geometry

Nitacii
Messages
4
Reaction score
0
Homework Statement
Given the current time dependent density ##\mathbf{J}## find the EM field inside and outside a sphere.
Relevant Equations
##\mathbf{J} = k(t) \delta(r-a) \sin \theta \hat{\boldsymbol{\phi}}##
So I tried to solve this using the Hertz potentials. I choose the magnetic one since this one corresponds to the magnetisation.
Before I start let me note that I denote a unit vector with a hat, while ##{x,y,z}## are the Cartesian coordinates and ##{r,\theta,\phi}## are the spherical coordinates.

First we need to relate the current to the magnetisation. Since for a surface current ##\mathbf{K}##, the magnetization satisfies the equation ##\mathbf{K} =\mathbf{M} \times \hat{\mathbf{n}}## we get
$$
\mathbf{M}|_{r =a} = k(t) \hat{\mathbf{z}}.
$$

The Hertz potential ##\boldsymbol{\Pi}_m## satisfies the equation

$$
\nabla^2 \boldsymbol{\Pi}_m - \frac{1}{c^2} \frac{\partial^2}{\partial t^2} \boldsymbol{\Pi}_m = -\mu_0 \mathbf{M}
$$
since the magnetiazion is non-zero only in the ##\hat{\mathbf{z}}## direction and utilizing the spherical symmetry of the problem I assumed that ##\boldsymbol{\Pi}_m = \Pi_z(r,t) \hat{\mathbf{z}}##, then the equation for the potential is

$$
\frac{1}{r} \frac{\partial^2}{\partial r^2} \left(r \Pi_z(r,t) \right) - \frac{1}{c^2} \frac{\partial^2}{\partial t^2 }\Pi_z(r,t) = -\mu_0 H(t) \delta (r-a).
$$

Now this is basically where I get stuck, I'm unable to solve this equations, I've tried combining the spherical Bessel functions so that the field is regular at the origin and zero at infinity (an outgoing wave), but I am not able to solve the equations above.

I assume when I would get the homogeneous equation for the potential I would then assume that the potential is continuous and impose the Maxwell equations by checking the boundary, note that ##\mathbf{A} = \nabla \times \boldsymbol{\Pi}## and ##\mathbf{B} = \nabla \times \nabla \times \boldsymbol{\Pi}_m##.

So the equations is basically what form of solutions satisfy the equation for the potential?

Edit: I've found an article which solves exactly the same problem avaiable on Arxiv, but they use Laplace transform and for me the derivation is quite unclear.
 
Last edited:
Physics news on Phys.org
Note that \mathbf{J} is zero everywhere except on r = a where it is not defined, but causes a jump in \mathbf{B} across r = a. Therefore you need to take the general solutions in a vacuum which are proportional to \sin \theta or \cos \theta, which are different for r < a and r > a, and adjust the arbitrary constants to satisfy the boundary conditions at r = a.
 
  • Like
Likes hutchphd and Nitacii
That's why I utilised the Hertz vectors, this automatically sets the correct dependence the on ##\theta##. If we reconstruct the vector potential ##\mathbf{A}## we get

$$
\mathbf{A} = \nabla \times \boldsymbol{\Pi} = - \sin \theta \frac{\partial}{\partial r} \Pi_z(r,t) \hat{\boldsymbol{\phi}} = A_\phi \hat{\boldsymbol{\phi}}.
$$

This seemed the correct reasoning because if we take then the equation for the vector potential ##\mathbf{A}## we get

$$
\frac{1}{r} \frac{\partial^2}{\partial r^2}\left(r A_\phi\right)+\frac{1}{r^2} \frac{\partial}{\partial \theta}\left[\frac{1}{\sin \theta} \frac{\partial}{\partial \theta}\left(\sin \theta A_\phi\right)\right]-\frac{1}{c^2} \frac{\partial^2}{\partial t^2} A_\phi =-\mu_0 J_\phi
$$
inputing the sine dependence on this we get

$$
\frac{\partial^2}{\partial r^2}(r A_\phi)-\frac{2}{r^2}(r A)-\frac{1}{c^2} \frac{\partial^2}{\partial t^2}(r A_\phi)= - \delta(r-a) \mu_0 k(t).
$$

this is the same dependence on the ##\theta## variable in the article I've just mentioned in the edit of the original question. Perhaps this is wrong...

Anyway I'm not sure what the general solution for ##\mathbf{\Pi}_m## is.

Perhpas I've didn't phrase the question clearly. I'm just asking what form of $\Pi_z$ do I have to use.
 
Last edited:
##|\Psi|^2=\frac{1}{\sqrt{\pi b^2}}\exp(\frac{-(x-x_0)^2}{b^2}).## ##\braket{x}=\frac{1}{\sqrt{\pi b^2}}\int_{-\infty}^{\infty}dx\,x\exp(-\frac{(x-x_0)^2}{b^2}).## ##y=x-x_0 \quad x=y+x_0 \quad dy=dx.## The boundaries remain infinite, I believe. ##\frac{1}{\sqrt{\pi b^2}}\int_{-\infty}^{\infty}dy(y+x_0)\exp(\frac{-y^2}{b^2}).## ##\frac{2}{\sqrt{\pi b^2}}\int_0^{\infty}dy\,y\exp(\frac{-y^2}{b^2})+\frac{2x_0}{\sqrt{\pi b^2}}\int_0^{\infty}dy\,\exp(-\frac{y^2}{b^2}).## I then resolved the two...
It's given a gas of particles all identical which has T fixed and spin S. Let's ##g(\epsilon)## the density of orbital states and ##g(\epsilon) = g_0## for ##\forall \epsilon \in [\epsilon_0, \epsilon_1]##, zero otherwise. How to compute the number of accessible quantum states of one particle? This is my attempt, and I suspect that is not good. Let S=0 and then bosons in a system. Simply, if we have the density of orbitals we have to integrate ##g(\epsilon)## and we have...
Back
Top