arXiv:physics/0611099v1 [physics.class-ph] 10 Nov 2006 a,wihi sal nrdcdjs sa nlgo the rule of empirical analog an be- an as as or forces Coulomb law just Coulomb magnetic the introduced electric the usually and is with which monopoles calculated consid- law, of are are ensemble poles poles tween permanent an The between as objects. ered forces magnetized the or magnets describe very are to Maxwell poles the convenient magnetic in however, Practically, or theory equations. electromagnetic of framework vie sfra possible as far as avoided of Again because electromagnet. reasonable equivalent roles. more the special with much no similarity is play the model magnet bar current poles, the of the model, ends this equiv- the In is or that magnet. electromagnet permanent an the distri- have to current we alent in coil, macroscopic shown a the as by mimicking bution surface By the along 1(c). Fig. bar which the current around macroscopic have circulates we currents, loop scopic 1(a). Fig. average. in spatial shown the canceled as through are end inside out end charges other one magnetic the at the of at pole Contributions pole north south magnetic a a have and we in magnet, dipoles bar magnetic a microscopic of ensemble oriented the ossec,w aeto have we consistency, hc r bandb oregann fmicroscopic of distributions coarse-graining source by and fields obtained are which space and time the provides also to It respect inversions. electrons. with of symmetries such spins correct rotations or mo- kinetic motions magnetic orbital with as microscopic associated lat- always and magnetic are The no far ments are so there field. found because natural magnetic monopoles more same is loop-current model the the ter exactly and In model give dipole loop. the model consider current fields, to far circulating of is a terms modeling as moment of magnetic way magnetic each the Another of to pair analogous dipole. is displaced It electric slightly polarities. a opposite modeled as with be monopoles i.e., can dipole, moment a magnetic magne- as A the along direction. oriented tization are which moments magnetic h s ftento fmgei oesol be should pole magnetic of notion the of use The micro- the coarse-grain we when hand, other the On omlyw r neetdi arsoi quantities, macroscopic in interested are we Normally microscopic of ensemble an is magnet permanent A eiaino antcCuobslwfrti,semi-infinit thin, for law Coulomb’s magnetic of Derivation edrv h antcCuoblwfo h ai principles force. basic Lorentz the the and from place equations law Maxwell monopoles Coulomb the magnetic magnetic for the that derive We to corresponds which law, ti hw httemgei oc ewe hn semi-infini thin, between force magnetic the that shown is It .INTRODUCTION I. derive RS,JpnSineadTcnlg gny oy 103-0028 Tokyo Agency, Technology and Science Japan CREST, 2 eas fisasnei the in absence its of because h antcCuoblaw Coulomb magnetic the 1 hnw coarse-grain we When . yt nvriy asr,Koo6581,Jpnand Japan 615-8510, Kyoto Katsura, University, Kyoto eateto lcrncSineadEngineering, and Science Electronic of Department 5 o logical For . Dtd eray2 2008) 2, February (Dated: aa Kitano Masao where h oc cigo charge a on acting force the ucin ihtee h clrptnilfrapoint a for potential scalar the these, charge With function. stoopst antcplsa ahed() fw adopt wal side we the If around current macroscopic circulating (a). the of end moment, consists magnetic picture each for at model loop poles pictur current the magnetic macroscopic opposite corresponding two the microscopic momentis dipole, of magnetic magnetic each ensemble a consider an we as is If magnet (b). permanent moments A magnetic 1: FIG. ent that note We be will works, previous the in roles rigorously. mentioned treated crucial not plays was which singularity currents, but field loop The the infinitesimal on provided. for directly be based will analysis laws detailed fundamental more a paper this oeta ot er g yChen by ago years forty than more monopoles. magnetic of of none notion force, the Lorentz contain the which and equations Maxwell’s from q aadlwcnb xrse sd as expressed be can law Savard etrfunction vector 2 o rvt,w nrdc clrfunction scalar a introduce we brevity, For h eiaino antcCuobslwwsgiven was law Coulomb’s magnetic of derivation The at a b (c) (b) (a) ∗ I URN EST O THIN A FOR DENSITY CURRENT II. r δ q 2 3 ( lcda h rgnis origin the at placed r tteedpit fec solenoid. each of points end the at d antcdipole magnetic is = ) G eslnisoesaCoulomb-type a obeys solenoids te feetoants,nml from namely electromagnetism, of F 0 ( 1 δ r ∇ ← ( G = ) x 2 G 1 ) 0 δ ( = fposition of ( 4 = y π SOLENOID ) 1 | q δ r − 1 ( | z q Japan , G , 2 stetredmninldelta dimensional three the is ) /ε 1 q G 1 0 and ) solenoids e at 1 r G φ ( H r ( = ( 1 r r = ) ∇ ( urn loop current 1 r ( = ) d = ,y z y, x, rmaohrcharge another from 1 3 · − G 4 n Nadeau and C π q/ε 1 r r | × r 2 ): = | .TeBiot- The ). G 3 0 ) . 1 δ G ( 3 G r 0 ( ( r ,where ), 0 r hold, ) ,and ), n a and 4 In . (c). l (1) e 2
dla r2
La
r1 κadla loops
r3 m
r4 Lb
FIG. 2: A pair of thin, semi-infinite solenoids, La and Lb
′ dC = Idl is a current moment (current I times length where dv′ is a volume element at r′ and ∇ = ∂/∂r′. dl) located at the origin. Using Eq. (4), it can be rewritten as As shown in Fig. 2, a thin solenoid can be constructed 3 as a stack of tiny loop currents at least in principle6. The Hm(r)= −(m · ∇)G1(r)+ mδ (r), (6) current density for an infinitesimal loop current place at ∇ 3 ∇ the origin is where the relations · G1(r) = δ (r) and × G1 = 0 have been used. As shown in Appendix, Eq. (6) can also 3 J m(r) = (−m × ∇)δ (r) (2) be derived from the Maxwell equations. It is well known that the field for an electric dipole where m is the magnetic moment for the loop current and p D r p ∇ G r 2 moment is ( )= −( · ) 1( ), which contains no its unit is (A m ) (See Appendix). To form a solenoid we delta-function terms unlike Eq. (6). This means that in stack these tiny loop currents along a curve La with a terms of near fields, the dipole and the current loop do constant line density κa (loops per unit length). Each not yield the same field. loop are aligned so that the direction of tangent of curve From Eq. (6), the strength of the magnetic field created and that of m coincide. by a line segment dla is The current density distribution for a line segment 3 dla ( k dm ) is dH(r) = (Cadla · ∇a)G1(r − ra)+ Cadlaδ (r − ra) 3 (7) dJ(r) = (−Cadla × ∇)δ (r − ra), (3)
D where ∇a = ∂/∂ra. The integration along a curve La where Ca = κam ∼ A m is the magnetic moment per unit yields the magnetic field created by the solenoid: length and characterizes the strength of the solenoid7, where m = m · dla/|dla| is the magnitude of m and ra H(r)= dH = Ca [G1(r − r1) − G1(r − r2)] represents the position of the segment dla. ZLa 3 + Ca δ (r − ra)dla, (8) III. FIELD BY A THIN, SEMI-INFINITE ZLa SOLENOID where r2 and r1 are the start and the end points of L . The second term of the right-hand side, which cor- Here we introduce an important formula. For a con- a responds to the magnetic flux confined in the solenoid, stant vector α and a vector field V (r), we have vanishes outside. For semi-infinite cases (r2 = ∞), we (α × ∇) × V − α × (∇ × V ) have ∇ ∇ = (α · )V − α( · V ). (4) 3 H(r)= CaG1(r − r1)+ Ca δ (r − ra)dla. (9) ZL With the Biot-Savard law and Eq. (2), we can find a the strength of magnetic field generated by a magnetic The first term of the right-hand side, which represents moment m placed at the origin as the magnetic field outside of the solenoid, is equivalent to the field for a monopole ga = µ0Ca located at r1: ′ ′ ′ Hm(r)= dv J m(r ) × G1(r − r ) Z g r − r B r G r r a 1 ∇ ( )= ga 1( − 1)= 3 . (10) = (−m × ) × G1(r), (5) 4π |r − r1| 3
The dimension of ga, V. COULOMB’S LAW BETWEEN TWO THIN, SEMI-INFINITE SOLENOIDS D H Vs/A ga = µ0Ca ∼ Am= Am=Vs=Wb (11) m m Now we can calculate the force (16) on solenoid b by the field (9) generated by solenoid a; correctly corresponds to that for the magnetic charge.
The magnetic flux ga confined along the solenoid fans out F b←a = CbB(r3)= Cbµ0H(r3) isotropically from the end point r . As seen in Fig. 2, a 1 = µ C C G (r − r ). (17) thin, semi-infinite solenoid can be viewed as a magnetic 0 b a 1 3 1 monopole located at the end. This equation corresponds to Coulomb’s law for magnetic charges ga at r1 and gb at r3:
gagb gagb r3 − r1 IV. MAGNETIC FORCE ACTING ON A F ← G r r 3 1 = 1( 3 − 1)= 3 . (18) SEMI-INFINITE SOLENOID µ0 4πµ0 |r3 − r1| These two equations are exactly the same but the for- From Eq. (2), we see that the Lorentz force acting on mer is for the total (integrated) force between currents a tiny loop current m placed at r in a magnetic field B of solenoids and the latter is for the force between mag- is netic point charges. It is interesting that the force is independent of the ′ ′ ′ F m = dv J m(r ) × B(r + r ) paths of either solenoids as long as their end points and Z strengths are kept constant. We should note that the ∇ = (m × ) × B(r). (12) derivation of the magnetic Coulomb law is not straight- forward. The properties of static, source-free fields must ∇ Using Eq. (4) and the conditions for magnetic field: · be utilized. ∇ −1 B = 0 (divergence-free) and × (µ0 B) = 0 (rotation- free), the expression can be modified as VI. TORQUE ON A SEMI-INFINITE SOLENOID F m = (m · ∇)B(r), (13) which is suitable for line-integral. Equation (12) can also Now we are interested in the torque due to the inte- grated magnetic force on a thin, semi-infinite solenoid. be modified as F m = ∇(m · B) with the divergence-free The torque N m with respect to the origin exerted on a condition only. The rotation-free condition is satisfied 4 only when the right hand side of the Maxwell-Amp`ere loop current m located at r is ∇ equation, × H = J + ∂D/∂t, vanishes. ′ ′ ′ ′ N m = dv (r + r ) × [J m(r ) × B(r + r )] Thus the magnetic force acting on a line element dlb Z at rb is ′ ′ ′ ′ = r × F m + dv [B(r + r ) × J m(r )] × r Z dF = (Cbdlb · ∇b)B(rb), (14) (0) (1) = N m + N m , (19) where Cb = κbm and ∇b = ∂/∂rb. Integration along a curve Lb yields the total force acting on the solenoid; where the first term corresponds to the torque due to the total force F m with respect to the origin and the second term corresponds to the torque with respect to the center F = dF = Cb (dlb · ∇b)B ZL ZL of the loop current. Using Eq. (13), the former can be b b modified as = Cb [B(r3) − B(r4)] , (15) (0) N m = r × F m = r × (m · ∇)B where r and r are the initial and end points of L , 4 3 b = (m · ∇)(r × B) − [(m · ∇)r] × B respectively. For semi-infinite cases [B(r = ∞) = 0], 4 ∇ we have = (m · )(r × B) − m × B, (20) with (m · ∇)r = m. The latter can be simplified as F = CbB(r3)= gbH(r3), (16) (1) ′ ′ ′ 3 ′ ′ −1 N m = dv B(r + r ) × (−m × ∇ )δ (r ) × r with H = µ B. This is equal to the force for a point Z 0 magnetic charge gb = µ0Cb placed at r3. ′ ′ ′ = r × [(m × ∇ ) × B(r + r )] ′ It is surprising that the sum of the forces acting on each r =0 B r r′ m ∇ ′ r′ part of the solenoid can be represented in terms only of + [ ( + ) × ( × )] × r′=0 ′ ′ ′ the magnetic field B(r3) at the end point r3. This is B ∇ m B m ∇ r = 0 + [( · ) − ( · ) ] × r′ =0 because the partial force is proportional to the (vector) gradient of the magnetic field. = m × B, (21) 4
′ ′ where we have utilized ∇ × r′ = 0 with ∇ = ∂/∂r′. . by a pair of small vectors, a and b. The center is located Finally, we have at the origin and current I is circulating along the edge. The current distribution on the four segments can be (0) (1) ∇ N m = N m + N m = (m · )[r × B(r)]. (22) approximately represented as The torque on a line element dl at r is b b 3 3 J I(a×b)(r)=I aδ (−b/2)+ bδ (a/2) dN = Cb(dlb · ∇b)[rb × B(rb)], (23) − aδ3(b/2) − bδ3(−a/2) (A.1) and the integration along the semi-infinite solenoid Lb terminated at r3 yields In the limit of |a|, |b| → 0 with m = I(a × b) being kept constant, it approaches N = dN = Cb (dlb · ∇b)[rb × B(rb)] ZL ZL 3 3 b b J m(r)= I −a(b · ∇)δ (r)+ b(a · ∇)δ (r) = C r × B(r )= r × F , (24) b 3 3 3 = [−I(a × b) × ∇]δ3(r) where F is the force (16) on the solenoid. Surprisingly, = (−m × ∇)δ3(r), (A.2) the torque exerted on a thin, semi-infinite solenoid coin- cides with that for a point magnetic charge gb = µ0Cb and Eq. (2) is obtained. The derivative of the delta func- place at r3. tion ∇δ3(r) serves as a differential operator when it is integrated together with other functions as in 1D cases: dxf(x)(d/dx)δ(x) = −(df/dx)(0) . For example, the VII. DISCUSSION magneticR force on m can be calculated as follows;
In this paper we only deals with thin solenoids. 3 F m = dv(−m × ∇)δ (r) × B(r) Solenoids with finite cross-section can be represented as a Z bundle of thin solenoids. Similarly a permanent magnet = dv −∇δ3(m · B)+ m(∇δ3 · B) can naturally be modeled as a bundle of thin solenoids. Z The degree of magnetization of magnets or magnetized ∇ ∇ objects is characterized by the quantity called macro- = [ (m · B) − m( · B)]r=0 D scopic magnetization, M ∼ A/m, which is defined as = [(m × ∇) × B] (0). (A.3) D volume density of magnetic moments, m ∼ A m2. The magnetization can also be represented as area-density of Next we derive Eq. (6) from Maxwell’s equations. The ∇ H thin solenoids. We remember that a thin solenoid is char- Amp`ere law for an infinitesimal loop current, × = D J m acterized by the current moment, C ∼ A m. Bundling ν can be written as D solenoids per unit area (∼ /m2), we can create the mag- ∇ × H(r) − mδ3(r) =0. (A.4) netization M = νC. In conclusion, we have shown that the force between We know that a rotation-free field can be represented as thin, semi-infinite solenoids obeys the Coulomb law that a gradient of some scalar field φ (r) as is for the equivalent magnetic point charges placed at m each end of solenoid. We have also shown that the torque H(r) − mδ3(r)= −∇φ (r). (A.5) exerted on the solenoid in a magnetic field coincides with m that for a corresponding magnetic point charge. It is Taking the divergence of each side, we have convenient to introduce magnetic charges or poles be- cause the magnetic Coulomb law can easily be applied 3 2 (∇ · m)δ (r)= ∇ φm(r), (A.6) for forces between them. But we should remember that without the justification given in this paper it is just a where ∇ · (µ0H) = 0 has been used. Comparing it with rough and ready method. 2 3 ∇ (∇ · m)G0(r)= −(∇ · m)δ (r), (A.7)
Acknowledgments 2 3 which is obtained from ∇ G0(r) = −δ (r), we find the solution of (A.6) to be This work is supported by the 21st Century COE program No. 14213201. φm(r)= −(∇ · m)G0(r). (A.8)
APPENDIX: CURRENT DENSITY AND FIELD With Eq. (A.5) and ∇G0 = −G1 we have Eq. (6): OF A TINY CURRENT LOOP 3 H(r)= −(m · ∇)G1(r)+ mδ (r), (A.9) The current density (2) for a tiny current loop can be derived as follows. We consider a parallelogram defined which can be confirmed to satisfy the Maxwell equations. 5
∗ Electronic address: [email protected] (1952). 1 J. D. Jackson, Classical Electrodynamics (Addison Wesley 6 In practice, a solenoid is made as a helix of a conducting and Sons, New York, 1998), 3rd ed. wire not as a stack of closed wire loops. The difference in 2 F. W. Warburton, “The magnetic pole, A useless concept,” current distribution, which can be represented as an addi- Am. Phys. Teacher 2, 1 (1934). tional current along the solenoid, can be made arbitrarily 3 H. S. C. Chen, “Note on the magnetic pole,” Am. J. Phys. small by reducing the wire current and increasing the num- 33, 563 (1965). ber of winding correspondingly. 4 D G. Nadeau: “Comment on Chen’s note on the magnetic 7 The relation A ∼ B means A and B are dimensionally pole,” Am. J. Phys. 34, 60 (1966). equivalent and is normally written as [A] = [B]. 5 J. Goldemberg, “An experimental verification of the Coulomb law for magnetic poles,” Am. J. Phys. 20, 591–592