Viewing message <a7v373$udo$1@rs04.hrz.uni-essen.de>
From: Urs Schreiber (Urs.Schreiber@uni-essen.de)
Subject: Re: Physically understanding the Dirac equation and
4D.
Newsgroups: sci.physics.research
Date: 2002-03-28 18:06:52 PST
"Roman Arce" <shampoo@fibertel.com.ar> schrieb im Newsbeitrag
news:3ec0ba0c.0203261529.4352333e@posting.google.com...
> Does anybody have any idea of how to get a more physical understanding
> of Dirac's equation? I mean I can mess with the Dirac gamma matrices
> and do the math (specially with the help of the software Mathematica)
> but always loosing complete contact with the physical reality. Having
> four-vectors made of 4 components and each component being just a
> number it's physically understandable to me but when you need to
> replace each one of the 4 numbers with a 4x4 matrix it becomes
> something that can be played mathematically but it's ununderstandable
> physically (to me at least).
> Can the Dirac equation be understood physically, graphically, or
I'm
> asking for something impossible like a human brain really watching
4
> dimensions graphically? With 4D for most things it's easy cos you
can
> always think of only 1 space dimension and you're only left with
2D
> which is obviously easy to understand graphically and therefore
> physically but with the Dirac matrices you have spinors so you're
> stucked with 3D and with time it's 4D so any graphical thing is no
> longer an option, but could we still have a physical understanding
> some other way? I'm sure that if our brains could 'see' 4D all this
> wouldn't be a problem and the physical understanding would be easy
(or
> at least possible).
> Any idea about how to 'see' the Dirac equation physically and any
help
> about this subject is welcomed.
Dirac operators are quantized connections.
That's due to Quillen(see [2]).This simple statement contains in it
the
key to all the intuition you may want about the meaning of Dirac
operators and Dirac equations in physics or elsewhere. I'll try
to give a first order approximation to an explanation of this
statement. But I will run out of time. To be continued...
First of all: I think your question is great. Almost all QM textbooks
give a treatment of the Dirac equation along the lines of Dirac's
original presentation. That was certainly a marvelous achievement,
then, but there has been a huge amount of progress - at least in the
mathematically oriented literature. Dirac operators are now recognized
as being intimately related to (semi-)Riemannian geometry and the
geometry of fiber bundles - that is: to the Geometry of Physics.
In fact, there is an approach due to Alain Connes [5] that goes the
other way round: Instead of starting with a geometry and deriving
a Dirac operator from there he starts with the algebra generated by
some function on a would-be manifold and an abstract operator
D. One can show that these ingredients alone define Riemannian
geometry as well as vast generalizations thereof. Connes identifies
the
Dirac operator formally with the inverse of the infinitesimal length
element
D = 1/ds .
This may look a little hand-wavy, but can be given rigorous meaning
(see [5]). It shows what status the Dirac operator really has: It is
as fundamental as the notion of distance.This approach is immensely
powerful. (Very recently it has been generalized to semi-Riemannian
geometry: math-ph/0110001.)
But you write:
> I can mess with the Dirac gamma matrices [...]
> but always loosing complete contact with the physical reality .
That's a pity, because Dirac theory can get you much closer to
"physical reality." But it's probably not unusual, either, since the
underlying structure of Clifford algebra and Dirac theory is hardly
visisble in matrix representation. The key to a better understanding
is not to work (or at least not to think) in terms of explicit matrices
but in terms of "abstract" Clifford algebra. But don't worry, you'll
soon realize that that's far less "abstract" than you might think!
Quite the opposite: Doing away with the matrices reveals a
notational calculus well adapted to describe geometric and thus
physical notions.
At this point I have to mention the school of thought called
"Geometric Algebra". This originates with David Hestenes, who,
as a student, was troubled with similar questions as you raised
in your post. To his great delight (he describes that in some paper,
but I forget which one) he discovered how everything begins
to make sense when Clifford algebra is properly dealt with in
physics. He became very enthusiastic about this discovery and
started to rewrite huge wads of physics in a notation based on
Clifford algebra. (It might have helped if he had been a little less
enthusiastic.) There are lots of things to know and read about
Dirac theory, but you can give yourself a head start by first of all
looking at Hestenes' original book [6]. (A new version is in
preparation and online available [8]. Also see [7].) I'd say when
you feel familiar with his presentation, *then* look at the classic
texts on spinors and geometry [1]. With math textbooks it is
always helpful to know beforehand what they are *really* talking
about, because sometimes they are too shy to tell you frankly!
:-)
I recommend to go to the Geometric Algebra/Geometric Calculus
homepage (I do not have the address at hand right now, search for
it via google) and have a look at the information to be found there.
The "geometric algebra"-school is often accused for overemphasizing
the importance of their work. That's because their practitioners are
so fond of the new insight they gain by rewriting existing results
in
a neater way. But for someone troubled by questions like you asked
above this is a great source of information.
Now that I have said all this in order to put your question and my
attempt at an answer in proper perspective, I should finally begin
with answering:
> Can the Dirac equation be understood physically, graphically,
Yes. But there may be several ways to look at it. One way is
this (from a recent thread called "Dirac for Dunces")
---- post from the thread "Dirac for Dunces":
"John Baez" <baez@galaxy.ucr.edu>, wrote in
news:a1b5fc$fcn$1@glue.ucr.edu...
[...]
> "But if he doesn't remember his tensors, how in the world will he
> ever understand spinors."
Apropos "understanding spinors":
Maybe it helps to chat about some fascinating and intuitively
accessible details (it sure helps me): Mathematically, there is a
close relation between spinors and reflections. The fascinating
implication for physics is that discrete checkerboard models,
where massive particles move around at the speed of light but
reverse direction (reflection!) every now and then, actually
describe *spinning* particles!
How does that come about?
The best way to talk about spinors and reflections is, once again,
Clifford algebra, so I'll stick to that. Adopting Danny Ross
Lunsford's ASCII notation for Clifford algebra I call
yi
the i-th unit vector with square
yi yi = - 1 .
(the negative sign is conventional, see below)
The (Clifford) product of two such unit vectors
yi yj ,
called a bivector, is, as the Wizard has already pointed out,
the parallelogram spanned by these two vectors. Since
this is supposed to change sign when its orientation is reversed
we have
yi yj = - yj yi (if i is different from j).
The first fact to note here is that every vector, say y0, carries
with it the instruction to reflect any other vector at the
(hyper-)plane orthogonal to y0 by applying y0 from both
sides:
A simple example makes this clear: Consider the
2-dimensional euclidean plane with unit vectors y0 and y1.
An arbitrary vector v can be written as a linear combination
of these
v = v0 y0 + v1 y1 .
If one multiplies y1 on both sides of v one gets
y0 v y0 = v0 (y0 y0 y0) + v1 (y0 y1 y0) .
But by the two rules of Clifford multiplication above it
follows that
y0 y0 y0 = -y0
and
y0 y1 y0 = - y0 y0 y1 = y1 .
So that finally
y0 v y0 = -v0 y0 + v1 y1,
which is v reflected at the "hyperplane" orthogonal to
y0, i.e. reflected at the line spanned by y1.
Now for rotations: Every rotation can be constructed
from two successive reflections. In particular, reflecting
a vector in succession at two intersecting hyperplanes
is tantamount to rotating this vector by an angle that is
exactly *twice* the angle between the hyperplanes.
(I have emphasized "twice" here, because it is exactly
this innocent and easily visualized factor of 2 that can be
traced back to be at the root of the not quite so easily
visualizable sign-change behavior of spinors under rotation
by 2pi, see below.)
Another simple example will easily illustrate this: Consider
again the Euclidean plane with orthogonal unit vectors y0
and y1 and the arbitrary vector v = v0 y0 + v1 y1.
Now watch how v behaves when it is reflected first
at y1 and then at y0:
y1 (y0 v y0) y1
= y1 (y0 (v0 y0 + v1 y1) y0) y1
= y1 (-v0 y0 + v1 y1) y1
= - v0 y0 - v1 y1
= -v
Hence, by these two reflections v is rotated by pi around the
origin, which is exactly twice the angle of pi/2 between y0 and y1.
If one instead wants to rotate v by an arbitrary angle, then one
has to exponentiate the action of the discrete pi rotation above:
v' = v rotated by the angle phi
= exp(phi/2 y1 y0) v exp(phi/2 y0 y1),
= R v R~
where the factor 1/2 is inserted to compensate the above
mentioned factor 2 and the tilde indicates that y1 y0 has
to be "reversed" to give y0 y1.
As I have mentioned in another recent post, the object
R, sometimes called a "rotor", which is "one
half" (the "left half" in this case) of a rotation operator, is
essentially a spinor. Why? Because under rotation it behaves
"just like" a vector, except that it picks up a sign when
rotated by 2pi. To see this, rotate the v' once more, this
time by an angle psi:
v'' = v' rotated by an angle psi
= exp(psi/2 y1 y0) v' exp(psi/2 y0 y1)
= exp(phi/2 y1 y0) R v R~ exp(psi/2 y0 y1)
= (exp(phi/2 y1 y0) R) v (exp(phi/2 y1 y0) R)~
Hence, while the *vector* v is rotated around by applying
rotation operators on *both* sides, the *spinor* R
is transformed under rotation by applying
exp(phi/2 y1 y0) only from the left. But because of the
factor 1/2 in the exponent of the latter (which was
introduced in order to cancel the factor of 2 that arose
because rotations were constructed from reflections)
the operator exp(phi/2 y1 y0) is equal to -1 for
a full rotation of phi = 2pi, not equal to 1. This does not
affect the transformation of the vector v, since it
is multiplied under 2pi rotations by -1 from both sides,
giving a total of 1. But it means that under a 2pi rotation
the object R transforms as
R -> -R .
Hence R is a spinor, qed. ... Well, ok, this is not really
precise and thus not really correct. I'll make it a little
more precise and thus a little more correct, but only so much
that I can discuss the physics that I actually set out to
discuss:
The somewhat more precise version, and I'll make a long story
short, is that which involves "minimal left ideals": The simple
point is that one wants the "algebraic spinor" R (as it is called)
to
only notice what is going on on its left side. Hence, loosely
speaking, one projects its right side onto some subspace and
forgets about it.
I'll explain that again in terms of an example: To get closer to
the physics let the Euclidean plane from the previous examples
now be equipped with the lorentzian metric (1,-1), i.e.
y0 y0 = 1, y1 y1 = -1 .
In order to project something we need a projector constructed
from y0 and y1 and I'll let that projector be
P = (1 + y0)/2
which implies
y0 P = P
P P = P ,
so that P is recognized as the projector onto the +1 eigenspace
of y0.
The point is that it is quite correct to call the rotator R from
above a spinor when it is applied on the projector P:
RP is a spinor.
But remember, the point is merely that RP is transformed
by acting on it from the left, the P is essentially only there
to make precise the idea that we do not care about what
happens on the right hand side.
But since now R is applied to a projector, several different
Rs give the same result RP. A simple inspection shows that the
space of all possible RPs is spanned by only two elements,
which can for example be chosen as follows:
r = (1 + y0 y1)/2 P
l = (1- y0 y1)/2 P
The reason why I call these elements "r" and "l" will arise
shortly.
Observe that from general considerations we know that
spinors in 2n dimensions have 2^n components. Incidentally,
in 1+1 dimensions considered here, 2^1 = 2 are two
components of a spinor, which complies with the two
components r and l found above. Hence, a spinor valued
wave function psi in 1+1 dimensions can be written as
psi = rho r + lambda l
where rho and lambda are complex coefficients that
give a weight (or rather "amplitude") to the components
r and l.
This finally puts us in position to do some interesting
physics: The Dirac equation in 1+1 Minkowski space reads
( y0 d0 + y1 d1 ) psi = i m psi .
(Where d0 and d1 are partial derivatives in the y0 and y1
direction (i.e. time and space direction, respectively)
acting on the components lambda and rho.)
As is often done, I multiply this equation by y0 on
both sides for convenience
( d0 + y0 y1 d1 )psi = i m y0 psi
and reorder a little to make the time evolution transparent:
d0 psi = - y0 y1 d1 psi + im y0 psi .
This is nice, because with the above discussion of reflections
and rotations we can now do what is most fun in physics: We
can point at each term in an equation and give it a fancy
physical interpretation. First observe that
y0 y1 r = r
y0 y1 l = -l
y0 r = l
y0 l = r .
y0 y1 is, by the above discussion, the generator of rotations
in the space-time plane, that is: of Lorenz boosts in the
y1 direction. y0 y1 r = r says that r is invariant under
such boosts. But this must mean that r is a lightlike object
pointing in the y1 direction. Similarly, l is a lightlike
object pointing in the opposite direction. This is consistent
with y0 r/l = l/r, which says, that under time reflection
r and l interchange roles.
One can now read off the little story that the above equation
is
telling:
In each time step (d0 psi) three different things happen:
The r component of psi is translated at the speed of light
a little in the -y1 direction (-d1 rho) while the l component is
translated at the speed of light a little in the +y1 direction
(d1 lambda). Furthermore, part of the r component is
switched to the l component (i m y0 r = i m l) and
vice versa (im y0 l = i m r).
This "story of the zig-zagging electron" has first been discovered
by Feynman for 1+1 dimensions, who called it the "chessboard
model", and has been generalized to 1+3 dimensions by
Ord, a few years ago (but in ugly matrix notation not in nice
Clifford algebra). Remarkably, this story contains no mention
of
spin, only of reflections, and yet it describes "spinning"
electrons.
I hope the above derivation makes this at first puzzling connection
a little transparent.
Although I have not done it yet, by means of the Clifford
algebra approach above it should be easy to generalize
the checkerboard model to higher dimensions and
arbitrary complex interaction terms. There are two important
points here:
The y0 yi which multiply the partial derivatives of the
free Dirac operator have eigenvalues +/- 1. Hence
the free Dirac operator explicitly says: Translate all +1
eigenstates in one direction and all -1 eigenstates in the
other (at the speed of light!), and then swap some
eigenstates (in higher dimensions).
The interaction terms, like i m y0 psi, in general
swap the "directional" eigenstates. By the above discussion
one can easily see how each term of a vector potential
A = Ai yi
affects the checkerboard model (i.e. what exactly is does
to a zig-zagging electron at each discrete instance of time.)
This is fun if interaction is gravity: Writing the Dirac equation
for curved space-time
(yi di + 1/4 yi omega_ijk yj yk )psi = i m psi
(where omega is the spin connection), one sees that the
connection influences the weights of the left and right
going components of the electron path. For example, the
simple case of 1+1 dimensions has the general affine
connection
omega
= yi omega_ijk yj yk
= omega_001 y1 + omega_101 y0
(because of antisymmetry in jk).
Multiplying by y0 as above gives the interaction
omega_001 y0 y1 + omega_101 .
which increases the left going amplitude by
omega_001 - omega_101
and the right going amplitude by
-omega_001 - omega_101 .
More information on the checkerboard model can be
found at
http://www.cbloom.com/physics/2d_dirac.html#new
(which apparently originates from old postings here on
s.p.r) and references given there.
You may have noticed that the interactions I discussed in my above
discussion, the vector potential and gravity, did non introduce kinks
into the
electron path, but merely influenced the amplitudes of the left and
right
moving components, thereby luring the electron in one direction or
another. But wouldn't one expect that a potential in the checkerboard
model
instead makes itself felt by directly influencing the probability for
reflections of the electron's velocity? Playing ping-pong with the
particle?
I would, and here is how:
Instead of a vector-potential, I now introduce a super-potential (as
in the Witten model of supersymmetric quantum mechanics). I can
provide the details if desired, but the result is that the 1+1 dimensional
massless Dirac operator then looks like
D = y0 d0 + y1 d1 + Y0 W,0 + Y1 W,1
where
y0 and y1 are the generators of Cl(1,1) as before
d0 and d1 are the partial derivatives
Y0 and Y1 are generators of *another* copy of Cl(1,1), anticommuting
with y0 and y1
W,0 and W,1 are the d0 and d1 derivatives of a scalar function W.
Because of the new generators, Y0 and Y1, this Dirac operator
does not act on a single spin bundle S any more, but on the product
of two spin bundles SxS*. This can be represented by introducing
a "vacuum" element |0> and the rule
y0 |0> = Y0 |0>
y1 |0> = Y1 |0> .
When the spinor states of my previous post
r/l = (1 +/- y0 y1) (1 + y0)/4
are applied to this vacuum state,
r/l -> (1 +/- y0 y1) (1 + y0) /4 |0>
they also span a representation space for the action of Y0 and Y1.
One finds:
Y0 (1 +/- y0 y1) (1 + y0) /4 |0>
= (1 +/- y0 y1) (1 - y0) /4 Y0 |0>
= (1 +/- y0 y1) (1 - y0) /4 y0 |0>
= - (1 +/- y0 y1) (1 - y0) /4 |0>
and
Y1 (1 +/- y0 y1) (1 + y0) /4 |0>
= (1 +/- y0 y1) (1 - y0) /4 Y1 |0>
= (1 +/- y0 y1) (1 - y0) /4 y1 |0>
= (y1 -/+ y0) (1 + y0) /4 |0>
Note that Y0 switches the projector (1+y0)/2 to its complement, hence
it annihilates the (1+y0)/2 particle that has been considered so far
and
instead creates a new kind of particle living in the (1-y0)/2 representation!
We'll also need the actions of y0 Y0 and y0 Y1. These are now readily
seen
to be
y0 Y0 (1 +/- y0 y1) (1 + y0) /4 |0>
= - y0 (1 +/- y0 y1) (1 - y0) /4 |0>
= (1 -/+ y0 y1) (1- y0)/4 |0>
and
y0 Y1 (1 +/- y0 y1) (1 + y0) /4 |0>
= y0 (y1 -/+ y0) (1 + y0) /4 |0>
= -(1 -/+ y0 y1) (1 + y0) /4 |0> .
Hence both, y0 Y0 and y1 Y1, reverse the direction of the particle.
The
former also switches the representation while the latter introduces
a sign.
Ok, so this is what the terms in the above Dirac operator do. As before,
I now write out the Hamiltonian, i.e. the generator of time translations:
D psi = 0
<=> (y0 d0 + y1 d1 + Y0 W,0 + Y1 W,1) psi = 0
<=> d0 psi = ( -y0 y1 d1 - y0 Y0 W,0 - y0 Y1 W,1) psi
The operators on the right each do something to the "electron" in
a small unit of time, namely:
y0 y1 d1 translates the r component in the -y1 direction and the
l component
in the +y1 direction a little bit, at the speed of light (as
before)
y0 Y0 W,0 annihilates the particle and creates a corresponding
particle
going in the opposite direction with amplitude W,0
y0 Y1 reverses the direction of the particle with amplitude W,1
Physical intuition is rather delighted:
There is a particle moving around at the speed of light, stochastically
reversing direction as on Feynman's checkerboard, but with a
probability
that varies with the strength of the spatial potential W,1: The
higher the
potential, the more the particle's path wiggles and the slower the
particle
propagates (in the mean). Also, because the probability to return when
heading
in a direction of lower potential decreases, the particle will (in
the mean)
travel
in the direction opposite to the potential's gradient. (There
is a subtlety
here,
because the proper potential is the square of W,1, but the qualitative
idea is
correct.) If, in addition, the potential depends on time, i.e. when
W,0 =/= 0 ,
so that the energy content of the system is not constant, the particle
may
"decay" into its partner representation and back again.
This is, however, not the story of the electron, anymore.
--------------
> or I'm
> asking for something impossible like a human brain really watching
4
> dimensions graphically?
No. Of course, you can have Dirac equations in all dimensions and
on all kinds of rather exotic geometries, and these may be difficult
to handle by themselves, but for ordinary setups there is an ordinary
visualization of the Dirac equation. It goes like this (see also [3]):
Let's first concentrate on the homogeneous case, e.g. when a
massless particle (like a neutrino) is described. The the Dirac
equation reads
D psi = 0
where
D = ym nabla_m
is the Dirac operator.
Here ym are the generators of the Clifford
algebra
{ym,yn} = -2 gmn
and gmn are the components of the metric on your manifold. In flat
Minkowski space this usually looks like
{ym,yn} = 2 eta_mn
where
eta = diag(+1,-1,-1,-1) .
All Dirac operators always look this way (in local coordinates).
nabla_m is a covariant derivative
nabla_m = @m + omega_m
where @m is the partial derivative by x^m and omega_m are the
components of a connection on some fiber bundle. Really, nabla_m
are components of the connection 1-form
nabla = dx^m (x) nabla_m
that, when applied to a tangent vector v = v^n @_n, returns the
operator v^n nabla_n. This, in turn, is a differential operator that,
when applied to any field, say S, returns the first order
approximation of the change of S when parallely transported along
v.
This is the reason why Quillen says that Dirac operators are
"quantized connections". All Dirac operators arise by replacing
in the connection 1-form
dx^m (x) nabla_m
the differential form dx^m by the Clifford generator ym. This
replacement is sometimes called "quantization", because
{dx^m, dx^n} = 0
and
{ym, yn} = 2 eta^mn
differ like commutators of classical phase space functions and
quantum operators. But "quantized connections" have nothing to
do with quantization in the sense of quantum physics, it's just a
formal analogy. Dirac operators have their place in classical physics,
too.
So
D = ym nabla_m
is a *very* general construction. Physics enters by specifying a
fiber bundle and a specific connection on that bundle so that
nabla gets a definite meaning.
Unfortunately, I am running somewhat out of steam and also out of
time. I'll leave it at that for the moment, which is, admittedly, a
very
incomplete state of things. But I'd be glad to discuss things further.
> Screw operators, screw differential equation, hail the path integral!
A Japanese proverb says (or so I have been told):
When you do not understand, mountains are mountains.
When you begin to understand, mountains are no longer mountains.
When you have understood, mountains will again be mountains.
----- references:
[1]
@book{
BerlineGetzlerVergne:1992,
author = {N. Berline and E. Getzler and M. Vergne},
title = {Heat Kernels and {D}irac Operators},
publisher = {Springer},
year = {1992}
}
@book{
LawsonMichelsohn:1989,
author = {H. Lawson and M. Michelsohn},
title = {Spin Geometry},
publisher = {{P}rinceton {U}niversity Press},
year = {1989}
}
@book{
BennTucker:1987,
author = {I. Benn and R. Tucker},
title = {An introduction to Spinors and Geometry},
publisher = {{A}dam {H}ilger},
year = {1987}
}
[1b]
@book{
Frankel:1997,
author = {T. Frankel},
title = {The Geometry of Physics},
publisher = {Cambridge University Press},
year = {1997}
}
[2]
@article{
RoepsdorffVehns:1999,
author = {G. Roepsdorff and Ch. Vehns},
title = {Generalized {D}irac Operators and Superconnections},
journal = {math-ph/9911006},
year = {1999}
}
[3]
@article{
RodriguesVazPavsic:1996,
author = {W. Rodrigues and J. Vaz and M. Pavsic},
title = {The {C}lifford bundle and the dynamics of the superparticle},
journal = {{B}anach Center Publications},
year = {1991},
volume = {37},
pages = {295-315}
}
[4]
@article{
RodriguesDeSouzaVazLounesto:1996,
author = {{W. Rodrigues, Jr.} and {Q. De Souza} and {J. Vaz,
Jr.} and P.
Lounesto},
title = {{D}irac-{H}estenes spinor fields in {R}iemann-{C}artan
spacetime},
journal = {\tt hep-th/9607073},
year = {1996}
}
[5]
@article{
Varilly:1997,
author = {J. Varilly},
title = {An introduction to Noncommutative Geometry},
journal = {\tt physics/9709045},
year = {1997}
}
@book{
Connes:1994,
author = {A: Connes},
title = {Noncommutative Geometry},
publisher = {Academic Press},
year = {1994}
}
[6]
@book{
Hestenes:1966,
author = {D. Hestenes},
title = {Space-time algebra},
publisher = {Gordon and Breach science publishers},
year = {1966}
}
[7]
@book{
HestenesSobczyk:1984,
author = {D. Hestenes and G. Sobczyk},
title = {{C}lifford Algebra to {G}eometric {C}alculus},
publisher = {D. Reidel Publishing Company},
year = {1984}
}
[8]
@book{
Hestenes:2002,
author = {D. Hestenes},
title = {Spacetime Calculus},
publisher = {(in preparation, draft available at {\tt
modelingts.la.asu.edu/HTML/STC.html})},
year = {2002}
}
--
Urs.Schreiber@uni-essen.de
From: Urs Schreiber (Urs.Schreiber@uni-essen.de)
Subject: Re: Dirac for Dunces
Newsgroups: sci.physics.research
Date: 2002-01-12 12:34:01 PST
"Urs Schreiber" <Urs.Schreiber@uni-essen.de> schrieb im Newsbeitrag
news:a1euhh$js0$1@rs04.hrz.uni-essen.de...
[...]
> Mathematically, there is a
> close relation between spinors and reflections. The fascinating
> implication for physics is that discrete checkerboard models,
> where massive particles move around at the speed of light but
> reverse direction (reflection!) every now and then, actually
> describe *spinning* particles!
[...]
> Actually, I am also too tired to continue writing this post,
> so I'll stop here.
Now that I'm fresh again, I can't refrain from adding one more interesting
detail. You may have noticed that the interactions I discussed in my
above
post, the vector potential and gravity, did non introduce kinks into
the
electron path, but merely influenced the amplitudes of the left and
right
moving components, thereby luring the electron in one direction or
another. But wouldn't one expect that a potential in the checkerboard
model
instead makes itself felt by directly influencing the probability for
reflections of the electron's velocity? Playing ping-pong with the
particle?
I would, and here is how:
Instead of a vector-potential, I now introduce a super-potential (as
in the Witten model of supersymmetric quantum mechanics). I can
provide the details if desired, but the result is that the 1+1 dimensional
massless Dirac operator then looks like
D = y0 d0 + y1 d1 + Y0 W,0 + Y1 W,1
where
y0 and y1 are the generators of Cl(1,1) as before
d0 and d1 are the partial derivatives
Y0 and Y1 are generators of *another* copy of Cl(1,1), anticommuting
with y0 and y1
W,0 and W,1 are the d0 and d1 derivatives of a scalar function W.
Because of the new generators, Y0 and Y1, this Dirac operator
does not act on a single spin bundle S any more, but on the product
of two spin bundles SxS*. This can be represented by introducing
a "vacuum" element |0> and the rule
y0 |0> = Y0 |0>
y1 |0> = Y1 |0> .
When the spinor states of my previous post
r/l = (1 +/- y0 y1) (1 + y0)/4
are applied to this vacuum state,
r/l -> (1 +/- y0 y1) (1 + y0) /4 |0>
they also span a representation space for the action of Y0 and Y1.
One finds:
Y0 (1 +/- y0 y1) (1 + y0) /4 |0>
= (1 +/- y0 y1) (1 - y0) /4 Y0 |0>
= (1 +/- y0 y1) (1 - y0) /4 y0 |0>
= - (1 +/- y0 y1) (1 - y0) /4 |0>
and
Y1 (1 +/- y0 y1) (1 + y0) /4 |0>
= (1 +/- y0 y1) (1 - y0) /4 Y1 |0>
= (1 +/- y0 y1) (1 - y0) /4 y1 |0>
= (y1 -/+ y0) (1 + y0) /4 |0>
Note that Y0 switches the projector (1+y0)/2 to its complement, hence
it annihilates the (1+y0)/2 particle that has been considered so far
and
instead creates a new kind of particle living in the (1-y0)/2 representation!
We'll also need the actions of y0 Y0 and y0 Y1. These are now readily
seen
to be
y0 Y0 (1 +/- y0 y1) (1 + y0) /4 |0>
= - y0 (1 +/- y0 y1) (1 - y0) /4 |0>
= (1 -/+ y0 y1) (1- y0)/4 |0>
and
y0 Y1 (1 +/- y0 y1) (1 + y0) /4 |0>
= y0 (y1 -/+ y0) (1 + y0) /4 |0>
= -(1 -/+ y0 y1) (1 + y0) /4 |0> .
Hence both, y0 Y0 and y1 Y1, reverse the direction of the particle.
The
former also switches the representation while the latter introduces
a sign.
Ok, so this is what the terms in the above Dirac operator do. As before,
I now write out the Hamiltonian, i.e. the generator of time translations:
D psi = 0
<=> (y0 d0 + y1 d1 + Y0 W,0 + Y1 W,1) psi = 0
<=> d0 psi = ( -y0 y1 d1 - y0 Y0 W,0 - y0 Y1 W,1) psi
The operators on the right each do something to the "electron" in
a small unit of time, namely:
y0 y1 d1 translates the r component in the -y1 direction and the
l component
in the +y1 direction a little bit, at the speed of light (as
before)
y0 Y0 W,0 annihilates the particle and creates a corresponding
particle
going in the opposite direction with amplitude W,0
y0 Y1 reverses the direction of the particle with amplitude W,1
Physical intuition is rather delighted:
There is a particle moving around at the speed of light, stochastically
reversing direction as on Feynman's checkerboard, but with a
probability
that varies with the strength of the spatial potential W,1: The
higher the
potential, the more the particle's path wiggles and the slower the
particle
propagates (in the mean). Also, because the probability to return when
heading
in a direction of lower potential decreases, the particle will (in
the mean)
travel
in the direction opposite to the potential's gradient. (There
is a subtlety
here,
because the proper potential is the square of W,1, but the qualitative
idea is
correct.) If, in addition, the potential depends on time, i.e. when
W,0 =/= 0 ,
so that the energy content of the system is not constant, the particle
may
"decay" into its partner representation and back again.
This is, however, not the story of the electron, anymore.
Urs.Schreiber@uni-essen.de
Viewing
message <a8c9d3$up0$1@rs04.hrz.uni-essen.de>
From: Urs Schreiber (Urs.Schreiber@uni-essen.de)
Subject: Re: Physically understanding the Dirac equation and
4D.
Newsgroups: sci.physics.research
Date: 2002-04-03 19:44:31 PST
"George Raetz" <bestwork.1@pcisys.net> schrieb im Newsbeitrag
news:uafo636qlr9sc7@corp.supernews.com...
[...]
> Concerning the generalization of the Feynman Checkerboard to 1+3 dimensions:
> In 1+1 dimensions, one follows the electron along a two dimensional
lattice,
> tallying each direction independently and weighing the paths
with the
> operator connecting the lattice directions number. (one thinks of
the factor
> i as
> connecting two lattice directions (1+i)i=i-1). That leads to
> a 4 component description.which fits nicely with the 1+1 dimensional
Dirac
> equation where the wave function is a complex 2-vector--again 4 numbers.
> To generalize to 3+1 dimensions it would seem you must keep
track of four
> directions and weigh the paths with rotation operators that connect
the
> lattice directions; and those operators could be complex quaternions
or 3d
> Clifford Algebra operators--8 component objects. In all a total of
4x8=32
> quantities.
> Far in excess of the 8 quantities of the complex Dirac 4-vector
wave
> function. It
> can be shown however, that the Dirac equation is a subset of this
process:
> http://www.pcisys.net/~bestwork.1/QRW/RandomWalk.htm
Thank you very much for this reference! I know of Tony Smith's very
interesting construction, but I am having problems seeing how it reproduces
the correct Dirac equation.This may very well be due to my own
insufficient understanding, though I have not found any proof
on his
site, either (is there one that I have missed?). A while ago we were
discussing
these matters here on s.p.r. I have appended a relevant post of mine
below. My would-be approach to generalizing the checkerboard model
to
1+3 dimensions was based on the idea that instead of *guessing* the
correct
sum-over-paths and checking if it reproduces the Dirac equation one
could try to *derive* the sum-over-paths from the Dirac equation by
looking
at the discrete approximation to the time evolution induced by it.
This
does work for the 1+1 dim case, but becomes difficult in 1+3 dimensions
since here there are no eigenstates of lightlike propagation in all
three
spatial directions (see appended post below for details).
Now you say you have shown that Tony Smith's construction reproduces
the correct Dirac equation as a subset of its possible solutions. That's
very interesting. As far as I know there is no preprint or journal-paper
doing
exactly that. This paper by Ord and McKeon
@article{
OrdMcKeon:1992,
author = {G. Ord and D. McKeon},
title = {On the {D}irac equation in $3+1$ dimensions},
journal = {Ann. Phys.},
year = {1992},
volume = {222},
pages = {244-253}
},
which you probably know of, considers only the special case of plane
waves, where |psi(x,y,z)> = |psi(x)> so that the 1+3 dim Dirac equation
reduces to a 1+1 dim like equation because of
y^m d_m |psi(x)> = y^x d_x |psi(x)>.
(That's my paraphrase of the central point in this paper.)
Even though, as the authors point out, every solutions can be obtained
by superposing plane waves, this restriction seems unsatisfactory.
Also, I gather that Louis Kauffman is currently working on the 1+3 dim
problem again. So it seems this is still considered an open problem,
officially.
Hence I have taken a look at the material you present on your site,
especially
that titled "The HyperDiamond Random Walk", found at
http://www.pcisys.net/~bestwork.1/QRW/the_flow_quaternions.htm ,
which is mostly new to me. This is very illuminating. But I'd kindly
ask
you to help me with a couple of problems I have with some of your claims.
I'll briefly summarize:
First you write down the difference equations of the HyperDiamond process
proposed by Tony Smith. Then you rewrite that as differential equations,
obtaining:
(1) d_k K_k = s_x K_l + s_y K_m + s_z K_n , (no sum over k)
where (k,l,m,n) runs over four permutations of the four directions
accessible from a vertex in the HyperDiamond lattice (termed r,g,b,y
in
your text).
To transform this into a form that resembles a Dirac equation, you make
the following ansatz for the amplitudes K:
(2) K_m = q_m Q
where Q is a single component quaternionic function and q_m are the
quaternionic versions of the vectors pointing along the HyperDiamond
edges at one vertex.
The main result, which is new to me and which I find very interesting,
is
that with this ansatz for K the *sum* of all four equations above reads
(up to a constant factor)
(3) q_m d_m Q = i E Q (sum over m) .
The left hand side is the perfectly covariant Dirac operator, which
you
write in the q_m basis. The right hand side looks almost like the right
hand
side of the Dirac equation, except for that quaternion number E, which
you
find to be the unit vector in the spatial (1,1,1) direction:
(4) E = (e_1 + e_2 + e_3)/sqrt(3) , E^2 = 1
While the left hand side is obtained directly by virtue of your ansatz,
it is
interesting that the right hand side has this simple form. It is suggestive
now
to play around with this until it looks completely like the Dirac equation.
But one thing worrying me at this point is that we started out with
*four*
equations (1) and are now looking merely at their *sum* (3). It is
non-trivial
to me that one does not loose information this way. I mean, suppose
you
solve (3), which gives you Q and thus K_m via (2). Will these
K_m
necessarily solve all four equations of (1) separately? After
staring at your
equations for a while I do not see how this should work. But
there is a lot
of symmetry in the problem, so I can imagine it may be true. Can you
show it?
Or else, can we get another pretty right-hand-side of (3) by, say,
adding the
first two and subtracting the last two equations of (1)?
This is related to another question: Can your ansatz (2) be translated
back
into a sum-over-paths prescription for Q? I mean, suppose I do the
sum
over paths following Tony Smith and end up with four quaternionic
functions K_m. What if these do not admit a Q such that (2) holds?
I see no simple reason why they should, in general. Do you?
But let's assume (4) and see how this may become the correct Dirac
equation:
With an eye on the usual chiral splitting, you suggest to use two Q's,
Q1, Q2, such that
E Q1 = Q2 .
You derive two coupled equations
q_m d_m Q1 = i Q2
q'_m d_m Q2 = i Q1 ,
where q' are the time reversed diamond edges. This looks very similar
to
the usual chiral splitting of the Dirac equation, but I do not quite
see that
it is really equivalent. I may have to study this closer. In equation
(26) you
say that it is equivalent for the choice
Q1 = {{psi1, psi3},{psi2, psi 4}}
Q2 = {{psi3, psi1},{psi4, psi 2}} .
This puzzles me. These Qi do not seem to satisfy E Q1 = Q2. Maybe I
am
confused, but I believe to have understood that in matrix representation
we
have
sqrt(3)E = sum of Pauli matrices = {{1,1+i},{1-i,-1}} .
Right? But this implies E Q1 =!= Q2 for general psi.
Please note that I am only trying to understand your proposal. If I
am asking stupid questions then please straighten me out!
--
Urs.Schreiber@uni-essen.de
-------------- copy of recent post relevant to this discussion: ----------
"Danny Ross Lunsford" <antimatter33@worldnet.att.net> schrieb im
Newsbeitrag
news:Poh98.5432$9%6.1085077@bgtnsc04-news.ops.worldnet.att.net...
>
> "Urs Schreiber" <Urs.Schreiber@uni-essen.de> schrieb:
>
> > As far as I know, nobody has yet given a *general* checkerboard
> > summing prescription for the Dirac propagator in 1+(d>1) dimensions.
>
> Look on Tony Smith's page here:
>
> http://www.innerx.net/personal/tsmith/FynCkb.html
Yes, I know this page and the corresponding lanl paper
http://xxx.lanl.gov/PS_cache/quant-ph/pdf/9503/9503015.pdf
and I know that Tony Smith does give a generalization of
Feynman's summing prescription from 1+1 to 1+3 dimensions.
But I have to say that I fail to see that this generalization
reproduces the Dirac propagator in 1+3 dimensions, and that I
did not find any proof that it does. Actually, I seem to have
convinced myself that it does not, but I may of course be quite
wrong. I therefore take this opportunity to state my understanding
of these matters.
First, I very briefly summarize (my understanding of) Tony Smith's
construction:
The starting point is the observation that the left |-> and right |+>
going states of the 1+1 dim checkerboard model can be labeled by
complex numbers
|-> ---> (1 + i)
|+> ---> (1 - i)
(up to a factor)
so that multiplication by the negative imaginary unit swaps
components:
(-i) (1 + i)/2 = (1 - i)/2
(-i) (1 - i)/2 = (1 + i)/2 .
Since the path-sum of the 1+1 dim model reads
phi =
sum over all possible paths of (-i eps m)^(number of bends of path)
=
sum over all possible paths of
product over all steps of one path of
-i eps m (if change of direction after this step
generated by i)
1 (otherwise)
this makes it look very natural to identify the imaginary unit
appearing in the sum over paths with the "generator" of kinks in
the path. To generalize this to higher dimensions, more square
roots of -1 are added, which gives the quaternion algebra in 1+3
dimensions. The two states |+> and |-> from above, which were
identified with complex numbers, are now generalized to four states
identified with the following quaternions (which can be identified
with vectors in M^4 indicating the direction in which a given path
is heading at one instant of time):
(1 + i + j + k)
(1 + i - j - k)
(1 - i + j - k)
(1 - i - j + k) ,
which again constitute a (minimal) left ideal of the algebra (meaning
that applying i,j, or k from the left on any linear combination of
these
four states gives another linear combination of these four states).
Hence, now i,j,k are considered as "generators" of kinks in three
spatial dimensions and the above summing prescription naturally
generalizes to
phi =
sum over all possible paths of
product over all steps of one path of
-i eps m (if change of direction after this step
generated by i)
-j eps m (if change of direction after this step
generated by j)
-k eps m (if change of direction after this step
generated by k)
1 (otherwise)
The physical amplitude is taken to be
A * e^(i alpha)
where A is the norm of phi and alpha the angle it makes with the
x0 axis.
As I said, this is merely my paraphrase of Tony Smith's proposal as
I understand it.
I fully appreciate that the above construction is a nice (very "natural")
generalization of the summing prescription of the 1+1 dim checkerboard
model. But if it is to describe real fermions propagating in physical
spacetime, this generalized path-sum has to reproduce the propagator
obtained from the Dirac equation in 1+3 dimensions, which we know
to correctly describe these fermions. Does it do that?
I currently doubt that it does and I will sketch what leads me to this
conclusion. But note that I have no proof, just some plausibility
considerations. (On the other hand, a simple computer simulation
would clarify the matter. Has such a simulation be done? Or are there
even exact solutions to Tony Smith's path-sum? In such a case the
rest of this post is rather obsolete.)
First, it is worth looking at how the 1+1 dim checkerboard model
actually follows from the 1+1 dim Dirac equation. I have said this
before in this thread, but I will do so again for completeness:
In 1+1 dim the free massive spinning particle is described by
D|psi> = 0
<=>
(y0 d0 + y1 d1 + im) |psi> = 0
<=>
-d0 |psi> = y0 y1 d1 + i y0 m |psi>
<=>
d0 |psi> = H/ih |psi>
H/ih = (- y0 y1 d1 - i y0 m)
(yi being basis elements of Cl(1,1), y0^2 = 1, y1^2 = -1).
In the context of discrete approximations to continuous dynamics
the partial derivatives di can be replaced by difference operators
[c d0] |psi(x0, x1)> = |psi(x0+c,x1)> - |psi(x0, x1)>
[c d1] |psi(x0, x1)> = |psi(x0,x1+c)> - |psi(x0, x1)>
(where c is some real factor and the discrete step size is set to
unity, epsilon=1, for convenience of notation).
The spinor wave function |psi> is sitting in some minimal left ideal
of the Clifford algebra, for example that generated by (1 + y0)/2:
|psi> = psi+ |+> + psi- |->
|+> = (1+ y0 y1)/2 (1 + y0)/2
|-> = (1- y0 y1)/2 (1 + y0)/2 .
Now cosider an initial condition with a particle sitting at the
origin
x0=0, x1=0 with only a |+> component:
|phi(x0=0)> = delta(x1) |+> .
The formal solution to D|phi> = 0 is
|phi(x0,x1)> = exp(x0 H/ih) delta(x1) |+> .
This says that after a discrete infinitesimal time step on finds:
|phi(x0=1,x1)>
= (1 + H/ih) delta(x1) |+>
= delta(x1) |+> + [-1 d1] delta(x1) |+> - im delta(x1)|->
= delta(x1) |+> + delta(x1-1) |+> - delta(x1) |+> - im delta(x1)|->
= delta(x1 - 1) |+> - im delta(x1) |-> ,
that is, the superposition of a particle at x1=1 in the |+> state and
a particle still at x1=0 but in the |-> state.
The next step gives
|phi(x0=2,x1)>
= delta(x1 - 2) |+> -im delta(x1 + 1)|->
- im delta(x1 - 1) |-> - delta(x1) |-> ,
and so on. This clearly shows that if one wants to know the amplitude
at some point (x0,x1) one can either continue this kind of time
propagation, or one can take the shortcut, recognize the pattern and
apply
Feynman's summation rule. The relation between both approaches is
exactly that expressed by the equation
(1+a)^N = sum_n binomial(N,n) a^n .
Now for 1+3 dimensions:
The Dirac equation reads:
D|psi> = 0
<=>
-d0 |psi> = (y0 y1 d1 + y0 y2 d2 + y0 y3 d3 + im y0) |psi>
<=>
d0|psi> = H/ih |psi>
H = y0 y1 d1 + y0 y2 d2 + y0 y3 d3 + im y0
and one can identify the quaternions as usual:
i <-> y0 y1
j <-> y0 y2
k <-> y0 y3 .
The first thing to observe is that one cannot construct a state that
is an eigenstate to more than one of these, since they anticommute.
In particular there is no state
delta(x1, x2, x3) | >
which evolves into
(1 + H_free) delta(x1, x2, x3) | >
= (1 + y0 y1 d1 + y0 y2 d2 + y0 y3 d3) delta(x1, x2, x3) | >
= delta(x1+1, x2+1, x3+1) | >
in a discrete time step. But this is what one would expect from
a path labeled (at one instant) by the quaternion
(1 + i + j + k)
in the quaternion path-sum discussed above.
I say "expect", because I realize that the exact relation of the
quaternion path sum to a local time evolution may be subtle due
to the final "projection" rule that assigns a complex number to
the quaternionic outcome of the path-sum.
Why don't we go the other way round and try to read off the
summation rule from the time generator? Let's see:
To represent spinors, choose the ideal generated by
(1 + y0)/2 .
Two component Weyl spinors are obtained by acting from the
left on
(1 + y0 y3)/2 (1 + y0)/2
or
(1 - y0 y3)/2 (1 + y0)/2 .
The two component Weyl spinors can be written
|eta>
= (eta1 + eta2 y0 y2) (1 + y0 y3)/2 (1 + y0)/2
= (eta+ (1 + y0 y2)/2 + eta- (1 - y0 y2)/2)
(1 + y0 y3)/2 (1 + y0)/2
|xi>
= (xi1 + xi2 y0 y2) (1 + y0 y3)/2 (1 + y0)/2
= (xi+ (1 + y0 y2)/2 + xi- (1 - y0 y2)/2)
(1 + y0 y3)/2 (1 + y0)/2
The mass term of the time evolution generator
i m y0
switches representations:
y0 (1 + y0 y3)/2 (1 + y0) /2 = (1 - y0 y3)/2 (1 + y0) .
Ok, what happens to a particle at the origin in the state
delta(x1,x2,x3) (1 + y0 y2)/2 (1 + y0 y3)/2 (1 + y0)/2
i.e. a particle with
eta+ = delta(x1,x2,x3)
and all other components zero.
One finds:
(1 + H) delta(x1,x2,x3) (1 + y0 y2)/2 (1 + y0 y3)/2 (1 +
y0)/2
=
delta(x1,x2-1,x3) (1 + y0 y2)/2 (1 + y0 y3)/2 (1 + y0)/2
+
(terms that I cannot bring into any pretty order right now).
In order to find a suitable summation prescription one has to
have a good idea how these
(terms that I cannot bring into any pretty order right now)
could be collected to yield some pattern.
Obviously, I do not have well established ideas on these matters,
just some half-baked ones.
> I'm skeptical of results in low dimensions as being of physical value
as
> models, even if they are mathematically interesting. The best example
is
> curvature, where the full thing only appears for N>3. (For N=3 the
Weyl
> tensor vanishes.)
That's true. And then, there are systems that actually are 1+(n<3)
dimensional, at least to a good approximation. One thing to note is
that we need not consider dynamics in spacetime but may instead be
looking at dynamics in configuration spaces of arbitrary dimension.
(Of course, one is really *always* considering dynamics in configuration
space. It just happens that for a single point particle physical spacetime
and configuration space coincide.) The Laplace operator describes
(quantum) motion of a single particle in spacetime, but it also describes
motion of the "configuration point" of some system in configuration
space.
The same applies to its square root, the Dirac operator. For
a single
particle, replacing its Laplace operator by the Dirac operator is known
to add spin to the particle. For more complex systems, replacing the
Laplace operator on configuration space by the Dirac operator on
configuration space is known to add "world-line supersymmetry" to the
system. This is, as you know, my current motivation for being interested
in the Dirac operator and its checkerboard model. In
http://www-stud.uni-essen.de/~sb0264/sqm.ps.zip
I have included results of numerical simulations of various 1+(n<3)
dim
checkerboard models for cosmological models that do have
1+(n<3) dimensional configuration space (classically). (One might
object
here, that while the classical configuration space is really 1+(n=small)
dimensional for homogeneous cosmological models, this cannot be exactly
true for quantum mechanics, but it is certainly a good approximation.)
For general interest, I have also included the case of constant
spatial superpotential, which is equivalent to the free massive electron.
Here it is interesting to observe how the zig-zag motion reappears
both
in the probability amplitude as well as in the conserved probability
current of a wave packet of finite size. This shows constructively
that
the zig-zag idea is not unique to the checkerboard model path sum.
> > > There *are* interpretational difficulties with the common conception
of
> > > antimatter, which we might eventually get into here if there
is any
> > > interest.
> > There is interest on my part. What do I have to do to make us get
> > into this?
>
> Later! I'm finishing up a paper...
Ah, ok, so what I have to do is... wait!
--
Urs.Schreiber@uni-essen.de
---------- end of forwarded post -----------------
Viewing
message <3CB48498.170F6607@uni-essen.de>
From: Urs Schreiber (Urs.Schreiber@uni-essen.de)
Subject: Re: Physically understanding the Dirac equation and
4D.
Newsgroups: sci.physics.research
Date: 2002-04-10 19:03:09 PST
Urs Schreiber wrote:
>
> "Urs Schreiber" <Urs.Schreiber@uni-essen.de> schrieb im Newsbeitrag
> news:a8c9d3$up0$1@rs04.hrz.uni-essen.de...
>
> [...]
>
> > Hence I have taken a look at the material you present on your site,
> > that titled "The HyperDiamond Random Walk", found at
> >
> > http://www.pcisys.net/~bestwork.1/QRW/the_flow_quaternions.htm
,
> >
> > which is mostly new to me. This is very illuminating. But I'd kindly
ask
> > you to help me with a couple of problems I have with some of your
claims.
>
> George Raetz <bestwork.1@pcisys.net> and I have been discussing
these
> points in private e-mail. It looks like there was actually one minor
misprint
> on that webpage, which has now been corrected, but that otherwise
the
> construction is all right. If that turns out to be true it'll
be rather
> exciting, I think.
Today I have taken the time to think about these matters
again. Looking at the new version of the above web site I have
to admit that I do not understand how the new modification
(mostly equation (26)) is supposed to work, either. Again,
that may well be my fault. But maybe it indicates that this
point might require a more detailed derivation. Even though
everythink in sight is quite elementary, the devil, as always,
is in the details.
But I am posting this in order to make a suggestion for a more
radical modification which should solve the above mentioned
problem as well as another one which has not been addressed
yet, but which is serious:
As stated, equation (15) is not covariant. That is because of
that quaternion E sitting on the left of the spinor Q in the
rhs of equation (15):
DQ = (iE)Q .
The Dirac operator D is covariant, but the unit quaternion E
on the rhs refers to a specific frame. Under a Lorentz
transformation L one finds
L DQ = iE LQ = L E' Q
<=>
DQ = E'Q
now with
E' = L~ E L .
instead of E.
This problem disappears when the unit quaternion E is brought
to the *right* of the spinor Q. What we would want is an
equation of the form
DQ = Q(iE) .
In fact, demanding that the spinor Q be an element of the
minimal left ideal generated by the primitive projector
P = (1+y0)(1+E)/4 ,
so that
Q = Q' P ,
one sees that
DQ = Q(iE)
almost looks like the the *Dirac-Lanczos equation*. (See
hep-ph/0112317, equation (5) or [1] equation (9.36)). To be
equivalent to the Dirac-Lanczos equation, and hence to be
correct, we need to require that
D = y0 @0 + y1 @1 + y2 @2 + y3 @3
instead of
... = @0 + e1 @1 + e2 @2 + e3 @3 .
All this amounts to sorting out in which particluar
representation we are actually working here.
In an attempt to address these issues, I now redo the steps
presented on
http://www.pcisys.net/~bestwork.1/QRW/the_flow_quaternions.htm
with some suitable modifications to arrive at the correct
Dirac-Lanczos equation (this is supposed to be a suggestion
subjected to discussion):
So consider a lattice in Minkoswki space generated by a unit
cell spanned by the four (Clifford) vectors
r = (y0 + y1 + y2 + y3)/2
g = (y0 + y1 - y2 - y3)/2
b = (y0 - y1 + y2 - y3)/2
y = (y0 - y1 - y2 + y3)/2 .
(yi are the generators of the Dirac algebra {yi,yj} =
diag(+1,-1,-1,-1)_ij.)
This is Tony Smith's "hyper diamond".
(Note that I use Clifford vectors instead of quaternions.)
Now consider a "Clifford algebra-weighted" random walk along
the edges of this lattice, which is described by four Clifford
valued "amplitudes":
Kr, Kg, Kb, Ky
and such that
@r Kr = k (Kg y2 y3 + Kb y3 y1 + Ky y1 y2)
@b Kb = k (Ky y2 y3 + Kr y3 y1 + Kg y1 y2)
@g Kg = k (Kr y2 y3 + Ky y3 y1 + Kb y1 y2)
@y Ky = k (Kb y2 y3 + Kg y3 y1 + Kr y1 y2) .
(This is geometrically motivated. The generators on the rhs
are those that rotate the unit vectors corresponding to the
amplitudes into each other. "k" is some constant.)
Note that I multiply the amplitudes from the *right* by the
generators of rotation, instead of multiplying them from the
left.
Next, assume that this coupled system of differential
equations is solved by a spinor Q
Q = Q' (1+y0)(1+iE)/4
E = (y2 y3 + y3 y1 + y1 y2)/sqrt(3)
with
Kr = r Q
Kg = g Q
Kb = b Q
Ky = y Q .
This ansatz for solving the above system by means of a single
spinor Q is, as I understand it, the central idea. But note
that I have here modified it on the technical side: Q is
explicitly an algebraic Clifford spinor in a definite minial
left ideal, E squares to -1, not to +1, and the Ki are
obtained from Q by premultiplying with the Clifford basis
vectors defined above.
Substituting this ansatz into the above coupled system of
differential equations one can form one covariant expression
by summing up all four equations:
(r @r + g @g + b @b + y @y) Q = k sqrt(3) Q E
The left hand side is immediate. To see that the right hand
side comes out as indicated simply note that
r + g + b + y = y0
and that
Q y0 = Q
by construction.
The above equation is the Dirac-Lanczos-Hestenes-Guersey
equation, the algebraic version of the equation describing the
free relativistic electron. The left hand side is the flat
Dirac operator
r @r + g @g + b @b + y @y = ym @m
and the right hand side, with
k = mc / (hbar sqrt(3)) ,
is equal to the mass term
i mc / hbar Q.
As usual, there are a multitude of ways to rewrite this. If
one wants to emphasize biquaternions then premultiplying
everything with y0 and splitting off the projector P on the
right of Q to express everything in terms of the, then also
biquaternionic, Q' (compare the definitions given above) gives
Lanczos' version (also used by Baylis and others).
I think this presentation improves a little on that given on
George Raetz's web site: The factor E on the right hand side
of the equation is no longer a nuisance but a necessity.
Everything is manifestly covariant (if one recalls that
algebraic spinors are manifestly covariant when nothing
non-covariand stands on their *left* side). The role of the
quaternionic structure is clarified, the construction itself
does not depend on it. Also, it is obvious how to generalize
to arbitrary dimensions. In fact, one may easily check that
for 1+1 dimensions the above scheme reproduces the Feynman
model.
While I enjoy this, there is still some scepticism in order as
long as a central questions remains to be clarified:
How much of the Ansatz
K(r,g,b,y) = (r,g,b,y) Q
is whishful thinking? For sure, every Q that solves the system
of coupled differential equations that describe the amplitude
of the random walk on the hyper diamond lattice also solves
the Dirac equation. But what about the other way round? Does
every Q that solves the Dirac equation also describe such a
random walk. Currently I do not see why it should.
--- references:
[1] W. Baylis, Clifford (Geometric) Algebras, Birkhaeuser
(1996)
--
Urs.Schreiber@uni-essen.de
Viewing
message <206f2305.0204081122.40a82726@posting.google.com>
Von:Urs Schreiber (Urs.Schreiber@uni-essen.de)
Betrifft:Re: What is Time?
Newsgroups:sci.physics.research
Datum:2002-04-09
09:46:15 PST
chaverondier wrote:
>
> "Urs Schreiber" <Urs.Schreiber@uni-essen.de> wrote:
> > In this model (which can be derived from the 1+1 dimensional Dirac
> > equation and is thus guaranteed to be relevant for "real" physics)
a
> > massive particle travels along locally lightlike paths on the line,
> > stochastically reversing direction with a certain probability.
Its spacetime
> > diagram looks like a long move on a checkerboard and its amplitude
> > can be found by pure combinatorics.
>
> In this model do massive particles move in some sort of brownian
h
> bar/2 action steps ?
The answer to your question is contained in these two papers:
T. Jacobson and L. Schulman, Quantum stochastics: The passage from a
relativistic to a non-relativistic path integral, J. Phys. A17 (1984),
p. 375-383
B. Gaveau, T. Jacobson, M. Kac, and L. Schulman, Relativistic
Extension of the Analogy between Quantum Mechanics and Brownian
Motion, Phys. Rev. Lett. 53,5 (1984), p. 419-422
As discussed in these papers, the Feynman checkerboard model is
governed by a Poisson process where the path has a constant
probability P per unit time to reverse its direction. Inserting
physical constants one easily finds that
P = mc^2 / hbar
so that the mean free distance that the path travels at the speed of
light in one direction before reversing is the compton wavelength
c (delta t) = c / P = hbar / mc ,
as might have been expected. As long as one observes at scales much
larger than the Compton wavelength the path looks like that of a
Brownian particle without correlation of succesive steps, and one
finds the usual relationship
(delta x)^2 / (delta t) = hbar / m
that governs non-relativistic quantum mechanics. But when looking at
a
resolution at about the Compton wavelength and better, it becomes
noticable that the path is not Brownian at all scales but that there
is a limiting speed c, which leads to correlations of succesive steps
over times of the order hbar / mc^2 . This is the regime of
relativistic quantum mechanics where
(delta x) / (delta t) = c .
Hence your question whether the particle path proceeds in intervals
of
quanta of action has to be addressed with respect to a "step size"
set
by the Compton wavelength. For a relativistic particle at rest the
time dependence looks like
exp(i (mc^2/hbar) t) ,
i.e the action over a time interval t is
S = mc^2 t .
Accordingly, on the mean free path length of the relativistic particle
as described by the Feynman chekcerboard model, the action increases
by
delta S = mc^2 hbar / mc^2 = hbar.
(Hm, maybe a closer look at these matters will reveal a further factor
of 1/2 or so.)
--
Urs.Schreiber@uni-essen.de