The Dirac Equation for Many-Electron Systems
Total Page:16
File Type:pdf, Size:1020Kb
THE DIRAC EQUATION FOR MANY-ELECTRON SYSTEMS BY K. S. VISWANATHAN, F.A.SC. (Memoir No. 122 of the Raman Research Institute, Bangalore-6) Received June 11, 1960 1. INTRODUCTION THE well-known relativistic wave equation for an electron' was derived by Dirac by factorising a second order equation into two linear equations. It is the object of the present note to point out that the Dirac equation is simply the eigenvalue equation of the magnitude of the momentum four- vector, and that one can derive it by expressing the magnitude of the momentum vector in terms of its four components. The above idea enables one to generalise the Dirac equation and obtain a relativistic equation for systems containing several electrons. In Section 3, we have given the wave equation (Equation 16) for a system composed of several particles and this is very similar in form to the Dirac equation. Relativistic wave equations for a system of two electrons have previously been given by Eddington, 2 Gaunt 3 and Breit' of which the one given by Breit is the most satisfactory. By replacing the velocities vl and v" of the electrons by the spin matrices — cal and — call in a Hamiltonian given by Darwin, 4 Breit was able to obtain an approximate wave equation for two-electron systems. It is shown in Section 4, that Equation (16) leads to the Breit equation when it is represented in the product space of the two electrons. 2. THE MOMENTUM FOUR-VECTOR Before proceeding further, we first state a result which was first proved by Weyl6 for a vector in a n-dimensional space and which we shall apply presently. (a) Lemma.—Let (x1, x2 , • .., x,,) be the co-ordinates of a vector r in an Eucledian space with reference to a system of orthogonal axes, and let r = (x1 2 + x22 + ... + xn 2)' denote the 'magnitude' or 'length' of the vector r. Then 35 36 K. S. VISWANATHAN r = E (v^xx w=^ ), (1) where the y's are elements of an abstract algebra satisfying the relations (y,y,+yyµ) = 2S, . The y's can be expressed as matrices and for the case n = 4, they are the Dirac matrices. (b) The Dirac Equation.—Let x 1 , x2, x3 and x4 (= ict) be the co-ordinates of a world point, and similarly let pl, P2, p3 and p4 (= iE/c) denote the com- ponents of the momentum four-vector. Besides x 4 and p4i we shall also use the symbols xo and po given by x4 = ixo and p 4 = ip o. Expressed as operators we have then p 4 = — ihZ/^x 4 and po = i^i^/^xo. Now the momentum four-vector P is a vector with constant magnitude imoc where mo is the rest mass of the particle. Apply now the result (1) to the vector P = (pl, P2, P3, p4). We then get 4 E yip = imoc. (2) If ) is an eigenstate of the momentum four-vector, we get from (2) the equation of an electron as 3 iy4po + z yip — imoc) ) = 0. (3) i=1 For the y's we now choose the following representations :- I 0 O — ivlO — la2 Y4= Y1=C ig1 0 ), v2—_ i0 0 I); 20), 0 — 1O3 Y3 - i93 O ) where I, a,, v2 and a3 denote respectively the unit matrix in two dimensions and the three Pauli matrices. Multiplying (3) to the left by — 1V4 1' we get 3 (loo + j' aipi + J3moc 1 0) = 0(4) which is the Dirac equation in its conventional form. When the electron is moving in a field, the components of the momentum four-vector are given by (pi + e/c Ai) (i = 0, 1, 2, 3) where Ao, Ai , A2 , A3 The Dirac Equation for Many-Electron Systems 37 4 are the scalar and vector potentials of the field. From the relations ^ (pi + e/c Ai) 2 = — m02c2, we see that the magnitude of the momentum four- vector is equal to imoc in this case also. Thus by replacing pi by (pi + e/c Ai) in (4), we get the equation for an electron moving in field as 3 ee j(PO + e Ao l + ai (pi + e Ai) + flmoc 10) = 0. (5) The left-hand side of (2) is similar in form to the expression of the length of a vector in terms of its direction cosines. The i's can thus be regarded as the representations in a matrix algebra of the direction cosines of the momentum four-vector. We have thus Yi= — Pz .^ or iyi= ui, (6 ) IPI where ui (i = 1, 2, 3, 4) are the components of the velocity four-vector. The matrices iyi thus represent the components of the velocity four-vector. Since ai= — i(Y4)-1 Yi(i=1, 2 ,3), we have ai= _ZPi = _xi P4 c or Xi = — Cai. (7) Similarly 2 ^3 = — 1 — 2(8) We thus get the well-known expressions for the components of the velocity of the particle without calculating the commutation relations of x1 , x2 and x 3 with the Hamiltonian. In finding out the dynamical variables that are the classical analogues of products (or quotients) of matrices, care should be taken to verify that only the algebraical rules that are common to both the matrix and the ordi- nary algebras are used. We have derived (4) by multiplying (2) by — iy4 1 making use of the relation (y 4)_' Y 4 = 1 which is common to both algebras, 38 K. S. VISWANATHAN though one could equally well derive (4) by multiplying (2) by — iy4. The classical analogue of a2 is thus — i (y 4) - ' yi and not — iy 4 yz. 3. SYSTEM OF MANY ELECTRONS We have seen that the Dirac equation can be derived by expressing the magnitude of the four-vector P + e/c A in terms of its components. Now (Pi + Ati) = mocui (i = 1, 2, 3, 4) mo z^, (9) / v/1 2, C2 where u1 , u 2i u 3 , u, are the components of the velocity four-vector. natural way to generalise the Dirac equation for a system of particles would be to consider a four-vector whose components are respectively E mocuz (i = 1, 2, 3, 4, or 0) here the sum is to be taken over all the particles'of the system. We shall denote the components of this vector by (Po + e/c A0), (P1 + e/c A,), (P2 + e/c A2) and (P3 + e/c A3). Now in the special theory of relativity the components of the total momentum are given by 7 P i = — -' f TIkdSk + Z rflocuz, c (10) where Tzk.. are the components of the energy-momentum tensor. In three- dimensional form, we can write for the total momentum of field plus charges f dV+Ep (10 a) and for the energy f W dV + 2 , (10 b) where S = C ExH is the Poynting vector and W = I (E" 2 + H2) The Dirac Equation for Many-Electron Systems 39 is the density of the field energy. For a system of charges, the electrostatic energy (U), apart from the self-energy terms, is equal to U= f WdV=E eAeB . rnB We thus see that the quantities — e (A0, Ai, A2, A3) stand respectively for the energy-momentum of the field. Let us denote the magnitude of P + e(c A by iP'. P' is equal to E moic only as a first approximation. We have in fact s 3 ti — P' 2 = — c$ (^' mi) z -f- ^ ( J mixik) $, (11) o k=1 i-1 where m0 nit - Thus - p, = — 2(,E Mj2 + 2 £ miml ,, r 3<1 +(,E mi2via + 2 mimivi . vi +<r or P12= E mi2c$ (1 — Vie, + 2 E mimlc2 ( i — V c ) ' (12) Now mimics 1 —v C2 ) =moi oic2 i_! ) J7ii 2 \( C21 _ vf\ J ( 2 . moimojc 1 + 2 2c2)' (13) 40 K. S. VISWANATHAN where vij = (vi — vi). (14) Substituting (13) in (12), we get P' 2 _ Em0i2c2 + 2 Em0im0jc 2 (1 +v 2 \ _ (E m0ic)2 + E m0im0jvij 2. Thus a second approximation for P' is given by n Emoimojvij P' ^ 2' moic + i- "(15) ' moic Applying now the results of Section 1 to the vector P + e/c A, we get the equation for the system of electrons as 3 (Po + A0 ) + ak (Pk + Ak) + 9P' I ) =0. (16) J k=1 4. REPRESENTATION IN THE PRODUCT SPACE AND THE BREIT EQUATION Equation (16) describes the system as a whole and can be considered to be the equation of motion for the centre of inertia of the system. Since we are operating in a four-dimensional space, the matrices ak and P are all four-dimensional. They refer to the entire-system of electrons and do not contain any labels of the individual particles. In practice, however, one needs equations that bring in explicitly the positions and spins of the individual electrons of the system. We shall see presently that Equation (16) can be transformed into one that contains explicit reference to the spins of the electrons if it is represented in the product space of the electrons. Since the y's represent the direction cosines of the momentum four- vector, their classical analogues are given by n ^' mizik = ¢ = 1(k = 1, 2, 3) Yk iP and (E imic) Y4 = p, (17 a) The Dirac Equation for Many-Electron Systems 41 Thus we have E mixik ak = — i (Y^)-1 Yk = — c E mi and (17 b) pp , _—P'2cEmi Thus we have 3 , L ak (Pk + c Ak)+flp k=1 _ l kL' (2' m2 tik) 2 — P12 1 c E mi i =— Af mic= —c f moi — C2 V(I ) ^ w 2 _ — moic 1 — c2 — C E miv42(18) icl V i=1 correct to terms of the order of v 4/c3 .