arxiv:physics/ v2 [physics.class-ph] 20 Jan 2004

Transcription

1 Potentials of a uniformly moving point charge in the Coulomb gauge arxiv:physics/374v [physics.class-ph] Jan 4 V Hnizdo National Institute for Occupational Safety and Health, Morgantown, West Virginia 655, USA vbh5@cdc.gov Abstract The Coulomb-gauge vector potential of a uniformly moving point charge is obtained by calculating the gauge function for the transformation between the Lorenz and Coulomb gauges. The expression obtained for the difference between the vector potentials in the two gauges is shown to satisfy a Poisson euation to which the inhomogeneous wave euation for this uantity can be reduced. The right-hand side of the Poisson euation involves an important but easily overlooked delta-function term that arises from a second-order partial derivative of the Coulomb potential of a point charge.. Introduction Gauge invariance is an important property of electrodynamics. Notwithstanding the thorough attention it has received in textbooks [ 4], the topic of gauge invariance seems to be still very much alive. Recently, articles have appeared on the resolution of apparent causality problems in the Coulomb gauge [5], the transformation from the Lorenz gauge to the Coulomb and some other gauges [6], the Coulomb-gauge vector potential in terms of the magnetic field [7], and the historical development of the whole concept of gauge invariance [8]. A discordant voice in this is a claim by Onoochin [9] that the electric field of a uniformly moving point charge comes out differently when it is calculated in the Lorenz and Coulomb gauges, which that author takes as evidence that the two gauges are not physically euivalent. Curiously, it is difficult to find in the literature an explicit expression for the Coulombgauge vector potential of a uniformly moving charge. To the present author s knowledge, such an expression has appeared only in a short comment [], where it is obtained directly from the well-known electric field and Coulomb-gauge scalar potential of the charge. The formulae of Jackson [6] for the transformation between the Lorenz and Coulomb gauges thus have come timely to provide an analytical check on the claim of Onoochin. The purpose of the present paper is to clear up the problem of the Coulomb-gauge vector potential This paper is written by V Hnizdo in his private capacity. No official support or endorsement by the Centers for Disease Control and Prevention is intended or should be inferred.

2 of a uniformly moving charge using Jackson s results as well as the method employed by Onoochin, in which the difference between the Coulomb- and Lorenz-gauge vector potentials in this problem is to be found by solving a Poisson euation. In section, we obtain the difference between the potentials in the Lorenz and Coulomb gauges for a uniformly moving point charge using the formalism of Jackson s paper [6], which guarantees that the two gauges will yield the same electric and magnetic fields. In section 3, we demonstrate that the expression obtained for the difference between the vector potentials in the two gauges satisfies the inhomogeneous wave euation for this difference, and in section 4 we show where and why the procedure used by Onoochin for solving that wave euation by reducing it to a Poisson euation went wrong. It will turn out that the right-hand side of this euation must include an easily overlooked delta-function term to yield the correct solution. Concluding remarks are made in section 5, and an appendix contains some calculational details.. Transformation from the Lorenz gauge to the Coulomb gauge The scalar potential V and the vector potential A of a uniformly moving point charge in the Lorenz gauge, defined by the condition A+ A/c t =, are well known. For a point charge moving with a constant velocity v = vˆx, the Lorenz-gauge potentials V L and A L are given by (we shall use the Gaussian system of units): V L (r, t) = (x vt) + γ (y + z ) A L (r, t) = v c V L(r, t) () where γ = ( v /c ) / and, for simplicity, the charge is assumed to pass through the origin x=y=z= at time t = (see, e.g., [3], section 9-3). In the Coulomb gauge, defined by the condition A =, the scalar potential V C of the charge takes a particularly simple form, V C (r, t) = (x vt) + y + z () as the scalar potential in this gauge is exactly the same as that of the instantaneous Coulomb interaction of electrostatics. On the other hand, the Coulomb-gauge vector potential A C is the retarded solution to a relatively complicated inhomogeneous wave euation A C = 4π c vδ(r vt) + V C c t (3) where = /c t is the d Alembertian operator and vδ(r vt) is the point-charge current density of the present problem (see, e.g., [], section 6.3). The gauge invariance of electrodynamics implies that there is a gauge function χ C that connects the Coulomb- and Lorenz-gauge potentials by V C = V L χ C c t A C = A L + χ C (4) which ensures that the Lorenz-gauge and Coulomb-gauge potentials will yield the same electric and magnetic fields. This is because the fields are generated from the potentials via the prescription

3 A(r, t) E(r, t) = V (r, t) c t B(r, t) = A(r, t) (5) and thus any electric-field and magnetic-field differences that could arise from the use of the different gauges are guaranteed to vanish: (V C V L ) (A C A L ) = χ C c t c t c t χ C = (6) (A C A L ) = χ C =. (7) Nevertheless, it should be instructive to demonstrate explicitly that it is indeed the case also in the present problem by finding the reuisite gauge function. Before we turn to this task, we give here for completeness the fields that the prescription (5) yields with the Lorenz-gauge potentials V L and A L of euation (): γ(r vt) E(r, t) = B(r, t) = v E(r, t). (8) [γ (x vt) + y + z ] 3/ c The same expressions for the electric and magnetic fields can be obtained also by Lorentztransforming the electrostatic Coulomb field of the charge from its rest frame to the laboratory frame. Jackson [6] has derived the following integral expression for the gauge function χ C in terms of the charge density: χ C (r, t) = c d 3 r R R/c dτρ(r, t τ) (9) where R = r r (a gauge function is defined to within an arbitrary additive constant, which we omit here). For a point charge moving with a constant velocity v along the x-axis, the charge density is This gives R/c ρ(r, t) = δ(r vt) = δ(x vt)δ(y)δ(z). () R/c dτρ(r, t τ) = δ(y )δ(z ) = v δ(y )δ(z ) R/c dτ δ[x v(t τ)] dτ δ[τ (t x /v)] = v δ(y )δ(z ){Θ[R/c (t x /v)] Θ[ (t x /v)]} = v δ(y )δ(z )[Θ(x x ) Θ(x vt)] () where Θ(x) is the Heaviside step function and [ x = x γ x vt + v (x vt) c + ] + z γ ) (y γ = () v /c. The gauge function (9) with the charge density () is thus 3

4 χ C (r, t) = c v vt dx Θ(x x ) Θ(x vt) (x x ) + y + z = c dx v x (x x ) + y + z = c [ arcsinh x vt v y + z arcsinh x x ]. (3) y + z Let us first check that the gauge function (3) yields the established difference V C V L between the scalar potentials in the two gauges. Using the identity (x x γ ) + y + z = v (x c (x vt) + vt) + + z γ ) (4) (y which follows from (), to simplify the result of the partial differentiation χ C /c t, we obtain V C V L = χ C c t = (x vt) + y + z (x vt) + (y γ + z ). (5) This is indeed the correct result [see euations () and ()]. Calculating the x-component of the difference A C A L between the vector potentials in the two gauges is now very simple because χ C / x = (/v) χ C / t on account of the dependence of the gauge function (3) on the variables x and t only through the combination x vt. We thus have A Cx A L x = χ C x = c χ C v c t = c v (V C V L ). (6) The y- and z-components of the difference A C A L = χ C are obtained by performing direct differentiations in a similar way to that of calculating the value (5) for the difference V C V L, yielding A C y A L y = χ C y = c y(x vt) v y + z (V C V L ) (7) A Cz A L z = χ C z = c v z(x vt) y + z (V C V L ). (8) These components have no singularities (they vanish at x vt = y = z = ). As A L y = A L z =, euations (7) and (8) also give the Coulomb-gauge components A C y and A C z themselves, respectively. It is instructive to perform the Lorentz transformation of the Coulomb-gauge fourpotential (V C, A C ) from the laboratory frame, where the charge moves with the constant velocity v = vˆx, to its rest frame (the primes denote the rest-frame uantities): V C = γ(v C va C x /c) (9) A C x = γ(a C x vv C /c) A C y = A C y A C z = A C z. () 4

5 We note first that with the x-component A C x of the vector potential in the Coulomb gauge given by A C x = c v (V C V L /γ ) () which follows from the expression (6) for the difference A C x A L x and the fact that A L x = (v/c)v L, euation (9) gives the rest-frame scalar potential V C as V C = γ V L = γ (x vt) + (y + z )/γ = x + y + z () where the Lorentz transformation of the coordinates, x = γ(x vt), y = y, z = z, is performed on the right-hand side. The rest-frame scalar potential V C is simply that of a point charge in electrostatics. We note also that because the Coulomb-gauge condition A = is not Lorentz invariant, the rest-frame vector potential A C is not divergenceless in the rest-frame variables x, y, z ; the potentials V C, A C are therefore no longer those of the Coulomb gauge. However, a direct calculation shows that the vector potential A C is irrotational, A C =, which expresses the fact that there is no magnetic field in the rest frame. Moreover, the vector potential A C is independent of the rest-frame time t, and thus the electric field in the rest frame is given only by E = V C, which yields correctly the electrostatic Coulomb field of a charge at rest. 3. Inhomogeneous wave euation for the vector-potential difference The difference A C A L must satisfy the inhomogeneous wave euation (A C A L ) = V C c t which is obtained by subtracting the wave euation for A L, (3) A L = 4π c vδ(r vt) (4) (see [], section 6.3) from the wave euation (3) for A C. It is straightforward to show that the x-component A C x A L x = (c/v)(v C V L ) of the difference A C A L indeed satisfies the x-component of the inhomogeneous wave euation (3): [ ] c v (V C V L ) = V C c t x. (5) The fact that euation (5) holds true follows directly from the wave euations and c v V C = V C c t x 4πc δ(x vt)δ(y)δ(z) (6) v V L = 4πδ(x vt)δ(y)δ(z). (7) 5

6 The wave euation (6) in turn holds true because of the facts that the d Alembertian = /c t, V C = 4πδ(x vt)δ(y)δ(z), and V C / t = v V C / x; the wave euation (7) embodies the fact that the Lorenz-gauge scalar potential V L is the retarded solution of the inhomogeneous wave euation with the right-hand side 4πδ(x vt)δ(y)δ(z). It is also straightforward to show that the y-component (7) and z-component (8) of the difference A C A L satisfy the inhomogeneous wave euation (3) by using the identity (fg) = g f + f g + f g f g c t c t and euation (5) in the evaluation of the reuisite derivatives; the singularities of the functions y(x vt)/(y +z ) and z(x vt)/(y +z ) at y = z = cannot introduce any deltafunction terms in the d Alembertians of the components (7) and (8) as the latter are functions with no singularities. 4. Poisson s euation for the vector-potential difference Onoochin [9] attempts to solve the x-component of the inhomogeneous wave euation (3) for the difference A C A L directly by reducing it to a Poisson euation, which is a method based on the fact that the space and time partial derivatives are not independent when the source term is moving uniformly (see [3], section 9-3). In the present case, the dependence on the variables x and t is only through the combination x vt, and thus the x-component of the difference A C A L can be written as A C x A L x = f(x vt, y, z), where the function f(x vt, y, z) satisfies an inhomogeneous wave euation f(x vt, y, z) = V C c t x = v (x vt) y z + 4πvδ(x vt)δ(y)δ(z) (9) c [(x vt) + y + z ] 5/ 3c that can be cast as a Poisson euation on the substitutions / x /c t = /γ x and x vt = χ/γ, where γ = ( v /c ) / : ( ) χ + y + f(χ/γ, y, z) = v z c χ /γ y z (χ /γ + y + z ) (8) + 4πvγδ(χ)δ(y)δ(z). (3) 5/ 3c The delta-function term on the right-hand side of euation (9) arises from the fact that V C /c t x = (v/c) V C / x and the delta-function identity [] x i x j r = 3x ix j r δ ij r 5 4π 3 δ ijδ(r) (3) where r = (x + x + x 3 )/, δ ij is the Kronecker delta symbol and δ(r) = δ(x )δ(x )δ(x 3 ) is the three-dimensional delta function. Strictly speaking, the first term on the right-hand side of euation (3) should be understood as lim a (3x i x j r δ ij )/(r + a ) 5/, where the limit is to be performed after an R 3 integration with a well-behaved test function. This limit is automatically implemented when the integration is done in spherical coordinates 6

7 and the integration over the angular variables is performed first (without any additional transformation of variables that would shift the origin, of course) [,]. Unlike the well-known delta-function identity (/r) = 4πδ(r), the identity (3) is needed only relatively rarely in electromagnetism, an example being the calculation of the fields of electric and magnetic dipoles from their potentials (see []; and [], euations (4.) and (5.64)). The identity (/r) = 4πδ(r) obviously reuires that the second-order partial derivatives r / x i, i =,, 3 are given by expressions like those of euation (3); the precise form of these can be seen most easily to follow from the limit a of x i x j r + a = 3x ix j r δ ij (r + a ) 5/ a δ ij (r + a ) 5/ (3) as here the limit a of the second term on the right-hand side is a representation of 4 3 πδ ijδ(r). The standard integral expression for the solution to Poisson s euation (3) is given by f = v [ χ dχ dy dz /γ y z 4πc (χ /γ + y + z ) 4π ] 5/ 3 γδ(χ )δ(y )δ(z ) (33) (χ χ ) + (y y ) + (z z ) where one recovers the original variables by putting χ = γ(x vt). Here, the delta-function term is readily integrated to yield f = v dχ dy dz χ /γ y z 4πc (χ /γ + y + z ) 5/ (χ χ ) + (y y ) + (z z ) vγ 3c χ + y + z. (34) A direct evaluation in closed form of the three-dimensional integral in (34) does not seem possible. However, this integral can be evaluated for the special case y = z =, and the result is (see Appendix): v dχ dy dz χ /γ y z 4πc (χ /γ + y + z ) 5/ (χ χ ) + y + z = vγ 3c χ = v 3c x vt. (35) For y = z =, the value of the second term on the right-hand side of (34) is exactly eual and opposite to the value of (35), and thus the solution (34) of the Poisson euation (3) vanishes at y = z =, as reuired by the result (6) that we obtained for the difference A C x A L x. Onoochin s calculation omits the delta-function term that arises from the second-order partial derivative of the Coulomb-gauge scalar potential V C, and he takes erroneously just the first term of (34) as the solution of the Poisson euation (3). According to euation (35), this incorrect solution is non-zero at y = z =, and, moreover, its partial time derivative 7

8 does not vanish there. As (V C V L )/ x = at y = z =, the result (35), if it were the true difference A Cx A Lx at y = z =, would lead there to a non-zero difference between the Coulomb- and Lorenz-gauge x-components of the electric field. This is Onoochin s evidence against the euivalence of the Lorenz and Coulomb gauges. (In his paper [9], Onoochin does not evaluate the integral in euation (34), but he has communicated to the present author a calculation of its value at y=z= that differs from the value given by euation (35) by a factor of 3 [3].) It is not difficult to show directly that f = c v χ /γ + y + z γ (36) χ + y + z is the vanishing-at-infinity solution of the Poisson euation (3). Using the relation (3) and the differentiation rule g(ax, y, z)/ x = a g(u, y, z)/ u u=ax, we have χ χ /γ + y + z = [ χ /γ y z γ (χ /γ + y + z ) 4π ] 5/ 3 γδ(χ)δ(y)δ(z) y (37) χ /γ + y + z = y χ /γ z (χ /γ + y + z ) 4π γδ(χ)δ(y)δ(z) (38) 5/ 3 z χ /γ + y + z = z χ /γ y (χ /γ + y + z ) 4π γδ(χ)δ(y)δ(z) (39) 5/ 3 which gives ( χ + y + z ( = ) χ /γ + y + z γ ) χ /γ y z (χ /γ + y + z ) 5/ 4π 3 ( + γ ) γδ(χ)δ(y)δ(z). (4) Using the relation (/r) = 4πδ(r), we have also ( χ + y + z ) γ = 4πγδ(χ)δ(y)δ(z). (4) χ + y + z Subtracting euation (4) from euation (4), multiplying the result by c/v and using the fact that γ = v /c, we obtain ( χ + y + z ) c v χ /γ + y + z γ χ + y + z = v χ /γ y z + 4πvγδ(χ)δ(y)δ(z) (4) c (χ /γ + y + z ) 5/ 3c which shows that the function f of euation (36) indeed satisfies the Poisson euation (3). Transforming back to the original variables through χ = γ(x vt), the function f becomes 8

9 f = c v (x vt) + y + z = c (x vt) + (y γ + z ) v (V C V L ) (43) which is the value (6) for the difference A C x A L x. The function f of euation (36) has to eual the integral solution (34), yielding for the integral in (34) a closed-form expression dχ dy dz χ /γ y z 4π (χ /γ + y + z ) 5/ (χ χ ) + (y y ) + (z z ) = c v ( χ /γ + y + z v 3c ) γ. (44) χ + y + z Direct numerical three-dimensional uadrature of this integral, with the first factor of the integrand regularized as mentioned above in connection with euation (3), resulted in values that were close to those of the closed-form expression on the right-hand side of (44). 5. Concluding remarks We found an explicit expression for the gauge function of the transformation between the Lorenz and Coulomb gauges for a uniformly moving point charge. The Coulomb-gauge potentials obtained using the gauge function are guaranteed to yield the same electric and magnetic fields as the well-known Lorenz-gauge potentials of the charge. The expression obtained for the difference between the vector potentials in the two gauges satisfies the Poisson euation to which the inhomogeneous wave euation for this difference reduces after a transformation of the variables. However, the right-hand side of the Poisson euation involves a delta-function term which has to be included in the integral expression for its solution. Although gauge invariance is a foregone conclusion in a gauge-invariant theory, we believe that an explicit demonstration of the euivalence of the Lorenz and Coulomb gauges in the basic problem of a uniformly moving charge was instructive. Acknowledgment The author acknowledges correspondence with V V Onoochin, whose determined objections against mainstream electrodynamics motivated the writing of this paper. Appendix After the transformation χ /γ = x, the integral (35) can be written as (vγ/c)i(γ, X), where I(γ, X) is the integral (we drop the primes on the integration variables): I(γ, X) = x y z dx dy dz (45) 4π (x + y + z ) γ 5/ (X x) + y + z where X = χ/γ and γ > are real parameters. Transforming from the Cartesian coordinates x, y, z to the spherical ones, r, θ, φ, we have 9

10 I(γ, X) = = π r dr dr sin θ dθ r (3 cos θ ) r 5 γ (X r cosθ) + r ( cos θ) dξ 3ξ r γ (X rξ) + r ( ξ ). (46) We integrate with respect to ξ first (this is needed to implement the limit a of the first term on the right-hand side of euation (3), see [,]): F(r, γ, X) = = ω 5 r 4 dξ 3ξ r [ γ (X rξ) + r ( ξ ) 3γω(A+B) + (+γ )(ω r 3γ X ) ln γω X+r ω r γ X γω X r +ω r γ X ] (47) where ω = γ A = (ω r 3γ X) X + r B = (ω r + 3γ X) X r. (48) The function F(r, γ, X) has the properties F(r, γ, X) = F(r, γ, X) lim F(r, γ, X) =. (49) r + It peaks at r = X, where its derivative with respect to r is discontinuous. The argument of the logarithm in (47) reduces to (γ X ωr)/(γ X +ωr) for r < X, and to (γ ω)/(γ+ω) = (γ ω) for r > X. The integration with respect to r is performed in two parts: I (γ, X) = I (γ, X) = X X dr F(r, γ, X) = ( + γ )( + γ ) + + γ ln(γ ω) 6γω 4 X ω 5 X (5) dr F(r, γ, X) = 3γ ω 4 X + γ ln(γ ω) ω 5 X (5) where again ω = (γ ) /. The whole integral I(γ, X) is thus I(γ, X) = I (γ, X) + I (γ, X) = 3γ X (5) and the integral (35) has the value (vγ/c)i(γ, X) = v/(3c X ) = vγ/(3c χ ). The results of this Appendix were obtained using the software system Mathematica [4], and were checked by performing numerical integrations. [] Jackson J D 999 Classical Electrodynamics 3rd edn (New York: Wiley) [] Griffiths D J 999 Introduction to Electrodynamics 3rd edn (Upper Saddle River, NJ: Prentice Hall)

11 [3] Panofsky W K H and Phillips M 96 Classical Electricity and Magnetism nd edn (Reading, MA: Addison-Wesley) [4] Landau L D and Lifshitz E M 976 The Classical Theory of Fields 4th revised English edn (Oxford: Pergamon) [5] Rohrlich F Causality, the Coulomb field, and Newton s law of gravitation Am. J. Phys Jefimenko O D Comment on Causality, the Coulomb field, and Newton s law of gravitation, by F. Rohrlich Am. J. Phys Rohrlich F Reply to Comment on Causality, the Coulomb field, and Newton s law of gravitation, by O. D. Jefimenko Am. J. Phys [6] Jackson J D From Lorenz to Coulomb and other explicit gauge transformations Am. J. Phys [7] Stewart A M 3 Vector potential of the Coulomb gauge Eur. J. Phys Hnizdo V 4 Comment on Vector potential of the Coulomb gauge Eur. J. Phys. 5 in press (3 arxiv.org e-print physics/398) [8] Jackson J D and Okun L B Historical roots of gauge invariance Rev. Mod. Phys [9] Onoochin V V On non-euivalence of Lorentz and Coulomb gauges within classical electrodynamics Ann. Fond. L. Broglie ( arxiv.org e-print physics/7) [] Labarthe J-J 999 The vector potential of a moving charge in the Coulomb gauge Eur. J. Phys. L3 L3 [] Frahm C P 983 Some novel delta-function identities Am. J. Phys [] Aguirregabiria J M, Hernándes A and Rivas M δ-function converging seuencies Am. J. Phys [3] Onoochin V V 3 personal communication [4] Wolfram S 999 The Mathematica Book 4th ed (Champaign, IL: Wolfram Media)