Showing posts with label vector fields. Show all posts
Showing posts with label vector fields. Show all posts

Friday, April 4, 2025

Lie Sphere Geometry Part 3: oriented circles

 This post is a continuation of Lie Sphere Geometry Part 2: unoriented circles. We will move to oriented circles now. Here I am trying to follow the exposition of this subject as it is described in Ch. 15. 1, Oriented circles in S3 of [1]. Instead of S3 I am taking S2 first. I say "I am trying to follow" instead of "I am following". The reason for that is that I am not entirely happy with my understanding of the exposition in this source. 

 I am not entirely happy with my understanding of the exposition in this source.


Of course there are other sources discussing the Lie Sphere Geometry, but each author has a different approach, different angle of approach, different main concepts, and I have chosen Jensen as best fitting to our purpose, at least for now. In this post we will discuss oriented circles on the sphere S2. Let us recall the definition from the previous post.

S2 = {xR3 x2 = 1}


For every mS2, and for every r∈[0,π] the unoriented circle Sr(m)is defined as

Sr(m) = {xS2 x·m = cos r}.

We have then

Sr(m)=Sπ-r(-m).


We now extend this definition of Sr(m) any r∈[0,2π]. Of course it makes sense, but Sr(m) = S2π-r(m). To remove this degeneracy, and following Ref. [1],  we define the oriented circle as follows:

Definition. An oriented circle in S2 is a circle Sr(m), together with a choice of  continuous unit normal vector field on it.

Ref. [1] expands this definition as follows.

We can use the radius to define an orientation of Sr(m) as follows. Recall that x∈Sr(m) if and only if x·m = cos(r). Let n∈TxS2 be the unit vector satisfying the equation

m = cos(r)x + sin(r)n.         (0)

Except for the cases r = aπ for any integer a, this uniquely determines n for each point x on the locus. The cases r = aπ, a any integer, are the point spheres m and -m, and there is no condition on n in

n(x) = (m - cos(r) x)/sin(r).                (1)

We consider point spheres to be in the set of all oriented spheres, but without orientation. Thus, the oriented spheres with
center m are parametrized by their radius r satisfying 0 ≤ r < 2π.

I am having some problem with the formula (1) as r approaches 0. The reasoning above was causing me some stomach problems. Therefore I suggest a somewhat different approach. The rest of this post is my own thinking about the subject.

So, let us consider the family of circles for a fixed m. In order to easily visualize thi family let us take m = (1,0,0). I will also use the letter t instead of r. That is because I want to treat the radius as a parameter of some kind of a dynamics. You will see it below.

The parametric equations x(t,ϕ) of the family are then:

x(t,ϕ) = cos(t),
y(t,ϕ) = sin(t)cos(ϕ),
z(t,ϕ) = sin(t)sin(ϕ).

Here ϕ is a parameter along the circle. For t=0 we have the point m=(1,0,0). We notice that

x(t,ϕ)2 =  x(t,ϕ)2+y(t,ϕ)2+z(t,ϕ)2 ≡1,

so that the circles of constant t are indeed all on S2.

Now, for a fixed ϕ calculate the partial derivative ∂x(t,ϕ)/∂t with respect to t. We will call it the "velocity vector field" and denote it n(t,ϕ):

n(t,ϕ) = ∂x(t,ϕ)/∂t,

or, explicitly

nx(t,ϕ) = - sin(t),
ny(t,ϕ) = cos(t)cos(ϕ),
nz(t,ϕ) = cos(t)sin(ϕ).

We see that n(t,ϕ)2 = 1, so we have a field of unit vectors. Moreover n(t,ϕ)is tangent to the sphere at x(t,ϕ), as it is tangent to the great circle of constant  ϕ.. Moreover, it is perpendicular to the circle of constant t. Let us calculate the numerator of (1) in our parametrization:

mx - cos(t)x(t,ϕ) = 1 -cos2(t) = sin2(t) = -sin(t)nx(t,ϕ) ,
my - cos(t)y(t,ϕ) = 0 - cos(t)sin(t)cos(ϕ) = - sin(t)ny(t,ϕ) ,
mz - cos(t)z(t,ϕ) = 0 - cos(t)sin(t)sin(ϕ). = - sin(t)nz(t,ϕ) .

Comparing with (1) we see that our n(t,ϕ) coincides with -n(x) of (1), but it makes a perfect sense also for t=0 and t=π.

Simulation of n(t,ϕ) for varying t, and ϕ=0, π/2,π,3π/2. 


To be continued ...

References

[1] G.R. Jensen et al., Surfaces in Classical Geometries, Springer 2016.


Thursday, October 10, 2024

The Quirks of Quaternions

The Spark of Curiosity

This post is inspired by a fascinating conversation I had with Igor Bayak and Bjab in the comments of my previous blog. Their input sparked my curiosity to dive deep into quaternions and vector fields on a three-dimensional sphere. 


These concepts aren’t just abstract math—they could be quite handy for my spinor studies. Plus, there’s something aesthetically satisfying about painting these mathematical fields. And let’s be honest, who doesn’t love a little visual beauty in math?

So, here’s a look at the world of quaternions through my curious eyes.


Quaternions 101: Meet i, j, and k

Let’s start with our three quirky quaternion friends: i, j, and k. These are the building blocks of quaternions. A general quaternion, x, can be written like this:

x=x1i+x2j+x3k+x4

Simple enough, right? Now, it’s time to turn these quaternions into something a little more structured—a matrix representation. (Cue dramatic music.)


Quaternion Matrix Magic

We’re going to multiply x by i, j, and k from the left and see what matrices pop out. This is where the magic happens. Let’s start with i.

When we multiply i by x, we get:

ix=x1+x2kx3j+x4i

Now, for the matrix interpretation:

  • At i, we have x^4, which gives us the first row of the matrix: {0, 0, 0, 1}.
  • At j, we have -x^3, giving the second row: {0, 0, -1, 0}.
  • At k, we have x^2, producing the third row: {0, 1, 0, 0}.
  • Finally, at unity, we get -x^1, completing the fourth row: {-1, 0, 0, 0}.

Putting it all together, the matrix that represents multiplication by i on the left is:

L1=(0001001001001000)L1 = \begin{pmatrix} 0 & 0 & 0 & 1 \\ 0 & 0 & -1 & 0 \\ 0 & 1 & 0 & 0 \\ -1 & 0 & 0 & 0 \end{pmatrix}

Matrix magic! (Applause, please.)


Now Multiply by j and k

Using the same process, we get matrices for multiplication by j and k on the left:

For j, we get:

L2=(0010000110000100)L2 = \begin{pmatrix} 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ -1 & 0 & 0 & 0 \\ 0 & -1 & 0 & 0 \end{pmatrix}

For k, we have:

L3=(0100100000010010)L3 = \begin{pmatrix} 0 & -1 & 0 & 0 \\ 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 0 & -1 & 0 \end{pmatrix}


What Happens on the Right Side?

Not to be left out (pun intended), we can also multiply quaternions from the right. When we do this, we get the following matrices:

For i on the right:

R1=(0001001001001000)R1 = \begin{pmatrix} 0 & 0 & 0 & 1 \\ 0 & 0 & 1 & 0 \\ 0 & -1 & 0 & 0 \\ -1 & 0 & 0 & 0 \end{pmatrix}

For j:

R2=(0010000110000100)R2 = \begin{pmatrix} 0 & 0 & -1 & 0 \\ 0 & 0 & 0 & 1 \\ 1 & 0 & 0 & 0 \\ 0 & -1 & 0 & 0 \end{pmatrix}

For k:

R3=(0100100000010010)R3 = \begin{pmatrix} 0 & 1 & 0 & 0 \\ -1 & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 0 & -1 & 0 \end{pmatrix}

So now we’ve got both the left and right multiplication matrices. Our quaternion friends are getting quite versatile.


Quaternions, SU(2), and a Three-Dimensional Sphere

Now things get even more interesting. Quaternions with a norm of 1 form a three-dimensional sphere. And not only that, but they form a group that is isomorphic to SU(2)—fancy math speak for “this group behaves like SU(2).”

To make this more concrete, consider the action of this group on the space of all quaternions. Let’s define an action u as:

u:xuxu: x \rightarrow ux

We now have a representation of this group acting on functions, written as:

(T(u)f)(x)=f(u1x)

Are you still with me? Good! Let’s move on.


Vector Fields: Expanding the Fun

Consider a one-parameter subgroup, say exp(ti). This subgroup generates a vector field, which we’ll call X1:

X1(f)=ddtf(exp(ti)x)t=0

When we calculate the action of X1 on the coordinate functions x^i, we get the components of X1:

X1(x)=L1jixjX^i (x) = - L1^i_j x^j

In the same way, we can derive the vector fields X2 and X3

Here are the results:


The Takeaway

Quaternions may sound intimidating at first, but once you break them down, they’re not only manageable—they’re downright fascinating. Their matrix representations, their connection to SU(2), and the rich vector fields they generate all have deep implications in both math and physics. And who knew matrices could be so much fun?

So, the next time you’re pondering the mysteries of the universe (or trying to impress someone at a party), just casually drop some knowledge about quaternion vector fields. It’s sure to be a hit!

Biolocation

  On Tuesday, December 23, Vlad Zhigalov (see e.g. here ) had a talk at the " Temporology " seminar hosted at Omsk.  He spoke abo...