
#1
Mar207, 05:08 AM

P: 886

The energy of a magnetic dipole in an external magnetic field is
[tex]U =  \mathbf{m} \cdot \mathbf{B}[/tex] Yet if I try to show this by calculating the energy in the fields, I get the wrong answer. I believe the problem is that I am making some arguements which neglect strange behavior at r=0, and thus am neglecting a factor. The math is short, so let me just write it out: magnetic field of a dipole [tex]\mathbf{B}(\mathbf{r}) = \frac {\mu_0} {4\pi r^3} \left(3(\mathbf{m}\cdot\hat{\mathbf{r}})\hat{\mathbf{r}}\mathbf{m}\right) + \frac{2\mu_0}{3}\mathbf{m}\delta^3(\mathbf{r})[/tex] energy in a magnetic field [tex]U = \frac{1}{2} \int d^3 r \ \frac{1}{\mu_0}B^2[/tex] So, for a magnetic dipole in a constant external field B_0, we have: [tex]U = \frac{1}{2} \int d^3 r \ \frac{1}{\mu_0}(\mathbf{B}_0 + \mathbf{B}_{dipole})^2[/tex] [tex]U = \frac{1}{2 \mu_0} \left[ \int d^3 r \ B_0^2 + \int d^3 r \ B_{dipole}^2 + \int d^3 r \ 2 \mathbf{B}_0 \cdot \mathbf{B}_{dipole} \right][/tex] The first two terms are independent of the orientations, and thus are just constants that can be ignored when calculating the energy dependence on the dipole's orientation with respect to the external field. So continuing... [tex]U = \frac{1}{\mu_0} \int d^3 r \ \mathbf{B}_0 \cdot \left( \frac {\mu_0} {4\pi r^3} \left(3(\mathbf{m}\cdot\hat{\mathbf{r}})\hat{\mathbf{r}}\mathbf{m}\right) + \frac{2\mu_0}{3}\mathbf{m}\delta^3(\mathbf{r}) \right)[/tex] [tex]U = \frac{2}{3} \mathbf{m} \cdot \mathbf{B}_0 + 4 \pi \int d^3 r \ \frac {1} {r^3} \left( 3(\mathbf{m}\cdot\hat{\mathbf{r}})(\mathbf{B}_0 \cdot \hat{\mathbf{r}})\mathbf{B}_0 \cdot \mathbf{m} \right)[/tex] Now, chosing spherical coordinates with the "z" axis parallel to B_0 (I'll be using the physics notation where theta is the angle from z and phi is the angle in the xy plane from the xaxis), we see [tex](\mathbf{m}\cdot\hat{\mathbf{r}})(\mathbf{B}_0 \cdot \hat{\mathbf{r}}) = (m_x \cos \phi \sin \theta + m_y \sin \phi \sin \theta + m_z \cos \theta) B_0 \cos \theta [/tex] Doing the integral over [tex]d\phi[/tex] will kill the m_x and m_y terms. So we have [tex]U = \mathbf{m} \cdot \mathbf{B}_0 \left[ \frac{2}{3} + 4 \pi \int d^3 r \ \frac {3\cos^2 \theta  1} {r^3} \right][/tex] However, since [tex]\int_0^{\pi} (3\cos^2 \theta  1) \sin \theta d\theta=0[/tex] this means: [tex]U = \frac{2}{3} \mathbf{m} \cdot \mathbf{B}_0[/tex] which is wrong. My arguements for showing that each "shell" at r doesn't contribute (doing the d\theta and d\phi integral yield zero) look correct to me, except at r=0. At r=0, theta and phi aren't even valid coordinates. And since I know the answer is wrong, it seems to suggest to me that the integrand is some kind of delta function (since it seems to contribute zero for all r accept at r=0 where is must be infinite such that the integral is a constant). But how do I go about calculating that contribution at r=0? 



#2
Mar207, 01:54 PM

P: 151





#3
Mar207, 03:18 PM

P: 886

[tex]\int_0^{\pi} (3 \cos^2 \theta  1) \sin \theta d \theta [/tex] which is indeed zero. Sorry about the typo. As I said above, I believe my calculations are correct except for r=0. I don't know how to correctly handle the r=0 point mathematically, so any help would be appreciated. Thanks. 



#4
Mar507, 04:22 PM

P: 151

Energy of dipole orientation
It seems you are calculating total field energy rather than just potential energy if dipole interaction. This (seemingly) makes a difference. I have only seem it verified be cacluating Del X A = B. No promises but, maybe you gotta recheck the whole scheme.
Hope I helped More later..... 



#5
Mar607, 03:43 AM

P: 886

I'm sorry if I am misunderstanding. Can you please elaborate? The point is that the energy of the system should be the same (up to a constant) no matter what method I use to calculate it. The two methods being compared here are [tex]U = \int \mathbf{F} \cdot \mathbf{dx}[/tex] and [tex]U=\frac{1}{2}\int (\epsilon_0 E^2 + \frac{1}{\mu_0}B^2)d^3r[/tex]. All the forces / interactions between objects in this thought problem are electromagnetic, so the two answers should be the same  they are not. So there is a problem here and I need help figuring it out. Can anyone else see what is causing the mismatch? 



#6
Mar607, 03:54 PM

P: 151

Ya I got the sign wrong (there is usually a  with the electrostatic potential). I was thinking vector potential rather than potential energy. Do you know that this stuff is in Jackson's electrodynamics? The coefficient of the delta function seems inconsistent with yours, but I don't think this solves the problem with your scheme.




#7
Mar607, 04:37 PM

P: 886





#8
Mar807, 01:29 AM

P: 151

The delta function discrepancy may be because of the edition that I have (old) of the good old JD's. I checked another text and that matches yours.
OK, what I am getting is as follows: First, the potential energy measured by m \dot B is the negative of the work required to rotate the dipole in the magnetic field B: That is, the energy of configuration depending only on orientation. Moreover, the two integrals of energy you wrote depend on different limits of integration and are not easy to relate. The integral F\dot dx is a path integral and the other is a volume integral. You seem to be looking for the energy of configuration and not the work done on a test charge (as is integral of F \dot dx). And, there is no integral theorem directly used here to relate the limits (boundary) of integration. More: There is another problem. It seems that we are assuming that the uniform B field is infinite. So, an infinite volume integral of the energy density of the uniform field squared would be infinite. Indeed, the whole integral in question diverges. However, applications of the Poynting's theorem of energy seem to be used only locally (in finite volumes). If not infinite, then what volume does one chose in this configuration? The volume does effect that integral you are retaining ( the one with B_0\dot B_{dipole}). So, I am beginning to believe that there is a discrepancy in your scheme. Only trying to help: If I didn't think it interesting I would even reply. Maybe this could go to another forum as a physics Q. (I would like to know how it turns out.) 



#9
Nov909, 04:16 PM

P: 105

This is actually quite an interesting problem. Remember that in the derivation of the energy of a magnetic field we get
[tex]U = \frac{1}{2 \mu_0} \left(\int_{V} d^3 r \ B^2\int_{\partial V} (\mathbf{A} \times \mathbf{B}) \cdot d\mathbf{a}\right) [/tex] over a volume V enclosing all the charge (with boundary [tex]\partial V[/tex]). Now we can always extend this out to infinity, and if our currents are finite the last term vanishes. However if we have a constant magnetic field this is not the case, so we have to include this term. But I find that it still doesn't give me the right answer. Now our constant magnetic field B_0 can be given the vector potential [tex]\mathbf{A}_0 = \frac{1}{2} \mathbf{B}_0 \times \mathbf{r}[/tex] and the dipole has potential [tex]\mathbf{A}_{\text{dipole}} = \frac{\mu_0}{4 \pi} \frac{\mathbf{m} \times \mathbf{\hat{r}}}{r^2}[/tex] Putting these together I get [tex](\mathbf{A}_0 \times \mathbf{B}_{\text{dipole}})\cdot \mathbf{\hat{r}} = \frac{1}{2}\frac{\mu_0}{4 \pi r^2}((\mathbf{m}\cdot\mathbf{\hat{r}})(\mathbf{B}_0\cdot\mathbf{\hat{r} })  (\mathbf{B}_0 \cdot \mathbf{m}))[/tex] [tex](\mathbf{A}_{\text{dipole}} \times \mathbf{B}_0)\cdot \mathbf{\hat{r}}=2 (\mathbf{A}_0 \times \mathbf{B}_{\text{dipole}})\cdot \mathbf{\hat{r}}.[/tex] Now noticing the integral over a spherical shell: [tex]\int d \Omega (\mathbf{m}\cdot\mathbf{\hat{r}})(\mathbf{B}_0\cdot\mathbf{\hat{r}}) = \frac{4 \pi}{3} (\mathbf{B}_0 \cdot \mathbf{m}) [/tex] it follows directly that the integral of the surface component over a spherical shell radius R is [tex]\int (\mathbf{A}_0 \times \mathbf{B}_{\text{dipole}} + \mathbf{A}_0 \times \mathbf{B}_{\text{dipole}}) = \frac{1}{2}\frac{\mu_0}{4 \pi R^2} R^2(4 \pi  \frac{4 \pi}{3}) (\mathbf{B}_0 \cdot \mathbf{m}) [/tex] [tex]= \frac{\mu_0}{3} (\mathbf{B}_0 \cdot \mathbf{m})[/tex] And here is my problem: this term doesn't correct it; it gives [tex]U=\frac{2}{3} (\mathbf{B}_0 \cdot \mathbf{m})  \frac{1}{6} (\mathbf{B}_0 \cdot \mathbf{m})[/tex] [tex] U = \frac{1}{2} (\mathbf{B}_0 \cdot \mathbf{m})[/tex] So either I dropped some factors in my calculation or there is something else wrong with using this formula in this situation. 


Register to reply 
Related Discussions  
"paradox" regarding energy of dipole orientation  Classical Physics  80  
Derivation of dipoledipole interaction energy  Advanced Physics Homework  2  
Dipole Dipole and Lewis Dot structes  Chemistry  1  
deriving the energy of a magnetic dipole in a magnetic field  Introductory Physics Homework  8  
deriving the energy of a magnetic dipole in a magnetic field  Classical Physics  1 