arxiv:math-ph/ v1 5 Mar 2007

Transcription

1 Submitted to J. Math. Phys. March 5, 2007 Conservation laws for the Maxwell-Dirac equations with a dual Ohm s law Nail H. Ibragimov and Raisa Khamitova Department of Mathematics and Science, Research Centre ALGA: Advances in Lie Group Analysis, Blekinge Institute of Technology, SE Karlskrona, Sweden Bo Thidé Swedish Institute of Space Physics, Ångström Laboratory, P. O. Box 537, SE Uppsala, Sweden Using a general theorem on conservation laws for arbitrary differential equations proved by Ibragimov, we have derived conservation laws for Dirac s symmetrized Maxwell-Lorentz equations under the assumption that both the electric and magnetic charges obey linear conductivity laws dual Ohm s law. We find that this linear system allows for conservation laws which are non-local in time. arxiv:math-ph/ v1 5 Mar 2007 I. INTRODUCTION PACS numbers: j, a,02.30.Jr,03.50.De,41.20.Jb In all areas of physics, conservation laws are essential since they allow us to draw conclusions of a physical system under study in an efficient way. Electrodynamics, in terms of the standard Maxwell electromagnetic equations for fields in vacuum, exhibit a rich set of symmetries to which conserved quantities are associated. Recently, there has been a renewed interest in the utilisation of such quantities. Here we use a theorem of Ibragimov 2006 to derive conservation laws for Dirac s symmetric version of the Maxwell-Lorentz microscopic equations, allowing for magnetic charges and magnetic currents, where the latter, just as electric currents, are assumed to be described by a linear relationship between the field and the current, i.e., an Ohm s law. The method of Ibragimov 2006 produces two new adjoint vector fields which fulfil Maxwell-like equations. In particular, we obtain conservation laws for the symmetrized electromagnetic field which are non-local in time. II. PRELIMINARIES A. Notation We will use the following notation see, e.g., Ref. Ibragimov, Let x = x 1,...,x n be independent variables and u = u 1,...,u m be dependent variables. The set of the first-order partial derivatives ui α = u α /x i will be denoted by u 1 = {ui α }, where α = 1,...,m and i, j,... = 1,...,n. The symbol D i denotes the total differentiation with respect to the variable x i : D i = x i + uα i u α + uα i j u j α + ui α jk u α jk + Electronic address: nib@bth.se Electronic address: bt@irfu.se Also at LOIS Space Centre, School of Mathematics and Systems Engineering, Växjö University, SE Växjö, Sweden We employ the usual convention of summation in repeated indices. Recall that a necessary condition for extrema of a variational integral V Lx,u,u 1 dx 1 with a Lagrangian Lx,u,u 1, depending on first-order derivatives, is given by the Euler-Lagrange equations δu α L L u α D i ui α = 0, α = 1,...,m. 2 We will understand by a symmetry of a certain system of differential equations a generator X = ξ i x,u x i + ηα x,u u α 3 of a continuous transformation group admitted by differential equations under consideration. A vector field C = C 1,...,C n is said to be a conserved vector for the differential equations 2 if the equation D i C i = 0 4 holds for any solution of Eq. 2. If one of the independent variables is time, e.g., x n = t, then the conservation law is often written in the form where E = de dt = 0, R n 1 C n x,ux,u 1 xdx 1 dx n 1. 5 Accordingly, C n is termed the density of the conservation law.

2 2 B. Basic conservation theorem We will employ the recent general theorem Ibragimov, 2006 on a connection between symmetries and conservation laws for arbitrary systems of sth-order partial differential equations F α x,u,u1,...,u s = 0, α = 1,...,m, 6 where F α x,u,u 1,...,u s involves n independent variables x = x 1,...,x n and m dependent variables u=u 1,...,u m, u=ux together with their derivatives up to an arbitrary order s. For our purposes, we formulate the theorem in the case of systems of first-order differential equations. Theorem II.1. See Ref. Ibragimov, 2006, Theorem 3.5. Let an operator 3 be a symmetry of a system of first-order partial differential equations F α x,u,u1 = 0, α = 1,...,m. 7 where v=v 1,...,v m. Then the quantities C i = v [ξ β i F β + ] η α ξ j u j α F β ui α, i=1,...,n, 8 furnish a conserved vector C=C 1,...,C n for the equations 7 considered together with the adjoint system where F α x,u,v,u1,v 1 = 0, α = 1,...,m, 9 δuα δ δu α = u α D i u α i, α = 1,...,m, and v = v 1,...,v m are are new dependent variables, i.e., v=vx. Remark II.1. The simultaneous system of equations 7 and 9 with 2m dependent variables u=u 1,...,u m, v=v 1,...,v m can be obtained as the Euler-Lagrange equations 2 with the Lagrangian L=v β F β x,u,u1,...,u s 10 Indeed, δv α = F α x,u,u1, 11 δu α = F α x,u,v,u1,v Remark II.2. The conserved quantities 8 can be written in terms of the Lagrangian 10 as follows: C i =Lξ i + η α ξ j u j α L ui α. 13 Remark II.3. If Eqs. 7 have r symmetries X 1,...,X r of the form 3, X µ = ξ i µ x,u x i + ηα µ x,u u α, µ = 1,...,r, then Eqs. 8 provide r conserved vectors C 1,...,C r with the components C i µ =Lξ i µ + ηµ α ξ µu j j α L ui α, µ = 1,...,r; i=1,...,n. III. ELECTROMAGNETIC EQUATIONS A. Basic equations and the Lagrangian Adopting Dirac s ideas on the existence of magnetic monopoles Dirac, 1931, one can formulate a symmetrized version of Maxwell s electromagnetic equations Schwinger, In SI units and in microscopic Lorentz form, these equations are [cf. Ref. Thidé, 2006, Eqs. 1.50]: E+ B t + µ 0j m = 0, B 1 E c 2 t µ 0j e = 0, E µ 0 c 2 ρ e = 0, B µ 0 ρ m = 0, together with the dual Ohm s law 14a 14b 14c 14d j e = σ e E, j m = σ m B. 15 where σ m and σ e are constant scalar rank zero quantities. The first equation in 15 is Ohm s law for electric currents. The second equation is a dual Ohm s law for magnetic currents, that, for symmetry reasons, was introduced in Ref. Thidé, 2006, Eqs. 2.60; see also Eq. 5 in Ref. Meyer-Vernet, 1982, Eq. 38 in Ref. Olesen, 1996, and its generalization Eq. 8 in Ref. Coceal et al., Now we substitute Eqs. 15 into Eqs. 14a 14b. The ensuing equations involve, along with the light velocity c, three other constants, σ e,σ m and µ 0. We eliminate two constants by setting t = ct, B=cB, σe = cµ 0 σ e, σ m = µ 0 c σ m, ρ e = c 2 µ 0 ρ e, ρ m = cµ 0 ρ m 16 and rewrite our basic Maxwell-Dirac equations 14, discarding tilde, as follows: E+ B t + σ mb=0, B E t σ ee=0, E ρ e = 0, B ρ m = 0, 17 The system 17 has eight equations for eight dependent variables: six coordinates of the electric and magnetic vector fields E=E 1,E 2,E 3 and B=B 1,B 2,B 3, respectively, and two scalar quantities, viz., the electric and magnetic monopole charge densities ρ e and ρ m.

3 3 Using the method of Ibragimov 2006 we write the Lagrangian 10 for Eqs. 17 in the following form: L=V E+ B t + σ mb +R e E ρ e + W B E t σ ee +R m B ρ m, 18 where V,W,R e,r m are adjoint variables we note in passing that V is a pseudovector and R m a pseudoscalar. With this Lagrangian we have: and B = E+ δv t + σ mb, E = B δw t σ ee, W = V+ σ e W R e, δe t V = W δb t + σ mv R m, = E ρ e, δr e = B ρ m, δr m = R e, δρ e = R m. δρ m It follows from Eqs that the Euler-Lagrange equations 2 for the Lagrangian 18 provide the electromagnetic equations 17 and the following adjoint equations for the new dependent variables V,W,R e,r m : V+ W t σ e W=0, W V t + σ mv=0, R e = 0, R m = Remark III.1. Let the spatial coordinates x 1,x 2,x 3 be x,y,z. For computing the variational derivatives /δe and /δb in Eqs. 20, it is convenient to use the coordinate representation of the Lagrangian 18, namely: L= V 1 E 3 y E2 z + B1 t + σ mb 1 +V 2 E 1 z E3 x + B2 t + σ mb 2 +V 3 E 2 x E1 y + B3 t + σ mb 3 +R e E 1 x + E2 y + E3 z ρ e 22 +W 1 B 3 y B 2 z E 1 t σ e E 1 +W 2 B 1 z B 3 x E 2 t σ e E 2 +W 3 B 2 x B 1 y E 3 t σ e E 3 +R m B 1 x+ B 2 y+ B 3 z ρ m. B. Symmetries Eqs. 17 are invariant under the translations of time t and the position vector x = x,y,z as well as the simultaneous rotations of the vectors x,e and B due to the vector formulation of Eqs. 17. These geometric transformations provide the following seven infinitesimal symmetries: X 0 = t, X 1 = x, X 2 = y, X 3 = z, X 12 = y x x y + E2 E 1 E1 E 2 + B2 B 1 B1 B 2, X 13 = z x x z + E3 E 1 E1 E 3 + B3 B 1 B1 B 3, X 23 = z y y z + E3 E 2 E2 E 3 + B3 B 2 B2 B The infinitesimal symmetries for the adjoint system 21 are obtained from 23 by replacing the vectors E and B by V and W, respectively. Moreover, since Eqs. 17 are homogeneous, they admit simultaneous dilations of all dependent variables with the generator where T = E E + B B + ρ e + ρ m, 24 ρ e ρ m E E = 3 i=1 E i 3 E i, B B = i=1 B i B i. Recall that the Maxwell equations in vacuum admit also the one-parameter group of Heaviside-Larmor-Rainich duality transformations E=Ecosα Bsinα, B=Esin α+ Bcosα 25 with the generator 3 X = E B B E i=1 E i B i Bi E i. Also recall that the mixing angle α in 25 is a pseudoscalar. It was shown in Ibragimov, 2006 that the group 25 provides the conservation of energy for the Maxwell equations. Let us clarify whether Eqs. 17 admit a similar group. Let therefore X = E B B E + ρ e ρ m. 26 ρ m ρ e The prolongation of the operator 26 is written X = E B B E + ρ e ρ m ρ m ρ e + E t B t + E x B x B t E t B x + E y B y B y E y + E z B z B z E x. 27 E z Reckoning shows that the operator 27 acts on the left-hand sides of Eqs. 17 as follows: X E+B t + σ m B = B E t σ m E, X B E t σ e E = E+B t + σ e B, X E ρ e = B ρ m, X B ρ m = E ρ e.

4 4 It follows that the operator 26 is admitted by Eqs. 17 only in the case IV. CONSERVATION LAWS A. Derivation of conservation laws σ m = σ e. 28 We will write the conservation law 4 in the form D t τ+div χ = 0, 29 where the pseudoscalar τ is the density of the conservation law 29, the pseudovector current χ =χ 1, χ 2, χ 3, and div χ χ = D x χ 1 +D y χ 2 +D z χ 3. Let us find the conservation law furnished by the symmetry 26 when the condition 28 is satisfied, σ m = σ e. Applying the formula 8 to the symmetry 26 and to the Lagrangian 18, we obtain the following density of the conservation law 29: Thus, τ = E L B t B L E t = E V+B W. τ = E V+B W. 30 The pseudovector χ is obtained likewise. For example, using the Lagrangian in the form 22, we have: χ 1 = E L B x B L E x = E 2 W 3 E 3 W 2 B 2 V 3 + B 3 V 2. The other coordinates of χ are computed likewise, and the final result is: χ =E W B V. 31 One can readily verify that 30 and 31 provide a conservation law for Eqs. 17 considered together with the adjoint equations 21. Indeed, using the well-known formula a b=b a a b and Eqs. 17 and 21, we obtain: Whence, D t τ=e t V+E V t + B t W+B W t = V B σ e E+E W+σ m V W E+σ m B B V σ e W, χ = E W B V = W E E W V B+B V. D t τ+div χ =σ m σ e E V B W. It follows again that the conservation law is valid only if σ m = σ e. Remark IV.1. The conservation law given by depends on solutions V,W of the adjoint system 21. However, substituting into Eqs. 30 and 31 any particular solution V,W of the adjoint system 21 with σ m = σ e, one obtains the conservation law for Eqs. 17 not involving V and W. Let us denote σ m = σ e = σ and take, e.g., the following simple solution of the adjoint system 21: V 1 = e σt, V 2 = V 3 = 0; W 1 = e σt, W 2 = W 3 = 0. Then Eqs yield: τ =E 1 + B 1 e σt χ 1 = 0 χ 2 =E 3 B 3 e σt χ 3 =B 2 E 2 e σt. Remark IV.2. The operator 26 generates the one-parameter group E=Ecosα Bsinα, B=Esin α+ Bcosα, ρ e = ρ e cosα ρ m sinα, ρ m = ρ e sinα+ ρ m cosα. where, again, the mixing angle α is a pseudoscalar. Remark IV.3. In the original variables used in Eqs. 14a 14b and 15, the operator 26 is written: X = 1 c E B cb E + cρ e 1 ρ m c ρ m. ρ e Applying similar calculations to the generator 24 of the dilation group provides the conservation law with τ = B V E W, χ =E V+B W. 32 This conservation law is valid for arbitrary σ m and σ e. Indeed, D t τ=b t V+B V t E t W E W t = V E+σ m B+B W+σ m V W B σ e E+E V σ e W = V E+B W W B+E V, χ = E V+ B W = V E E V + W B B W. Hence, D t τ+div χ = 0. Let us find the conservation law provided by the symmetry X 0 = /t from 23. Formula 8 yields: τ =L E t L E t B t L B t =L+E t W B t V. Since the Lagrangian L given by 18 vanishes on the solutions of Eqs. 17, we can take or τ = E t W B t V 33 τ = W [ B σ e E]+V [ E+σ m B]. 34

5 5 Let us calculate the pseudovector χ. Formula 8 yields: χ 1 = E t L E x B t L B x. Using the Lagrangian in the form 22, we have: χ 1 = E 2 t V 3 + E 3 t V 2 B 2 t W 3 + B 3 t W 2. The other coordinates of χ are computed similarly, and the final result is χ =V E t +W B t. 35 Thus, the time translational invariance of Eqs. 17 leads to the conservation law 29 with τ and χ given by 34 and 35, respectively. Remark IV.4. Let us substitute in Eqs. 34 and 35 the following simple solution of the adjoint system cf. Remark IV.1: V 1 = e σ mt, V 2 = V 3 = 0; W 1 = e σ et, W 2 = W 3 = 0. Then Eqs yield: τ =B 3 y B 2 z σ e E 1 e σ et +E 3 y E 2 z + σ m B 1 e σ mt, χ 1 = 0, χ 2 = E 3 t eσ mt B 3 t eσ et, χ 3 = E 2 t eσ mt + B 2 t eσ et. The conservation law provided by the symmetry X 1 = /x from 23 has the following density: τ = E x L E t B x L B t = E x W B x V. For the pseudovector χ the formula 8 yields: χ 1 =L E x L E x B x L B x. Using the Lagrangian in the form 22, we have: χ 1 =L E 2 x V 3 + E 3 x V 2 B 2 x W 3 + B 3 x W 2. The other coordinates of χ are calculated similarly: χ 2 = E 1 xv 3 E 3 xv 1 + B 1 xw 3 B 3 xw 1, χ 3 = E 1 xv 2 + E 2 xv 1 B 1 xw 2 + B 2 xw 1. We can ignorelin χ 1 since D x L=0on solutions of Eqs. 17 and 21, and the final result is χ =V E x +W B x. 36 Replacing x by y and z we obtain the following conservation laws corresponding to X 2 = /y and X 3 = /z, respectively: and τ = E y W B y V, χ =V E y +W B y τ = E z W B z V, χ =V E z +W B z. Applying Formula 8 to the symmetry X 12 and to the Lagrangian 22, we obtain the following density of the conservation law: τ =E 2 L Et 1 E 1 L Et 2 +xe y ye x L E t + B 2 L Bt 1 B 1 L Bt 2 +xb y yb x L B t =W 2 E 1 W 1 E 2 +ye x xe y W V 2 B 1 V 1 B 2 yb x xb y V. The densities of the conservation laws for X 13 and X 23 are and τ = W 3 E 1 W 1 E 3 +ze x xe z W τ = W 3 E 2 W 2 E 3 +ze y ye z W V 3 B 1 V 1 B 3 zb x xb z V V 3 B 2 V 2 B 3 zb y yb z V, respectively. Finally, the densities of conservation laws corresponding to the rotation generators X i j can be written as one vector: τ = W E+W x E V B V x B, 37 where x=x,y,z. The operator X 12 provides the following pseudovector χ: χ 1 = V 3 E 1 ye 2 xv 3 E 3 xv 2 +xe 2 yv 3 E 3 yv 2 W 3 B 1 yb 2 x W 3 B 3 x W 2 +xb 2 y W 3 B 3 y W 2, χ 2 = V 3 E 2 + ye 1 xv 3 E 3 xv 1 xe 1 yv 3 E 3 yv 1 W 3 B 2 + yb 1 xw 3 B 3 xw 1 xb 1 yw 3 B 3 yw 1, χ 3 =V 1 E 1 +V 2 E 2 ye 1 x V 2 E 2 x V 1 +xe 1 y V 2 E 2 y V 1 +W 1 B 1 +W 2 B 2 yb 1 xw 2 B 2 xw 1 +xb 1 yw 2 B 2 yw 1. The pseudovector χ for the operator X 13 has the following form: χ 1 =V 2 E 1 ze 2 x V 3 E 3 x V 2 +xe 2 z V 3 E 3 z V 2 +W 2 B 1 zb 2 x W 3 B 3 x W 2 +xb 2 z W 3 B 3 z W 2, χ 2 = V 1 E 1 V 3 E 3 + ze 1 x V 3 E 3 x V 1 xe 1 z V 3 E 3 z V 1 W 1 B 1 W 3 B 3 + zb 1 xw 3 B 3 xw 1 xb 1 zw 3 B 3 zw 1, χ 3 =V 2 E 3 ze 1 xv 2 E 2 xv 1 +xe 1 z V 2 E 2 z V 1 +W 2 B 3 zb 1 x W 2 B 2 x W 1 +xb 1 z W 2 B 2 z W 1. The operator X 23 provides the following pseudovector χ: χ 1 =V 2 E 2 +V 3 E 3 ze 2 yv 3 E 3 yv 2 +ye 2 z V 3 E 3 z V 2 +W 2 B 2 +W 3 B 3 zb 2 y W 3 B 3 y W 2 +yb 2 z W 3 B 3 z W 2, χ 2 = V 1 E 2 + ze 1 y V 3 E 3 y V 1 ye 1 z V 3 E 3 z V 1 W 1 B 2 + zb 1 yw 3 B 3 yw 1 yb 1 zw 3 B 3 zw 1, χ 3 = V 1 E 3 ze 1 y V 2 E 2 y V 1 +ye 1 z V 2 E 2 z V 1 W 1 B 3 zb 1 y W 2 B 2 y W 1 +yb 1 z W 2 B 2 z W 1.

6 6 B. Two-solution representation of conservation laws The conserved quantities obtained in Subsection IV.A involve solutions V, W of the adjoint equations 21. It may be useful for applications to give an alternative representation of the conserved quantities in terms of the electric and magnetic vector fields E,B only. We suggest here one possibility based on the observation that one can satisfy the adjoint system 21 by letting Vx,t=Bx, t, Wx,t=Ex, t, R e x,t= Ex, t ρ e x, t, R m x,t= Bx, t ρ m x, t, 38 where Ex,s,Bx,s solve Eqs. 17 with s = t. Indeed, employing the substitution 38 and the notation s = t we have V+ W t σ e W = Bx,s+ Ex,s s W V t + σ mv = Ex,s Bx,s s R e = Ex,s ρ e x,s=0, R m = Bx,s ρ m x,s. s t σ e Ex,s, s t + σ mbx,s, Hence, the adjoint equations 21 reduce to 17: Ex,s+ Bx,s + σ m Bx,s=0, s Bx,s Ex,s σ e Ex,s=0, s Ex,s ρ e x,s=0, Bx,s ρ m x,s= Let Ex,t,Bx,t and E x,t,b x,t be any two solutions of the electromagnetic equations 17. Substituting in 38 the solution E,B, we obtain the two-solution representations of the conservation laws. For example, the conservation law given by has in this representation the following coordinates: τ = Ex,t B x, t+bx,t E x, t, χ =[Ex,t E x, t] [Bx,t B x, t]. 41 In particular, if the solutions Ex,t,B,t are identical, 41 provides the one-solution representation: τ = Ex,t Bx, t+bx,t Ex, t, χ =[Ex,t Ex, t] [Bx,t Bx, t]. 42 All other conservation laws can be treated likewise, e.g., the conservation law given by 33 and 35 has the following two-solution representation: τ = E t x,t E x, t B t x,t B x, t, χ =[B x, t E t x,t]+[e x, t B t x,t]. Acknowledgments 43 One of the authors B. T. gratefully acknowledges the financial support from the Swedish Governmental Agency for Innovation Systems VINNOVA. References N. H. Ibragimov, J. Math. Anal. Appl. in press 2006, doi: /j.jmaa N. H. Ibragimov, Elementary Lie group analysis and ordinary differential equations John Wiley & Sons, Chichester, 1999, 2nd ed. P. A. M. Dirac, Proc. Roy. Soc. London 133, J. Schwinger, Science 165, B. Thidé, Electromagnetic Field Theory Upsilon Books, Uppsala, Sweden, 2006, URL N. Meyer-Vernet, Am. J. Phys. 50, P. Olesen, Phys. Lett. B 366, O. Coceal, W. A. Sabra, and S. Thomas, Europhys. Lett. 35,