Arxiv:Physics/0611099V1
Total Page:16
File Type:pdf, Size:1020Kb
Derivation of magnetic Coulomb’s law for thin, semi-infinite solenoids Masao Kitano∗ Department of Electronic Science and Engineering, Kyoto University, Katsura, Kyoto 615-8510, Japan and CREST, Japan Science and Technology Agency, Tokyo 103-0028, Japan (Dated: February 2, 2008) It is shown that the magnetic force between thin, semi-infinite solenoids obeys a Coulomb-type law, which corresponds to that for magnetic monopoles placed at the end points of each solenoid. We derive the magnetic Coulomb law from the basic principles of electromagnetism, namely from the Maxwell equations and the Lorentz force. I. INTRODUCTION from Maxwell’s equations and the Lorentz force, none of which contain the notion of magnetic monopoles. The derivation of magnetic Coulomb’s law was given A permanent magnet is an ensemble of microscopic 3 4 magnetic moments which are oriented along the magne- more than forty years ago by Chen and Nadeau . In tization direction. A magnetic moment can be modeled this paper a more detailed analysis based directly on the as a dipole, i.e., as a slightly displaced pair of magnetic fundamental laws will be provided. The field singularity monopoles with opposite polarities. It is analogous to the for infinitesimal loop currents, which plays crucial roles electric dipole. Another way of modeling is to consider but was not mentioned in the previous works, will be each magnetic moment as a circulating current loop. In treated rigorously. terms of far fields, the dipole model and the loop-current model give exactly the same magnetic field. The lat- ter model is more natural because there are no magnetic II. CURRENT DENSITY FOR A THIN monopoles found so far and microscopic magnetic mo- SOLENOID ments are always associated with kinetic rotations such as orbital motions or spins of electrons. It also provides For brevity, we introduce a scalar function G0 and a correct symmetries with respect to the time and space vector function G1 of position r = (x,y,z): inversions. 1 r Normally we are interested in macroscopic quantities, G0(r)= , G1(r)= . (1) which are obtained by coarse-graining of microscopic 4π|r| 4π|r|3 fields and source distributions1. When we coarse-grain ∇ ∇ 3 the oriented ensemble of microscopic magnetic dipoles in We note that G0 = −G1 and · G1 = δ (r) hold, 3 a bar magnet, we have a magnetic north pole at one end where δ (r)= δ(x)δ(y)δ(z) is the three dimensional delta and a south pole at the other end as shown in Fig. 1(a). function. With these, the scalar potential for a point Contributions of the magnetic charges inside are canceled charge q placed at the origin is φ(r) = (q/ε0)G0(r), and out through the spatial average. the force acting on a charge q1 at r1 from another charge r F ← G r r On the other hand, when we coarse-grain the micro- q2 at 2 is 1 2 = (q1q2/ε0) 1( 1 − 2). The Biot- H C G r scopic loop currents, we have macroscopic current which Savard law can be expressed as d = d × 1( ), where circulates around the bar along the surface as shown in Fig. 1(c). By mimicking the macroscopic current distri- bution by a coil, we have an electromagnet that is equiv- magnetic dipole current loop alent to the permanent magnet. In this model, the poles, or the ends of bar magnet play no special roles. Again the current model is much more reasonable because of the similarity with the equivalent electromagnet. The use of the notion of magnetic pole should be arXiv:physics/0611099v1 [physics.class-ph] 10 Nov 2006 avoided as far as possible2 because of its absence in the framework of electromagnetic theory or in the Maxwell equations. Practically, however, magnetic poles are very (a) (b) (c) convenient to describe the forces between permanent magnets or magnetized objects. The poles are consid- FIG. 1: A permanent magnet is an ensemble of microscopic ered as an ensemble of monopoles and the forces be- magnetic moments (b). If we consider each magnetic moment tween poles are calculated with the magnetic Coulomb as a magnetic dipole, the corresponding macroscopic picture law, which is usually introduced just as an analog of the is two opposite magnetic poles at each end (a). If we adopt electric Coulomb law or as an empirical rule5. For logical the current loop model for magnetic moment, the macroscopic consistency, we have to derive the magnetic Coulomb law picture consists of circulating current around the side wall (c). 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.