Wednesday, October 16, 2024

The Spin Chronicles (Part 3): Spin frames

What’s a Spin Axis, Anyway?

This post was inspired by a question from Bjab, who asked: "What do you mean by ‘spin axis’?" Well, buckle up because I’m about to explain that! I’m not entirely thrilled with the answer I have right now, but hey—it works (kind of). Hopefully, in the near future, I'll come up with something that satisfies both Bjab and my restless curiosity. For now we have to be happy with the standard concept of the spin structure - representing the universe as a set of interlocked rotations - described below.

Spin structure

Where Do Spinors Live?

Imagine a strange and mysterious land—let’s call it V—a two-dimensional complex vector space. It’s not just any space, though. It has "structure." Picture it equipped with a hermitian non-degenerate scalar product that we’ll denote by (u,v)(u,v), and a bilinear antisymmetric non-degenerate form we’ll call ϵ\epsilon. Sounds fancy, right? These two should be related in a specific way, and here’s how:

There exists a special kind of basis for this space, which we’ll call “special orthonormal.” It satisfies:

  1. (e1,e1)=(e2,e2)=1(e_1,e_1) = (e_2,e_2) = 1,
  2. (e1,e2)=0(e_1, e_2) = 0,
  3. ϵ(e1,e2)=1.

“Special” indeed. The first natural question is: Can this magical basis even exist?


The Magic of C²: Relax, It’s All Possible

Yes, it exists! Here’s why. Let’s choose any basis in V. Using this, we can identify V with C2C^2 (our familiar two-dimensional complex space). In C2C^2, the scalar product (u,v)(u,v) is defined as:

(u,v)=u1∗v1+u2∗v2

where the star (*) represents complex conjugation. As for ϵ\epsilon, we define it as:

ϵ(u,v)=u1v2u2v1

Now, using a standard basis e1=(10)e_1 = \begin{pmatrix} 1 \\ 0 \end{pmatrix}, e2=(01)e_2 = \begin{pmatrix} 0 \\ 1 \end{pmatrix}, it’s easy to check that our special orthonormal basis conditions hold. So, crisis averted! We don’t need to worry about whether this basis exists—it does. Phew.

Oh, and here’s a fun fact: if (e1,e2)(e_1, e_2) is a special orthonormal basis, then (−e1,e2)(-e_1, -e_2) is also a special orthonormal basis. Isn’t symmetry beautiful?


Back to the Mysterious World of V

So, we’ve got our scalar product, our fancy ϵ\epsilon, and we’re cruising in V—home to the enigmatic spinors. But are there other “special orthonormal” bases? And if so, how many?

Let’s say e1,e2e'_1, e'_2 is another basis. There exists some matrix AA such that:

ei=ejAij(i,j=1,2)

(Yes, that’s a summation over “j”.)

Now, here’s the neat part: e1,e2e'_1, e'_2 will also be special orthonormal if and only if the matrix AA is unitary with determinant 1. Matrices like this form a group called SU(2).

By the way, SU(2) is isomorphic to the group of unit quaternions. I won’t dive into that now, but it’s a pretty cool connection between algebra and geometry.


The Spinor Connection

Let’s denote by O(V) the space of all these special orthonormal bases. The group SU(2) acts on O(V) from the right, both transitively and freely. The elements of O(V) are what we call spinor bases, and the elements of V are simply spinors. If uu is a spinor, and you change the basis using an SU(2) matrix AA, the components of the spinor change according to this formula:

ui=Ajiuj

Got it? Good! Now, let’s connect this with something real, like the lab—where actual experiments happen (sometimes with coffee stains on the data sheets).


What About Our Lab?

To bridge this to the physical world, we need some more algebra (don’t worry, just a little more). Quaternions are often used for handling 3D rotations, but since we’re sticking with matrices, here’s what we do in quantum mechanics:

For each point x=(x1,x2,x3)x = (x^1, x^2, x^3) in R3R^3, associate the following 2x2 Hermitian matrix, σ(x)\sigma(x):

σ(x)=(x3x1ix2x1+ix2x3)

Now, for any SU(2) matrix AA, it turns out:

Aσ(x)A=σ(R(A)x)

where R(A)R(A) is a 3x3 real rotation matrix, an element of the special orthogonal group SO(3). In fact, the map from AR(A)is a homomorphism. And every element of SO(3) can be generated this way. The catch? If AA and BB are two SU(2) elements, then R(A)=R(B)R(A) = R(B) if and only if B=±AB = \pm A.

Hence, SU(2) is a double cover of SO(3). It’s like the universe giving you two chances to rotate.

Note: With quaternions it would look more "natural". To each x we would associate a pure imaginary quaternion 

q(x) = xi+yj+zk, 

and for each unit quaternion u, we would have 

uq(x)u* = q(R(u)x)

where R(u) is in SO(3). 


The Spin Structure: So What About That Spin Axis?

Let’s bring this back to the real world—our lab, where we have Cartesian coordinates. Any two right-handed orthogonal systems (say, different setups in the same room) are related by a unique rotation matrix from SO(3). The set of all such frames forms a space O(L), which SO(3) acts on.

Now, here’s the kicker: nature somehow maps each element of O(V) (the spinor world) to an element of O(L) (our lab’s coordinate system) in a way that respects the rotations. We call this map the spin structure.

In other words, when you rotate an element of O(V) by an SU(2) matrix A, the corresponding lab element rotates, too, By and SO(3) matrix R(A). It’s like a cosmic interlock between spinors and real-world rotations.

Note (for experts): This is a baby version of a spin structure, good for a 3D flat space or, more generally, for a space with a distant parallelism. But this baby version is all that we need for now.


What’s Next?

All this brings us closer to answering the question: what do I mean by a spin axis? Stay tuned, because we’ll tackle that mystery in the next post! We will also take under the loop the spin structure - the most important point here.

P.S. I am experimenting with the mathematical notation on this blog. Sometimes I mess up, then correct what can be corrected.

P.S. I have uploaded the pdf version of this post to academia.edu.

18 comments:

  1. Re experiments.

    In the paragraph from the words "In fact, the map" everything started to fall apart. What a person sees is different from what can be copied to the clipboard. This is very strange. The engine cannot handle processing, perhaps due to syntax errors?

    It ate whole ordinary words, causing complete lack of meaning.

    ReplyDelete
    Replies
    1. Thank you. That probably means that is better to show just pdf?

      Delete
    2. Working directly with the HTML code is a nightmare for me!

      Delete
    3. I understand what I did wrong this time. Tomorrow's post should be better.

      Delete
    4. "I understand what ... ".
      That's great.

      Delete
    5. Thanks for the PDF, there're just small differences compared to what can be seen in the post, but it makes the text much more readable.

      In A Linear Algebraic Introduction to Quaternions file you shared the link to several posts ago, it said that Hamilton claimed and fervently defended that the quaternion product uq itself was the rotation operation and not the "conjugation" uqu*, like written in the post and what seemed to be closer to correct interpretation after some evaluation by hand.

      During the "disentanglement" some interesting things popped up that might relate to the spin structure and spin axis measurement topics at hand, and all comments on the texts below and eventual corrections are highly appreciated and welcomed.

      A bit of a background for the discussion, so we're all on the same page, sort of.

      A quaternion q=(q0+iq1), where q0 and q1 are its real/scalar and imaginary/vector parts, respectively, has a polar form and unique decomposition as,
      q = ||q|| * (cosϕ_q + i*n_q*sinϕ_q) = ||q|| * exp(i*n_q*ϕ_q) = ||q|| * U_q,
      where ||q||^2 = q0^2 + q1^2 is its norm, angle ϕ_q defined as cosϕ_q = q0/||q|| and sinϕ_q = |q1|/||q||, n_q = q1/|q1| is its unit vector or axis, and U_q = q/||q|| = exp(i*n_q*ϕ_q) is its unique unit quaternion, ||U_q||=1.
      In that respect, a pure imaginary vector quaternion q'= (0+i*q1),
      q' = i*q1 = i*|q1|*n_q1 = |q1| * (cosϕ_q1 + i*n_q1*sinϕ_q1) = |q1| * exp(i*n_q1*ϕ_q1) = |q1|*U_q1,
      shares the same unit vector axis with its general/origin quaternion q, n_q1 = n_q, which means that q can be written as q = q0 + i*n_q*|q1| = q0 + |q1|*U_q1, where U_q1 = i*n_q is imaginary unit vector of q1 or its unique unit vector/imaginary quaternion. In other words, all quaternions with their axis equal to n_q have the same U_q1.
      Using polar decomposition of q1, it's easily seen that ϕ_q1=pi/2, as cosϕ_q1=0 and sinϕ_q1=1 regardless of the value |q1| or direction of axis n_q, meaning that every pure imaginary vector quaternion q' has the form,
      q' = i*q1 = |q1|*U_q1 = |q1| * exp(i*n_q*pi/2).

      Defining the vector product, that is the product of two pure imaginary vector quaternions q' and p' as,
      q' * p' = i*q1 * i*p1 = i^2 * q1 * p1 = i^2 * (q1.p1 - i*(q1×p1)),
      where "." and "×" are Gibbs' defined and usually known dot and cross vector products in vector analysis, respectively, for a product of two general quaternions q=(q0+q') and p=(p0+p') we obtain,
      q * p = q0*p0 - q1.p1 + i*q0*p1 + i*p0*q1 + i*(q1×p1)
      = q0*p0 - |q1|*|p1|*cos(ϕ_(q1,p1)) + q0*|p1|*U_p1 + p0*|q1|*U_q1 + i*|q1|*|p1|*sin(ϕ_(q1,p1))*n⊥(q1,p1),
      where ϕ_(q1,p1) is the angle between q1 and p1 lying in the plane spanned by those vectors, and n⊥(q1,p1) is the unit vector pointing in the direction, using right-hand rule, perpendicular to that plane.

      Now that we hopefully all have the basic ropes in our hands and stuff under our belts, lets see in following comment what can be deduced from multiplication of sort of collinear quaternions, that is those quaternions sharing the same imaginary unit 3d vector or axis, n_q, like q, U_q, q' and of course U_q1.

      Delete
    6. Lets start with multiplication of pure imaginary quaternions, like U_q1 and q' for example. Since we're looking at quaternions with the same axis n_q, we don't need to worry about the fact that quaternion product is not commutative because that comes into play only for those quaternions with non-collinear imaginary 3d vector axes.
      Applying U_q1 to q', that is to q1 itself, we get,
      U_q1 * q' = i*n_q * i*n_q*|q1| = |q1| * U_q1 * U_q1 =
      = i^2 * |q1| (n_q*n_q) = -|q1|*(n_q.n_q - i*(n_q×n_q)) = -|q1|,
      which can be interpreted, taking into consideration that ϕ_q1=pi/2 and thus q'=U_q1*|q1| = |q1|*exp(i*n_q*pi/2), in a sense that U_q1 or imaginary 3d unit vector "pulls" from the 4th dimension of reals/scalars, which is perpendicular to imaginary 3d space R^3, the axis n_q perpendicular to that dimension, as for the pi/2 angle, and when applied again on the same axis, it then "pushes" it towards for pi/2 in the same direction of original action. Since all imaginary vector quaternions in general and their axes in particular are perpendicular to the real 4th dimension of scalars, it's natural that all pure imaginary vector quaternions have its ϕ_q'=pi/2 regardless of their specific axes n_q. That also suggests, as there are infinitely many possibilities to draw an axis through the point, that there might be infinitely many real directions or scalar "dimensions" of 4th dimensional character perpendicular to our imaginary 3d space.

      When U_q1 is applied again, acting or rotating the axis n_q in the same direction as the one in which it got pulled out from the 4th real dimension, the axis then points in the opposite direction in the 4th scalar dimension, as for the minus sign at the end, that is the scalar or real got rotated once for pi/2 to become imaginary vector i*q1=i*|q1|n_q, and then again for 90° to make a complete angle of pi and so become sort of its original reflection or opposite, -|q1|. In that sense, multiplication of imaginary vectors with themselves, yields q'*q' = |q1|^2 * U_q1*U_q1 = -|q1|^2.

      The situation clears itself a bit more when we take a look at what U_q1 does when acting on its origin general quaternion q,
      U_q1 * q = i*n_q * (q0 + q') = q0 * i*n_q - |q1| = -|q1| + i*n_q*q0,
      indicating that while it "pushed" imaginary vector q1 into the 4th dimensional scalar in the opposite direction, at the same time it "pulled" the scalar quaternion part q0 from the real 4th dimension into imaginary 3d space, pointing in its characteristic imaginary 3d direction n_q. In mathematical jargon, we could say that it sort of exchanged the values of real and imaginary parts, if not for that minus sign, but when applied again, as indicated in previous passage, it turns its origin general quaternions completely to their opposites,
      U_q1 * U_q1 * q = -1 * q.

      Delete
    7. Other three possible combinations of multiplying general with their vector quaternions, that is q*q1, U_q*U_q1 and U_q*q1, can be easily evaluated/deduced using polar decomposition, and all three, that is all four of them in total, have identical resulting unit quaternion of the form U_q*U_q1, so lets see its exact form on the example above,
      U_q1 * q = -|q1| + i*n_q*q0 = ||q|| * exp(i*n_q*pi/2) * exp(i*n_q*ϕ_q) =
      = ||q|| * (-|q1|/||q|| + i*n_q * q0/||q||) = ||q|| * exp(i*n_q*(pi/2 + ϕ_q)) =
      = ||q||*(-sin(ϕ_q) + i*n_q*cos(ϕ_q)) = ||q||*(cos(ϕ_q + pi/2) + i*n_q*sin(ϕ_q + pi/2),
      where, for the first term in the third line, relations for cos(ϕ_q) and sin(ϕ_q) from previous comment were used. Like discussed in the passage above, action of a unit vector on its original general quaternion just "adds" an additional pi/2 to its original characteristic angle ϕ_q which is evaluated with respect to the real 4th scalar dimension, going counter-clockwise as indicated by the relation for its real part, Re(q) = q0 = ||q||*cos(ϕ_q).

      In that sense, Hamilton was not actually completely wrong when he argued that multiplication of quaternions acts as a rotation; it sort of does, but not only strictly in imaginary 3d space R^3, but in relation to complete 4 dimensional quaternion space, including also the real 4th dimension of scalars.
      To complete this part of multiplication of "collinear" quaternions, and to drive the point/discussion home that this multiplication sort of speaks about rotation with respect to the 4th scalar dimension, lets check an example of a general quaternion product. Using the expression from the previous comment, for q*q we get,
      q * q = q0^2 - |q1|^2 + 2*q0*|q1| * n_q*i = ||q||^2 * U_q * U_q =
      = ||q||^2 * (q0^2/||q||^2 - |q1|^2/||q||^2 + 2*(q0*|q1|/||q||^2)) * n_q*i =
      = ||q||^2 * (cos^2(ϕ_q) - sin^2(ϕ_q) + 2*cos(ϕ_q)*sin(ϕ_q) * n_q*i) =
      = ||q||^2 * (cos(2*ϕ_q) + sin(2*ϕ_q) * n_q*i) = ||q||^2 * exp(i*n_q*(2*ϕ_q)),
      showing that the angle of quaternion squared, or its unit quaternion squared U_q*U_q, is in fact 2*ϕ_q, that is application of a unit quaternion to its origin general quaternion "moves" the origin quaternion for an additional value of its characteristic angle ϕ_q, "from" the 4th scalar dimension and "into" the imaginary 3d space R^3. It can be reasonably deduced from the discussion above that this movement or rotation "from" 4th dimension, perpendicular to 3d space R^3, turns into movement "to" when the overall value of rotation angle value reaches pi/2.

      Next comes the general application of quaternion multiplication, like an arbitrary (U_q * p) for example, where axes n_q and n_p are non-collinear, but before diving into that a bit more complex part, our imaginary 3d world takes its "toll", so its time for a dinner and an energy refill, sort of. :)

      Delete
  2. В квантовой механике применяются комплексные волновые функции, но для прояснения алгебры этих функций можно обойтись и без комплексных чисел. Именно об этом ваш пример с кватернионами. Наверно, дополнительно можно было бы привести пример с 4-мерными действительными кососимметрическими матрицами. Но с точки зрения физики важно дать интерпретацию этих матриц.

    ReplyDelete
    Replies
    1. Ну раз сказал а, то уже скажу и бэ. На заметку физикам: эрмитовы матрицы, чисто мнимые кватернионы и кососимметрические матрицы изоморфны алгебре Ли векторных полей Киллинга 4-мерного евклидова пространства, которая эквивалентна алгебре Ли линейных векторных полей, касательных к трёхмерным сферам с центром в нулевой точке.

      Delete
    2. Do you get (1.1.6) from your matrices? Or you get something else?

      Delete
    3. Вы меня спросили, что я ещё (кроме алгебры кватернионов) получаю из этих матриц и я вам ответил, что ничего. Однако, это не означает, что из алгебры кватернионов нельзя получить алгебру Ли. Наоборот каждой алгебре Клиффорда соответствует своя алгебра Ли. Во введении об этом я просто не пишу.

      Delete
    4. So, how did you get (1.1.6)? Where from?

      Delete
  3. SO(3) aka SU(2) in multivector form is XZ, YZ, XY which has a transverse-longitudinal structure of sorts. For quaternions aka Cl(2) it would be X, Y, XY plus the scalar aka also a transverse-longitudinal (plus scalar/timelike) structure. Spinors are kind of vector-like in the vector-half spinor-half spinor triality sense at 8-dim which 4-dim might be a dimensional reduction from.

    I tend to think of everything in a multivector sense and via Tony Smith I have cellular automata pictures in mind as an analogy for multivectors. The transverse-longitudinal-timelike structure of spacetime is plastered all over the place for multivectors.

    ReplyDelete
  4. "For each point
    x=(x1,x2,x3) in R3, associate the following 2x2 Hermitian matrix,
    σ(x):
    σ(x)=(x3 x1−ix2, x1+ix2 −x3)"

    Could we use Anti-Hermitian matrix instead?
    Like iσ(x) ?

    ReplyDelete
  5. We could, indeed. It is sufficient to multiply by "i" or by "-i". Then we would have "quaternion algebra." But here, for applications in QM, I prefer Hermitian matrices.

    ReplyDelete
  6. Аркадиуш, позвольте полюбопытствовать,- вы сами рисуете картинки к своим постам?

    ReplyDelete
    Replies
    1. No. It is the computer that is painting them. I am only providing instructions to the computer what to paint and how.

      Delete

Thank you for your comment..

The Spin Chronicles (Part 3): Spin frames

What’s a Spin Axis, Anyway? This post was inspired by a question from Bjab, who asked: "What do you mean by ‘spin axis’?" Well, bu...