rotation

Slicing of the 3D Mandelbrot set, and analysis.

February 8, 2021 math and physics play , , , , , , ,

Some slices.

As followup to:

here is a bit more experimentation with Paraview slice filtering. This time I saved all the point data with the escape time counts, and rendered it with a few different contours

The default slice filter places the plane in the x-y orientation:

but we can also tilt it in the Paraview render UI

and suppress the contour view to see just the slice

As a very GUI challenged user, I don’t find the interface particularly intuitive, but have at least figured out this one particular slicing task, which is kind of cool.  It’s impressive that the UI can drive interesting computational tasks without having to regenerate or reload any of the raw data itself.  This time I was using the MacOSX Paraview client, which is nicer looking than the Windows version, but has some weird glitches in the file dialogues.

Analysis.

The graphing play above shows some apparent rotational symmetry our vector equivalent to the Mandelbrot equation
\begin{equation}
\Bx \rightarrow \Bx \Be_1 \Bx + \Bc.
\end{equation}
It was not clear to me if this symmetry existed, as there were artifacts in the plots that made it appear that there was irregularity. However, some thought shows that this irregularity is strictly due to sampling error, and perhaps also due to limitations in the plotting software, as such an uneven surface is probably tricky to deal with.

To see this, here are the first few iterations of the Mandlebrot sequence for an arbitary starting vector \( \Bc \).
\begin{equation}
\begin{aligned}
\Bx_0 &= \Bc \\
\Bx_1 &= \Bc \Be_1 \Bc + \Bc \\
\Bx_2 &= \lr{ \Bc \Be_1 \Bc + \Bc } \Be_1 \lr{ \Bc \Be_1 \Bc + \Bc } + \Bc \Be_1 \Bc + \Bc.
\end{aligned}
\end{equation}

Now, what happens when we rotate the starting vector \( \Bc \) in the \( y-z \) plane. The rotor for such a rotation is
\begin{equation}
R = \exp\lr{ e_{23} \theta/2 },
\end{equation}
where
\begin{equation}
\Bc \rightarrow R \Bc \tilde{R}.
\end{equation}
Observe that if \( \Bc \) is parallel to the x-axis, then this rotation leaves the starting point invariant, as \( \Be_1 \) commutes with \( R \). That is
\begin{equation}
R \Be_1 \tilde{R} =
\Be_1 R \tilde{R} = \Be_1.
\end{equation}
Let \( \Bc’ = R \Bc \tilde{R} \), so that
\begin{equation}
\Bx_0′ = R \Bc \tilde{R} = R \Bx_0 \tilde{R} .
\end{equation}
\begin{equation}
\begin{aligned}
\Bx_1′
&= R \Bc \tilde{R} \Be_1 R \Bc \tilde{R} + R \Bc \tilde{R} \\
&= R \Bc \Be_1 \Bc \tilde{R} + R \Bc \tilde{R} \\
&= R \lr{ \Bc \Be_1 \Bc R + \Bc } \tilde{R} \\
&= R \Bx_1 \tilde{R}.
\end{aligned}
\end{equation}
\begin{equation}\label{eqn:m2:n}
\begin{aligned}
\Bx_2′
&= \Bx_1′ \Be_1 \Bx_1′ + \Bc’ \\
&= R \Bx_1 \tilde{R} \Be_1 R \Bx_1 \tilde{R} + R \Bc \tilde{R} \\
&= R \Bx_1 \Be_1 \Bx_1 \tilde{R} + R \Bc \tilde{R} \\
&= R \lr{ \Bx_1 \Be_1 \Bx_1 + \Bc } \tilde{R} \\
&= R \Bx_2 \tilde{R}.
\end{aligned}
\end{equation}

The pattern is clear. If we rotate the starting point in the y-z plane, iterating the Mandelbrot sequence results in precisely the same rotation of the x-y plane Mandelbrot sequence. So the apparent rotational symmetry in the 3D iteration of the Mandelbrot vector equation is exactly that. This is an unfortunately boring 3D fractal. All of the interesting fractal nature occurs in the 2D plane, and the rest is just a consequence of rotating that image around the x-axis. We get some interesting fractal artifacts if we slice the rotated Mandelbrot image.

A couple more reciprocal frame examples.

December 14, 2020 math and physics play , , , , , , , , , , , , ,

[If mathjax doesn’t display properly for you, click here for a PDF of this post]

This post logically follows both of the following:

  1. Curvilinear coordinates and gradient in spacetime, and reciprocal frames, and
  2. Lorentz transformations in Space Time Algebra (STA)

The PDF linked above above contains all the content from this post plus (1.) above [to be edited later into a more logical sequence.]

More examples.

Here are a few additional examples of reciprocal frame calculations.

Problem: Unidirectional arbitrary functional dependence.

Let
\begin{equation}\label{eqn:reciprocal:2540}
x = a f(u),
\end{equation}
where \( a \) is a constant vector and \( f(u)\) is some arbitrary differentiable function with a non-zero derivative in the region of interest.

Answer

Here we have just a single tangent space direction (a line in spacetime) with tangent vector
\begin{equation}\label{eqn:reciprocal:2400}
\Bx_u = a \PD{u}{f} = a f_u,
\end{equation}
so we see that the tangent space vectors are just rescaled values of the direction vector \( a \).
This is a simple enough parameterization that we can compute the reciprocal frame vector explicitly using the gradient. We expect that \( \Bx^u = 1/\Bx_u \), and find
\begin{equation}\label{eqn:reciprocal:2420}
\inv{a} \cdot x = f(u),
\end{equation}
but for constant \( a \), we know that \( \grad a \cdot x = a \), so taking gradients of both sides we find
\begin{equation}\label{eqn:reciprocal:2440}
\inv{a} = \grad f = \PD{u}{f} \grad u,
\end{equation}
so the reciprocal vector is
\begin{equation}\label{eqn:reciprocal:2460}
\Bx^u = \grad u = \inv{a f_u},
\end{equation}
as expected.

Problem: Linear two variable parameterization.

Let \( x = a u + b v \), where \( x \wedge a \wedge b = 0 \) represents spacetime plane (also the tangent space.) Find the curvilinear coordinates and their reciprocals.

Answer

The frame vectors are easy to compute, as they are just
\begin{equation}\label{eqn:reciprocal:1960}
\begin{aligned}
\Bx_u &= \PD{u}{x} = a \\
\Bx_v &= \PD{v}{x} = b.
\end{aligned}
\end{equation}
This is an example of a parametric equation that we can easily invert, as we have
\begin{equation}\label{eqn:reciprocal:1980}
\begin{aligned}
x \wedge a &= – v \lr{ a \wedge b } \\
x \wedge b &= u \lr{ a \wedge b },
\end{aligned}
\end{equation}
so
\begin{equation}\label{eqn:reciprocal:2000}
\begin{aligned}
u
&= \inv{ a \wedge b } \cdot \lr{ x \wedge b } \\
&= \inv{ \lr{a \wedge b}^2 } \lr{ a \wedge b } \cdot \lr{ x \wedge b } \\
&=
\frac{
\lr{b \cdot x} \lr{ a \cdot b }

\lr{a \cdot x} \lr{ b \cdot b }
}{ \lr{a \wedge b}^2 }
\end{aligned}
\end{equation}
\begin{equation}\label{eqn:reciprocal:2020}
\begin{aligned}
v &= -\inv{ a \wedge b } \cdot \lr{ x \wedge a } \\
&= -\inv{ \lr{a \wedge b}^2 } \lr{ a \wedge b } \cdot \lr{ x \wedge a } \\
&=
-\frac{
\lr{b \cdot x} \lr{ a \cdot a }

\lr{a \cdot x} \lr{ a \cdot b }
}{ \lr{a \wedge b}^2 }
\end{aligned}
\end{equation}
Recall that \( \grad \lr{ a \cdot x} = a \), if \( a \) is a constant, so our gradients are just
\begin{equation}\label{eqn:reciprocal:2040}
\begin{aligned}
\grad u
&=
\frac{
b \lr{ a \cdot b }

a
\lr{ b \cdot b }
}{ \lr{a \wedge b}^2 } \\
&=
b \cdot \inv{ a \wedge b },
\end{aligned}
\end{equation}
and
\begin{equation}\label{eqn:reciprocal:2060}
\begin{aligned}
\grad v
&=
-\frac{
b \lr{ a \cdot a }

a \lr{ a \cdot b }
}{ \lr{a \wedge b}^2 } \\
&=
-a \cdot \inv{ a \wedge b }.
\end{aligned}
\end{equation}
Expressed in terms of the frame vectors, this is just
\begin{equation}\label{eqn:reciprocal:2080}
\begin{aligned}
\Bx^u &= \Bx_v \cdot \inv{ \Bx_u \wedge \Bx_v } \\
\Bx^v &= -\Bx_u \cdot \inv{ \Bx_u \wedge \Bx_v },
\end{aligned}
\end{equation}
so we were able to show, for this special two parameter linear case, that the explicit evaluation of the gradients has the exact structure that we intuited that the reciprocals must have, provided they are constrained to the spacetime plane \( a \wedge b \). It is interesting to observe how this structure falls out of the linear system solution so directly. Also note that these reciprocals are not defined at the origin of the \( (u,v) \) parameter space.

Problem: Quadratic two variable parameterization.

Now consider a variation of the previous problem, with \( x = a u^2 + b v^2 \). Find the curvilinear coordinates and their reciprocals.

Answer

\begin{equation}\label{eqn:reciprocal:2100}
\begin{aligned}
\Bx_u &= \PD{u}{x} = 2 u a \\
\Bx_v &= \PD{v}{x} = 2 v b.
\end{aligned}
\end{equation}
Our tangent space is still the \( a \wedge b \) plane (as is the surface itself), but the spacing of the cells starts getting wider in proportion to \( u, v \).
Utilizing the work from the previous problem, we have
\begin{equation}\label{eqn:reciprocal:2120}
\begin{aligned}
2 u \grad u &=
b \cdot \inv{ a \wedge b } \\
2 v \grad v &=
-a \cdot \inv{ a \wedge b }.
\end{aligned}
\end{equation}
A bit of rearrangement can show that this is equivalent to the reciprocal frame identities. This is a second demonstration that the gradient and the algebraic formulations for the reciprocals match, at least for these special cases of linear non-coupled parameterizations.

Problem: Reciprocal frame for generalized cylindrical parameterization.

Let the vector parameterization be \( x(\rho,\theta) = \rho e^{-i\theta/2} x(\rho_0, \theta_0) e^{i \theta} \), where \( i^2 = \pm 1 \) is a unit bivector (\(+1\) for a boost, and \(-1\) for a rotation), and where \(\theta, \rho\) are scalars. Find the tangent space vectors and their reciprocals.

fig. 1. “Cylindrical” boost parameterization.

Note that this is cylindrical parameterization for the rotation case, and traces out hyperbolic regions for the boost case. The boost case is illustrated in fig. 1 where hyperbolas in the light cone are found for boosts of \( \gamma_0\) with various values of \(\rho\), and the spacelike hyperbolas are boosts of \( \gamma_1 \), again for various values of \( \rho \).

Answer

The tangent space vectors are
\begin{equation}\label{eqn:reciprocal:2480}
\Bx_\rho = \frac{x}{\rho},
\end{equation}
and

\begin{equation}\label{eqn:reciprocal:2500}
\begin{aligned}
\Bx_\theta
&= -\frac{i}{2} x + x \frac{i}{2} \\
&= x \cdot i.
\end{aligned}
\end{equation}
Recall that \( x \cdot i \) lies perpendicular to \( x \) (in the plane \( i \)), as illustrated in fig. 2. This means that \( \Bx_\rho \) and \( \Bx_\theta \) are orthogonal, so we can find the reciprocal vectors by just inverting them
\begin{equation}\label{eqn:reciprocal:2520}
\begin{aligned}
\Bx^\rho &= \frac{\rho}{x} \\
\Bx^\theta &= \frac{1}{x \cdot i}.
\end{aligned}
\end{equation}

fig. 2. Projection and rejection geometry.

Parameterization of a general linear transformation.

Given \( N \) parameters \( u^0, u^1, \cdots u^{N-1} \), a general linear transformation from the parameter space to the vector space has the form
\begin{equation}\label{eqn:reciprocal:2160}
x =
{a^\alpha}_\beta \gamma_\alpha u^\beta,
\end{equation}
where \( \beta \in [0, \cdots, N-1] \) and \( \alpha \in [0,3] \).
For such a general transformation, observe that the curvilinear basis vectors are
\begin{equation}\label{eqn:reciprocal:2180}
\begin{aligned}
\Bx_\mu
&= \PD{u^\mu}{x} \\
&= \PD{u^\mu}{}
{a^\alpha}_\beta \gamma_\alpha u^\beta \\
&=
{a^\alpha}_\mu \gamma_\alpha.
\end{aligned}
\end{equation}
We find an interpretation of \( {a^\alpha}_\mu \) by dotting \( \Bx_\mu \) with the reciprocal frame vectors of the standard basis
\begin{equation}\label{eqn:reciprocal:2200}
\begin{aligned}
\Bx_\mu \cdot \gamma^\nu
&=
{a^\alpha}_\mu \lr{ \gamma_\alpha \cdot \gamma^\nu } \\
&=
{a^\nu}_\mu,
\end{aligned}
\end{equation}
so
\begin{equation}\label{eqn:reciprocal:2220}
x = \Bx_\mu u^\mu.
\end{equation}
We are able to reinterpret \ref{eqn:reciprocal:2160} as a contraction of the tangent space vectors with the parameters, scaling and summing these direction vectors to characterize all the points in the tangent plane.

Theorem 1.1: Projecting onto the tangent space.

Let \( T \) represent the tangent space. The projection of a vector onto the tangent space has the form
\begin{equation}\label{eqn:reciprocal:2560}
\textrm{Proj}_{\textrm{T}} y = \lr{ y \cdot \Bx^\mu } \Bx_\mu = \lr{ y \cdot \Bx_\mu } \Bx^\mu.
\end{equation}

Start proof:

Let’s designate \( a \) as the portion of the vector \( y \) that lies outside of the tangent space
\begin{equation}\label{eqn:reciprocal:2260}
y = y^\mu \Bx_\mu + a.
\end{equation}
If we knew the coordinates \( y^\mu \), we would have a recipe for the projection.
Algebraically, requiring that \( a \) lies outside of the tangent space, is equivalent to stating \( a \cdot \Bx_\mu = a \cdot \Bx^\mu = 0 \). We use that fact, and then take dot products
\begin{equation}\label{eqn:reciprocal:2280}
\begin{aligned}
y \cdot \Bx^\nu
&= \lr{ y^\mu \Bx_\mu + a } \cdot \Bx^\nu \\
&= y^\nu,
\end{aligned}
\end{equation}
so
\begin{equation}\label{eqn:reciprocal:2300}
y = \lr{ y \cdot \Bx^\mu } \Bx_\mu + a.
\end{equation}
Similarly, the tangent space projection can be expressed as a linear combination of reciprocal basis elements
\begin{equation}\label{eqn:reciprocal:2320}
y = y_\mu \Bx^\mu + a.
\end{equation}
Dotting with \( \Bx_\mu \), we have
\begin{equation}\label{eqn:reciprocal:2340}
\begin{aligned}
y \cdot \Bx^\mu
&= \lr{ y_\alpha \Bx^\alpha + a } \cdot \Bx_\mu \\
&= y_\mu,
\end{aligned}
\end{equation}
so
\begin{equation}\label{eqn:reciprocal:2360}
y = \lr{ y \cdot \Bx^\mu } \Bx_\mu + a.
\end{equation}
We find the two stated ways of computing the projection.

Observe that, for the special case that all of \( \setlr{ \Bx_\mu } \) are orthogonal, the equivalence of these two projection methods follows directly, since
\begin{equation}\label{eqn:reciprocal:2380}
\begin{aligned}
\lr{ y \cdot \Bx^\mu } \Bx_\mu
&=
\lr{ y \cdot \inv{\Bx_\mu} } \inv{\Bx^\mu} \\
&=
\lr{ y \cdot \frac{\Bx_\mu}{\lr{\Bx_\mu}^2 } } \frac{\Bx^\mu}{\lr{\Bx^\mu}^2} \\
&=
\lr{ y \cdot \Bx_\mu } \Bx^\mu.
\end{aligned}
\end{equation}

End proof.

Lorentz transformations in Space Time Algebra (STA)

December 12, 2020 math and physics play , , , , , , , , , , , , , , , , , ,

[If mathjax doesn’t display properly for you, click here for a PDF of this post]

Motivation.

One of the remarkable features of geometric algebra are the complex exponential sandwiches that can be used to encode rotations in any dimension, or rotation like operations like Lorentz transformations in Minkowski spaces. In this post, we show some examples that unpack the geometric algebra expressions for Lorentz transformations operations of this sort. In particular, we will look at the exponential sandwich operations for spatial rotations and Lorentz boosts in the Dirac algebra, known as Space Time Algebra (STA) in geometric algebra circles, and demonstrate that these sandwiches do have the desired effects.

Lorentz transformations.

Theorem 1.1: Lorentz transformation.

The transformation
\begin{equation}\label{eqn:lorentzTransform:580}
x \rightarrow e^{B} x e^{-B} = x’,
\end{equation}
where \( B = a \wedge b \), is an STA 2-blade for any two linearly independent four-vectors \( a, b \), is a norm preserving, that is
\begin{equation}\label{eqn:lorentzTransform:600}
x^2 = {x’}^2.
\end{equation}

Start proof:

The proof is disturbingly trivial in this geometric algebra form
\begin{equation}\label{eqn:lorentzTransform:40}
\begin{aligned}
{x’}^2
&=
e^{B} x e^{-B} e^{B} x e^{-B} \\
&=
e^{B} x x e^{-B} \\
&=
x^2 e^{B} e^{-B} \\
&=
x^2.
\end{aligned}
\end{equation}

End proof.

In particular, observe that we did not need to construct the usual infinitesimal representations of rotation and boost transformation matrices or tensors in order to demonstrate that we have spacetime invariance for the transformations. The rough idea of such a transformation is that the exponential commutes with components of the four-vector that lie off the spacetime plane specified by the bivector \( B \), and anticommutes with components of the four-vector that lie in the plane. The end result is that the sandwich operation simplifies to
\begin{equation}\label{eqn:lorentzTransform:60}
x’ = x_\parallel e^{-B} + x_\perp,
\end{equation}
where \( x = x_\perp + x_\parallel \) and \( x_\perp \cdot B = 0 \), and \( x_\parallel \wedge B = 0 \). In particular, using \( x = x B B^{-1} = \lr{ x \cdot B + x \wedge B } B^{-1} \), we find that
\begin{equation}\label{eqn:lorentzTransform:80}
\begin{aligned}
x_\parallel &= \lr{ x \cdot B } B^{-1} \\
x_\perp &= \lr{ x \wedge B } B^{-1}.
\end{aligned}
\end{equation}
When \( B \) is a spacetime plane \( B = b \wedge \gamma_0 \), then this exponential has a hyperbolic nature, and we end up with a Lorentz boost. When \( B \) is a spatial bivector, we end up with a single complex exponential, encoding our plane old 3D rotation. More general \( B \)’s that encode composite boosts and rotations are also possible, but \( B \) must be invertible (it should have no lightlike factors.) The rough geometry of these projections is illustrated in fig 1, where the spacetime plane is represented by \( B \).

Projection and rejection geometry.

fig 1. Projection and rejection geometry.

 

What is not so obvious is how to pick \( B \)’s that correspond to specific rotation axes or boost directions. Let’s consider each of those cases in turn.

Theorem 1.2: Boost.

The boost along a direction vector \( \vcap \) and rapidity \( \alpha \) is given by
\begin{equation}\label{eqn:lorentzTransform:620}
x’ = e^{-\vcap \alpha/2} x e^{\vcap \alpha/2},
\end{equation}
where \( \vcap = \gamma_{k0} \cos\theta^k \) is an STA bivector representing a spatial direction with direction cosines \( \cos\theta^k \).

Start proof:

We want to demonstrate that this is equivalent to the usual boost formulation. We can start with decomposition of the four-vector \( x \) into components that lie in and off of the spacetime plane \( \vcap \).
\begin{equation}\label{eqn:lorentzTransform:100}
\begin{aligned}
x
&= \lr{ x^0 + \Bx } \gamma_0 \\
&= \lr{ x^0 + \Bx \vcap^2 } \gamma_0 \\
&= \lr{ x^0 + \lr{ \Bx \cdot \vcap} \vcap + \lr{ \Bx \wedge \vcap} \vcap } \gamma_0,
\end{aligned}
\end{equation}
where \( \Bx = x \wedge \gamma_0 \). The first two components lie in the boost plane, whereas the last is the spatial component of the vector that lies perpendicular to the boost plane. Observe that \( \vcap \) anticommutes with the dot product term and commutes with he wedge product term, so we have
\begin{equation}\label{eqn:lorentzTransform:120}
\begin{aligned}
x’
&=
\lr{ x^0 + \lr{ \Bx \cdot \vcap } \vcap } \gamma_0
e^{\vcap \alpha/2 }
e^{\vcap \alpha/2 }
+
\lr{ \Bx \wedge \vcap } \vcap \gamma_0
e^{-\vcap \alpha/2 }
e^{\vcap \alpha/2 } \\
&=
\lr{ x^0 + \lr{ \Bx \cdot \vcap } \vcap } \gamma_0
e^{\vcap \alpha }
+
\lr{ \Bx \wedge \vcap } \vcap \gamma_0.
\end{aligned}
\end{equation}
Noting that \( \vcap^2 = 1 \), we may expand the exponential in hyperbolic functions, and find that the boosted portion of the vector expands as
\begin{equation}\label{eqn:lorentzTransform:260}
\begin{aligned}
\lr{ x^0 + \lr{ \Bx \cdot \vcap} \vcap } \gamma_0 e^{\vcap \alpha}
&=
\lr{ x^0 + \lr{ \Bx \cdot \vcap} \vcap } \gamma_0 \lr{ \cosh\alpha + \vcap \sinh \alpha} \\
&=
\lr{ x^0 + \lr{ \Bx \cdot \vcap} \vcap } \lr{ \cosh\alpha – \vcap \sinh \alpha} \gamma_0 \\
&=
\lr{ x^0 \cosh\alpha – \lr{ \Bx \cdot \vcap} \sinh \alpha} \gamma_0
+
\lr{ -x^0 \sinh \alpha + \lr{ \Bx \cdot \vcap} \cosh \alpha } \vcap \gamma_0.
\end{aligned}
\end{equation}
We are left with
\begin{equation}\label{eqn:lorentzTransform:320}
\begin{aligned}
x’
&=
\lr{ x^0 \cosh\alpha – \lr{ \Bx \cdot \vcap} \sinh \alpha} \gamma_0
+
\lr{ \lr{ \Bx \cdot \vcap} \cosh \alpha -x^0 \sinh \alpha } \vcap \gamma_0
+
\lr{ \Bx \wedge \vcap} \vcap \gamma_0 \\
&=
\begin{bmatrix}
\gamma_0 & \vcap \gamma_0
\end{bmatrix}
\begin{bmatrix}
\cosh\alpha & – \sinh\alpha \\
-\sinh\alpha & \cosh\alpha
\end{bmatrix}
\begin{bmatrix}
x^0 \\
\Bx \cdot \vcap
\end{bmatrix}
+
\lr{ \Bx \wedge \vcap} \vcap \gamma_0,
\end{aligned}
\end{equation}
which has the desired Lorentz boost structure. Of course, this is usually seen with \( \vcap = \gamma_{10} \) so that the components in the coordinate column vector are \( (ct, x) \).

End proof.

Theorem 1.3: Spatial rotation.

Given two linearly independent spatial bivectors \( \Ba = a^k \gamma_{k0}, \Bb = b^k \gamma_{k0} \), a rotation of \(\theta\) radians in the plane of \( \Ba, \Bb \) from \( \Ba \) towards \( \Bb \), is given by
\begin{equation}\label{eqn:lorentzTransform:640}
x’ = e^{-i\theta} x e^{i\theta},
\end{equation}
where \( i = (\Ba \wedge \Bb)/\Abs{\Ba \wedge \Bb} \), is a unit (spatial) bivector.

Start proof:

Without loss of generality, we may pick \( i = \acap \bcap \), where \( \acap^2 = \bcap^2 = 1 \), and \( \acap \cdot \bcap = 0 \). With such an orthonormal basis for the plane, we can decompose our four vector into portions that lie in and off the plane
\begin{equation}\label{eqn:lorentzTransform:400}
\begin{aligned}
x
&= \lr{ x^0 + \Bx } \gamma_0 \\
&= \lr{ x^0 + \Bx i i^{-1} } \gamma_0 \\
&= \lr{ x^0 + \lr{ \Bx \cdot i } i^{-1} + \lr{ \Bx \wedge i } i^{-1} } \gamma_0.
\end{aligned}
\end{equation}
The projective term lies in the plane of rotation, whereas the timelike and spatial rejection term are perpendicular. That is
\begin{equation}\label{eqn:lorentzTransform:420}
\begin{aligned}
x_\parallel &= \lr{ \Bx \cdot i } i^{-1} \gamma_0 \\
x_\perp &= \lr{ x^0 + \lr{ \Bx \wedge i } i^{-1} } \gamma_0,
\end{aligned}
\end{equation}
where \( x_\parallel \wedge i = 0 \), and \( x_\perp \cdot i = 0 \). The plane pseudoscalar \( i \) anticommutes with \( x_\parallel \), and commutes with \( x_\perp \), so
\begin{equation}\label{eqn:lorentzTransform:440}
\begin{aligned}
x’
&= e^{-i\theta/2} \lr{ x_\parallel + x_\perp } e^{i\theta/2} \\
&= x_\parallel e^{i\theta} + x_\perp.
\end{aligned}
\end{equation}
However
\begin{equation}\label{eqn:lorentzTransform:460}
\begin{aligned}
\lr{ \Bx \cdot i } i^{-1}
&=
\lr{ \Bx \cdot \lr{ \acap \wedge \bcap } } \bcap \acap \\
&=
\lr{\Bx \cdot \acap} \bcap \bcap \acap
-\lr{\Bx \cdot \bcap} \acap \bcap \acap \\
&=
\lr{\Bx \cdot \acap} \acap
+\lr{\Bx \cdot \bcap} \bcap,
\end{aligned}
\end{equation}
so
\begin{equation}\label{eqn:lorentzTransform:480}
\begin{aligned}
x_\parallel e^{i\theta}
&=
\lr{
\lr{\Bx \cdot \acap} \acap
+
\lr{\Bx \cdot \bcap} \bcap
}
\gamma_0
\lr{
\cos\theta + \acap \bcap \sin\theta
} \\
&=
\acap \lr{
\lr{\Bx \cdot \acap} \cos\theta

\lr{\Bx \cdot \bcap} \sin\theta
}
\gamma_0
+
\bcap \lr{
\lr{\Bx \cdot \acap} \sin\theta
+
\lr{\Bx \cdot \bcap} \cos\theta
}
\gamma_0,
\end{aligned}
\end{equation}
so
\begin{equation}\label{eqn:lorentzTransform:500}
x’
=
\begin{bmatrix}
\acap & \bcap
\end{bmatrix}
\begin{bmatrix}
\cos\theta & – \sin\theta \\
\sin\theta & \cos\theta
\end{bmatrix}
\begin{bmatrix}
\Bx \cdot \acap \\
\Bx \cdot \bcap \\
\end{bmatrix}
\gamma_0
+
\lr{ x \wedge i} i^{-1} \gamma_0.
\end{equation}
Observe that this rejection term can be explicitly expanded to
\begin{equation}\label{eqn:lorentzTransform:520}
\lr{ \Bx \wedge i} i^{-1} \gamma_0 =
x –
\lr{ \Bx \cdot \acap } \acap \gamma_0

\lr{ \Bx \cdot \acap } \acap \gamma_0.
\end{equation}
This is the timelike component of the vector, plus the spatial component that is normal to the plane. This exponential sandwich transformation rotates only the portion of the vector that lies in the plane, and leaves the rest (timelike and normal) untouched.

End proof.

Problems.

Problem: Verify components relative to boost direction.

In the proof of thm. 1.2, the vector \( x \) was expanded in terms of the spacetime split. An alternate approach, is to expand as
\begin{equation}\label{eqn:lorentzTransform:340}
\begin{aligned}
x
&= x \vcap^2 \\
&= \lr{ x \cdot \vcap + x \wedge \vcap } \vcap \\
&= \lr{ x \cdot \vcap } \vcap + \lr{ x \wedge \vcap } \vcap.
\end{aligned}
\end{equation}
Show that
\begin{equation}\label{eqn:lorentzTransform:360}
\lr{ x \cdot \vcap } \vcap
=
\lr{ x^0 + \lr{ \Bx \cdot \vcap} \vcap } \gamma_0,
\end{equation}
and
\begin{equation}\label{eqn:lorentzTransform:380}
\lr{ x \wedge \vcap } \vcap
=
\lr{ \Bx \wedge \vcap} \vcap \gamma_0.
\end{equation}

Answer

Let \( x = x^\mu \gamma_\mu \), so that
\begin{equation}\label{eqn:lorentzTransform:160}
\begin{aligned}
x \cdot \vcap
&=
\gpgradeone{ x^\mu \gamma_\mu \cos\theta^b \gamma_{b 0} } \\
&=
x^\mu \cos\theta^b \gpgradeone{ \gamma_\mu \gamma_{b 0} }
.
\end{aligned}
\end{equation}
The \( \mu = 0 \) component of this grade selection is
\begin{equation}\label{eqn:lorentzTransform:180}
\gpgradeone{ \gamma_0 \gamma_{b 0} }
=
-\gamma_b,
\end{equation}
and for \( \mu = a \ne 0 \), we have
\begin{equation}\label{eqn:lorentzTransform:200}
\gpgradeone{ \gamma_a \gamma_{b 0} }
=
-\delta_{a b} \gamma_0,
\end{equation}
so we have
\begin{equation}\label{eqn:lorentzTransform:220}
\begin{aligned}
x \cdot \vcap
&=
x^0 \cos\theta^b (-\gamma_b)
+
x^a \cos\theta^b (-\delta_{ab} \gamma_0 ) \\
&=
-x^0 \vcap \gamma_0

x^b \cos\theta^b \gamma_0 \\
&=
– \lr{ x^0 \vcap + \Bx \cdot \vcap } \gamma_0,
\end{aligned}
\end{equation}
where \( \Bx = x \wedge \gamma_0 \) is the spatial portion of the four vector \( x \) relative to the stationary observer frame. Since \( \vcap \) anticommutes with \( \gamma_0 \), the component of \( x \) in the spacetime plane \( \vcap \) is
\begin{equation}\label{eqn:lorentzTransform:240}
\lr{ x \cdot \vcap } \vcap =
\lr{ x^0 + \lr{ \Bx \cdot \vcap} \vcap } \gamma_0,
\end{equation}
as expected.

For the rejection term, we have
\begin{equation}\label{eqn:lorentzTransform:280}
x \wedge \vcap
=
x^\mu \cos\theta^s \gpgradethree{ \gamma_\mu \gamma_{s 0} }.
\end{equation}
The \( \mu = 0 \) term clearly contributes nothing, leaving us with:
\begin{equation}\label{eqn:lorentzTransform:300}
\begin{aligned}
\lr{ x \wedge \vcap } \vcap
&=
\lr{ x \wedge \vcap } \cdot \vcap \\
&=
x^r \cos\theta^s \cos\theta^t \lr{ \lr{ \gamma_r \wedge \gamma_{s}} \gamma_0 } \cdot \lr{ \gamma_{t0} } \\
&=
x^r \cos\theta^s \cos\theta^t \gpgradeone{
\lr{ \gamma_r \wedge \gamma_{s} } \gamma_0 \gamma_{t0}
} \\
&=
-x^r \cos\theta^s \cos\theta^t \lr{ \gamma_r \wedge \gamma_{s}} \cdot \gamma_t \\
&=
-x^r \cos\theta^s \cos\theta^t \lr{ -\gamma_r \delta_{st} + \gamma_s \delta_{rt} } \\
&=
x^r \cos\theta^t \cos\theta^t \gamma_r

x^t \cos\theta^s \cos\theta^t \gamma_s \\
&=
\Bx \gamma_0
– (\Bx \cdot \vcap) \vcap \gamma_0 \\
&=
\lr{ \Bx \wedge \vcap} \vcap \gamma_0,
\end{aligned}
\end{equation}
as expected. Is there a clever way to demonstrate this without resorting to coordinates?

Problem: Rotation transformation components.

Given a unit spatial bivector \( i = \acap \bcap \), where \( \acap \cdot \bcap = 0 \) and \( i^2 = -1 \), show that
\begin{equation}\label{eqn:lorentzTransform:540}
\lr{ x \cdot i } i^{-1}
=
\lr{ \Bx \cdot i } i^{-1} \gamma_0
=
\lr{\Bx \cdot \acap } \acap \gamma_0
+
\lr{\Bx \cdot \bcap } \bcap \gamma_0,
\end{equation}
and
\begin{equation}\label{eqn:lorentzTransform:560}
\lr{ x \wedge i } i^{-1}
=
\lr{ \Bx \wedge i } i^{-1} \gamma_0
=
x –
\lr{\Bx \cdot \acap } \acap \gamma_0

\lr{\Bx \cdot \bcap } \bcap \gamma_0.
\end{equation}
Also show that \( i \) anticommutes with \( \lr{ x \cdot i } i^{-1} \) and commutes with \( \lr{ x \wedge i } i^{-1} \).

Answer

This problem is left for the reader, as I don’t feel like typing out my solution.

The first part of this problem can be done in the tedious coordinate approach used above, but hopefully there is a better way.

For the last (commutation) part of the problem, here is a hint. Let \( x \wedge i = n i \), where \( n \cdot i = 0 \). The result then follows easily.

PHY2403H Quantum Field Theory. Lecture 22: Dirac sea, charges, angular momentum, spin, U(1) symmetries, electrons and positrons. Taught by Prof. Erich Poppitz

December 17, 2018 phy2403 , , , , , , , ,

This post is a synopsis of the material from the second last lecture of QFT I. I missed that class, but worked from notes kindly provided by Emily Tyhurst, and Stefan Divic, filling in enough details that it made sense to me.

[Click here for an unabrided PDF of my full notes on this day’s lecture material.]

Topics covered include

  • The Hamiltonian action on single particle states showed that the Hamiltonian was an energy eigenoperator
    \begin{equation}\label{eqn:qftLecture22:140}
    H \ket{\Bp, r}
    =
    \omega_\Bp \ket{\Bp, r}.
    \end{equation}
  • The conserved Noether current and charge for spatial translations, the momentum operator, was found to be
    \begin{equation}\label{eqn:momentumDirac:260}
    \BP =
    \int d^3 x
    \Psi^\dagger (-i \spacegrad) \Psi,
    \end{equation}
    which could be written in creation and anhillation operator form as
    \begin{equation}\label{eqn:momentumDirac:261}
    \BP = \sum_{s = 1}^2
    \int \frac{d^3 q}{(2\pi)^3} \Bp \lr{
    a_\Bp^{s\dagger}
    a_\Bp^{s}
    +
    b_\Bp^{s\dagger}
    b_\Bp^{s}
    }.
    \end{equation}
    Single particle states were found to be the eigenvectors of this operator, with momentum eigenvalues
    \begin{equation}\label{eqn:momentumDirac:262}
    \BP a_\Bq^{s\dagger} \ket{0} = \Bq (a_\Bq^{s\dagger} \ket{0}).
    \end{equation}
  • The conserved Noether current and charge for a rotation was found. That charge is
    \begin{equation}\label{eqn:qftLecture22:920}
    \BJ = \int d^3 x \Psi^\dagger(x) \lr{ \underbrace{\Bx \cross (-i \spacegrad)}_{\text{orbital angular momentum}} + \inv{2} \underbrace{\mathbf{1} \otimes \Bsigma}_{\text{spin angular momentum}} } \Psi,
    \end{equation}
    where
    \begin{equation}\label{eqn:qftLecture22:260}
    \mathbf{1} \otimes \Bsigma =
    \begin{bmatrix}
    \Bsigma & 0 \\
    0 & \Bsigma
    \end{bmatrix},
    \end{equation}
    which has distinct orbital and spin angular momentum components. Unlike NRQM, we see both types of angular momentum as components of a single operator. It is argued in [3] that for a particle at rest the single particle state is an eigenvector of this operator, with eigenvalues \( \pm 1/2 \) — the Fermion spin eigenvalues!
  • We examined two \( U(1) \) global symmetries. The Noether charge for the “vector” \( U(1) \) symmetry is
    \begin{equation}\label{eqn:qftLecture22:380}
    Q
    =
    \int \frac{d^3 q}{(2\pi)^3} \sum_{s = 1}^2
    \lr{
    a_\Bp^{s \dagger} a_\Bp^s

    b_\Bp^{s \dagger}
    b_\Bp^s
    },
    \end{equation}
    This charge operator characterizes the \( a, b \) operators. \( a \) particles have charge \( +1 \), and \( b \) particles have charge \( -1 \), or vice-versa depending on convention. We call \( a \) the operator for the electron, and \( b \) the operator for the positron.
  • CPT (Charge-Parity-TimeReversal) symmetries were also mentioned, but not covered in class. We were pointed to [2], [3], [4] to start studying that topic.

References

[1] C. Doran and A.N. Lasenby. Geometric algebra for physicists. Cambridge University Press New York, Cambridge, UK, 1st edition, 2003.

[2] Dr. Michael Luke. Quantum Field Theory., 2011. URL https://www.physics.utoronto.ca/~luke/PHY2403F/References_files/lecturenotes.pdf. [Online; accessed 05-Dec-2018].

[3] Michael E Peskin and Daniel V Schroeder. An introduction to Quantum Field Theory. Westview, 1995.

[4] Dr. David Tong. Quantum Field Theory. URL http://www.damtp.cam.ac.uk/user/tong/qft.html.

PHY2403H Quantum Field Theory. Lecture 9: Unbroken and spontaneously broken symmetries, Higgs Lagrangian, scale invariance, Lorentz invariance, angular momentum quantization. Taught by Prof. Erich Poppitz

October 11, 2018 phy2403 , , , , , , , , , , , , , , ,

[Click here for a PDF of this post with nicer formatting (and a Mathematica listing that I didn’t include in this blog post’s latex export)]

DISCLAIMER: Very rough notes from class, with some additional side notes.

These are notes for the UofT course PHY2403H, Quantum Field Theory I, taught by Prof. Erich Poppitz fall 2018.

Last time

We followed a sequence of operations

  1. Noether’s theorem
  2. \( \rightarrow \) conserved currents
  3. \( \rightarrow \) charges (classical)
  4. \( \rightarrow \) “correspondence principle”
  5. \( \rightarrow \hatQ \)
  • Hermitian operators
  • “generators of symmetry”
    \begin{equation}\label{eqn:qftLecture9:20}
    \hatU(\alpha) = e^{i \alpha \hatQ}
    \end{equation}
    We found
    \begin{equation}\label{eqn:qftLecture9:40}
    \hatU(\alpha) \phihat \hatU^\dagger(\alpha) = \phihat + i \alpha \antisymmetric{\hatQ}{\phihat} + \cdots
    \end{equation}

Example: internal symmetries:

(non-spacetime), such as \( O(N)\) or \( U(1) \).

In QFT internal symmetries can have different “\underline{modes of realization}”.

    [I]

  1. “Wigner mode”. These are also called “unbroken symmetries”.
    \begin{equation}\label{eqn:qftLecture9:60}
    \hatQ \ket{0} = 0
    \end{equation}
    i.e. \( \hatU(\alpha) \ket{0} = 0 \).
    Ground state invariant. Formally \( :\hatQ: \) annihilates \( \ket{0} \).
    \( \antisymmetric{\hatQ}{\hatH} = 0 \) implies that all eigenstates are eigenstates of \( \hatQ \) in \( U(1) \). Example from HW 1
    \begin{equation}\label{eqn:qftLecture9:80}
    \hatQ = \text{“charge” under \( U(1) \)}.
    \end{equation}
    All states have definite charge, just live in QU.
  2. “Nambu-Goldstone mode” (Landau-ginsburg). This is also called a “spontaneously broken symmetry”\footnote{
    First encounter example (HWII, \( SU(2) \times SU(2) \rightarrow SU(2) \)). Here a \( U(1) \) spontaneous broken symmetry.}.
    \( H \) or \( L \) is invariant under symmetry, but ground state is not.

fig. 1. Mexican hat potential.

fig. 2. Degenerate Mexican hat potential ( v = 0)

Example:
\begin{equation}\label{eqn:qftLecture9:100}
\LL = \partial_\mu \phi^\conj \partial^\mu \phi – V(\Abs{\phi}),
\end{equation}
where
\begin{equation}\label{eqn:qftLecture9:120}
V(\Abs{\phi}) = m^2 \phi^\conj \phi + \frac{\lambda}{4} \lr{ \phi^\conj \phi }^2.
\end{equation}
When \( m^2 > 0 \) we have a Wigner mode, but when \( m^2 < 0 \) we have an issue: \( \phi = 0 \) is not a minimum of potential.
When \( m^2 < 0 \) we write
\begin{equation}\label{eqn:qftLecture9:140}
\begin{aligned}
V(\phi)
&= – m^2 \phi^\conj \phi + \frac{\lambda}{4} \lr{ \phi^\conj \phi}^2 \\
&= \frac{\lambda}{4} \lr{
\lr{ \phi^\conj \phi}^2 – \frac{4}{\lambda} m^2 } \\
&= \frac{\lambda}{4} \lr{
\phi^\conj \phi – \frac{2}{\lambda} m^2 }^2 – \frac{4 m^4}{\lambda^2},
\end{aligned}
\end{equation}
or simply
\begin{equation}\label{eqn:qftLecture9:780}
V(\phi)
=
\frac{\lambda}{4} \lr{ \phi^\conj \phi – v^2 }^2 + \text{const}.
\end{equation}
The potential (called the Mexican hat potential) is illustrated in fig. 1 for non-zero \( v \), and in
fig. 2 for \( v = 0 \).
We choose to expand around some point on the minimum ring (it doesn’t matter which one).
When there is no potential, we call the field massless (i.e. if we are in the minimum ring).
We expand as
\begin{equation}\label{eqn:qftLecture9:160}
\phi(x) = v \lr{ 1 + \frac{\rho(x)}{v} } e^{i \alpha(x)/v },
\end{equation}
so
\begin{equation}\label{eqn:qftLecture9:180}
\begin{aligned}
\frac{\lambda}{4}
\lr{\phi^\conj \phi – v^2}^2
&=
\lr{
v^2 \lr{ 1 + \frac{\rho(x)}{v} }^2
– v^2
}^2 \\
&=
\frac{\lambda}{4}
v^4 \lr{ \lr{ 1 + \frac{\rho(x)}{v} }^2 – 1 } \\
&=
\frac{\lambda}{4}
v^4
\lr{
\frac{2 \rho}{v} + \frac{\rho^2}{v^2}
}^2.
\end{aligned}
\end{equation}

\begin{equation}\label{eqn:qftLecture9:200}
\partial_\mu \phi =
\lr{
v \lr{ 1 + \frac{\rho(x)}{v} } \frac{i}{v} \partial_\mu \alpha
+ \partial_\mu \rho
} e^{i \alpha}
\end{equation}

so
\begin{equation}\label{eqn:qftLecture9:220}
\begin{aligned}
\LL
&= \Abs{\partial \phi^\conj}^2 – \frac{\lambda}{4} \lr{ \Abs{\phi^\conj}^2 – v^2 }^2 \\
&=
\partial_\mu \rho \partial^\mu \rho + \partial_\mu \alpha \partial^\mu \alpha \lr{ 1 + \frac{\rho}{v} }

\frac{\lambda v^4}{4} \frac{ 4\rho^2}{v^2} + O(\rho^3) \\
&=
\partial_\mu \rho \partial^\mu \rho
– \lambda v^2\rho^2
+
\partial_\mu \alpha \partial^\mu \alpha \lr{ 1 + \frac{\rho}{v} }.
\end{aligned}
\end{equation}
We have two fields, \( \rho \) : a massive scalar field, the “Higgs”, and a massless field \( \alpha \) (the Goldstone Boson).

\( U(1) \) symmetry acts on \( \phi(x) \rightarrow e^{i \omega } \phi(x) \) i.t.o \( \alpha(x) \rightarrow \alpha(x) + v \omega \).
\( U(1) \) global symmetry (broken) acts on the Goldstone field \( \alpha(x) \) by a constant shift. (\(U(1)\) is still a symmetry of the Lagrangian.)

The current of the \( U(1) \) symmetry is:
\begin{equation}\label{eqn:qftLecture9:240}
j_\mu = \partial_\mu \alpha \lr{ 1 + \text{higher dimensional \( \rho \) terms} }.
\end{equation}

When we quantize
\begin{equation}\label{eqn:qftLecture9:260}
\alpha(x) =
\int \frac{d^3p}{(2\pi)^3 \sqrt{ 2 \omega_p }} e^{i \omega_p t – i \Bp \cdot \Bx} \hat{a}_\Bp^\dagger +
\int \frac{d^3p}{(2\pi)^3 \sqrt{ 2 \omega_p }} e^{-i \omega_p t + i \Bp \cdot \Bx} \hat{a}_\Bp
\end{equation}
\begin{equation}\label{eqn:qftLecture9:280}
j^\mu(x) = \partial^\mu \alpha(x) =
\int \frac{d^3p}{(2\pi)^3 \sqrt{ 2 \omega_p }} \lr{ i \omega_\Bp – i \Bp } e^{i \omega_p t – i \Bp \cdot \Bx} \hat{a}_\Bp^\dagger +
\int \frac{d^3p}{(2\pi)^3 \sqrt{ 2 \omega_p }} \lr{ -i \omega_\Bp + i \Bp } e^{-i \omega_p t + i \Bp \cdot \Bx} \hat{a}_\Bp.
\end{equation}

\begin{equation}\label{eqn:qftLecture9:300}
j^\mu(x) \ket{0} \ne 0,
\end{equation}
instead it creates a single particle state.

Examples of symmetries

In particle physics, examples of Wigner vs Nambu-Goldstone, ignoring gravity the only exact internal symmetry in the standard module is
\( (B\# – L\#) \), believed to be a \( U(1) \) symmetry in Wigner mode.

Here \(B\#\) is the Baryon number, and \( L\# \) is the Lepton number. Examples:

  • \( B(p) = 1 \), proton.
  • \( B(q) = 1/3 \), quark
  • \( B(e) = 1 \), electron
  • \( B(n) = 1 \), neutron.
  • \( L(p) = 1 \), proton.
  • \( L(q) = 0 \), quark.
  • \( L(e) = 0 \), electron.

The major use of global internal symmetries in the standard model is as “approximate” ones. They become symmetries when one neglects some effect( “terms in \( \LL \)”).
There are other approximate symmetries (use of group theory to find the Balmer series).

Example from HW2:

QCD in limit
\begin{equation}\label{eqn:qftLecture9:320}
m_u = m_d = 0.
\end{equation}
\( m_u m_d \ll m_p \) (the products of the up-quark mass and the down-quark mass are much less than a composite one (name?)).
\( SU(2)_L \times SU(2)_R \rightarrow SU(2)_V \)

EWSB (Electro-Weak-Symmetry-Breaking) sector

When the couplings \( g_2, g_1 = 0 \). (\( g_2 \in SU(2), g_1 \in U(1) \)).

Scale invariance

\begin{equation}\label{eqn:qftLecture9:340}
\begin{aligned}
x &\rightarrow e^{\lambda} x \\
\phi &\rightarrow e^{-\lambda} \phi \\
A_\mu &\rightarrow e^{-\lambda} A_\mu
\end{aligned}
\end{equation}
Any unitary theory which is scale invariant is also \underline{conformal} invariant. Conformal invariance means that angles are preserved.
The point here is that there is more than scale invariance.

We have classical internal global continuous symmetries.
These can be either

  1. “unbroken” (Wigner mode)
    \begin{equation}\label{eqn:qftLecture9:360}
    \hatQ\ket{0} = 0.
    \end{equation}
  2. “spontaneously broken”
    \begin{equation}\label{eqn:qftLecture9:380}
    j^\mu(x) \ket{0} \ne 0
    \end{equation}
    (creates Goldstone modes).
  3. “anomalous”. Classical symmetries are not a symmetry of QFT.
    Examples:

    • Scale symmetry (to be studied in QFT II), although this is not truly internal.
    • In QCD again when \( \omega_\Bq = 0 \), a \( U(1\) symmetry (chiral symmetry) becomes exact, and cannot be preserved in QFT.
    • In the standard model (E.W sector), the Baryon number and Lepton numbers are not symmetries, but their difference \( B\# – L\# \) is a symmetry.

Lorentz invariance.

We’d like to study the action of Lorentz symmetries on quantum states. We are going to “go by the book”, finding symmetries, currents, quantize, find generators, and so forth.

Under a Lorentz transformation
\begin{equation}\label{eqn:qftLecture9:400}
x^\mu \rightarrow {x’}^\mu = {\Lambda^\mu}_\nu x^\nu,
\end{equation}
We are going to consider infinitesimal Lorentz transformations
\begin{equation}\label{eqn:qftLecture9:420}
{\Lambda^\mu}_\nu \approx
{\delta^\mu}_\nu + {\omega^\mu}_\nu
,
\end{equation}
where \( {\omega^\mu}_\nu \) is small.
A Lorentz transformation \( \Lambda \) must satisfy \( \Lambda^\T G \Lambda = G \), or
\begin{equation}\label{eqn:qftLecture9:800}
g_{\mu\nu} = {{\Lambda}^\alpha}_\mu g_{\alpha \beta} {{\Lambda}^\beta}_\nu,
\end{equation}
into which we insert the infinitesimal transformation representation
\begin{equation}\label{eqn:qftLecture9:820}
\begin{aligned}
0
&=
– g_{\mu\nu} +
\lr{ {\delta^\alpha}_\mu + {\omega^\alpha}_\mu }
g_{\alpha \beta}
\lr{ {\delta^\beta}_\nu + {\omega^\beta}_\nu } \\
&=
– g_{\mu\nu} +
\lr{
g_{\mu \beta}
+
\omega_{\beta\mu}
}
\lr{ {\delta^\beta}_\nu + {\omega^\beta}_\nu } \\
&=
– g_{\mu\nu} +
g_{\mu \nu}
+
\omega_{\nu\mu}
+
\omega_{\mu\nu}
+
\omega_{\beta\mu}
{\omega^\beta}_\nu.
\end{aligned}
\end{equation}
The quadratic term can be ignored, leaving just
\begin{equation}\label{eqn:qftLecture9:840}
0 =
\omega_{\nu\mu}
+
\omega_{\mu\nu},
\end{equation}
or
\begin{equation}\label{eqn:qftLecture9:860}
\omega_{\nu\mu} = – \omega_{\mu\nu}.
\end{equation}
Note that \( \omega \) is a completely antisymmetric tensor, and like \( F_{\mu\nu} \) this has only 6 elements.
This means that the
infinitesimal transformation of the coordinates is
\begin{equation}\label{eqn:qftLecture9:440}
x^\mu \rightarrow {x’}^\mu \approx x^\mu + \omega^{\mu\nu} x_\nu,
\end{equation}
the field transforms as
\begin{equation}\label{eqn:qftLecture9:460}
\phi(x) \rightarrow \phi'(x’) = \phi(x)
\end{equation}
or
\begin{equation}\label{eqn:qftLecture9:760}
\phi'(x^\mu + \omega^{\mu\nu} x_\nu) =
\phi'(x) + \omega^{\mu\nu} x_\nu \partial_\mu\phi(x) = \phi(x),
\end{equation}
so
\begin{equation}\label{eqn:qftLecture9:480}
\delta \phi = \phi'(x) – \phi(x) =
-\omega^{\mu\nu} x_\nu \partial_\mu \phi.
\end{equation}

Since \( \LL \) is a scalar
\begin{equation}\label{eqn:qftLecture9:500}
\begin{aligned}
\delta \LL
&=
-\omega^{\mu\nu} x_\nu \partial_\mu \LL \\
&=

\partial_\mu \lr{
\omega^{\mu\nu} x_\nu \LL
}
+
(\partial_\mu x_\nu) \omega^{\mu\nu} \LL \\
&=
\partial_\mu \lr{

\omega^{\mu\nu} x_\nu \LL
},
\end{aligned}
\end{equation}
since \( \partial_\nu x_\mu = g_{\nu\mu} \) is symmetric, and \( \omega \) is antisymmetric.
Our current is
\begin{equation}\label{eqn:qftLecture9:520}
J^\mu_\omega
=

\omega^{\mu\nu} x_\mu \LL
.
\end{equation}
Our Noether current is
\begin{equation}\label{eqn:qftLecture9:540}
\begin{aligned}
j^\nu_{\omega^{\mu\rho}}
&= \PD{\phi_{,\nu}}{\LL} \delta \phi – J^\mu_\omega \\
&=
\partial^\nu \phi\lr{ – \omega^{\mu\rho} x_\rho \partial_\mu \phi } + \omega^{\nu \rho} x_\rho \LL \\
&=
\omega^{\mu\rho}
\lr{
\partial^\nu \phi\lr{ – x_\rho \partial_\mu \phi } + {\delta^{\nu}}_\mu x_\rho \LL
} \\
&=
\omega^{\mu\rho} x_\rho
\lr{
-\partial^\nu \phi \partial_\mu \phi + {\delta^{\nu}}_\mu \LL
}
\end{aligned}
\end{equation}
We identify
\begin{equation}\label{eqn:qftLecture9:560}

{T^\nu}_\mu =
-\partial^\nu \phi \partial_\mu \phi + {\delta^{\nu}}_\mu \LL,
\end{equation}
so the current is
\begin{equation}\label{eqn:qftLecture9:580}
\begin{aligned}
j^\nu_{\omega_{\mu\rho}}
&=
-\omega^{\mu\rho} x_\rho
{T^\nu}_\mu \\
&=
-\omega_{\mu\rho} x^\rho
T^{\nu\mu}
.
\end{aligned}
\end{equation}
Define
\begin{equation}\label{eqn:qftLecture9:600}
j^{\nu\mu\rho} = \inv{2} \lr{ x^\rho T^{\nu\mu} – x^{\mu} T^{\nu\rho} },
\end{equation}
which retains the antisymmetry in \( \mu \rho \) yet still drops the parameter \( \omega^{\mu\rho} \).
To check that this makes sense, we can contract
\( j^{\nu\mu\rho} \) with \( \omega_{\rho\mu} \)
\begin{equation}\label{eqn:qftLecture9:880}
\begin{aligned}
j^{\nu\mu\rho} \omega_{\rho\mu}
&= -\inv{2} \lr{ x^\rho T^{\nu\mu} – x^{\mu} T^{\nu\rho} }
\omega_{\mu\rho} \\
&=
-\inv{2} x^\rho T^{\nu\mu}
\omega_{\mu\rho}
– \inv{2} x^{\mu} T^{\nu\rho}
\omega_{\rho\mu} \\
&=
-\inv{2} x^\rho T^{\nu\mu}
\omega_{\mu\rho}
– \inv{2} x^{\rho} T^{\nu\mu}
\omega_{\mu\rho} \\
&=
– x^{\rho} T^{\nu\mu}
\omega_{\mu\rho},
\end{aligned}
\end{equation}
which matches \ref{eqn:qftLecture9:580} as desired.

Example. Rotations \( \mu\rho = ij \)

\begin{equation}\label{eqn:qftLecture9:620}
\begin{aligned}
J^{0 i j} \epsilon_{ijk}
&=
\inv{2} \lr{ x^i T^{0j} – x^{j} T^{0i} } \epsilon_{ijk} \\
&=
x^i T^{0j} \epsilon_{ijk}.
\end{aligned}
\end{equation}
Observe that this has the structure of \( (\Bx \cross \Bp)_k \), where \( \Bp \) is the momentum density of the field.
Let
\begin{equation}\label{eqn:qftLecture9:640}
L_k \equiv Q_k = \int d^3 x J^{0ij} \epsilon_{ijk}.
\end{equation}
We can now quantize and build a generator
\begin{equation}\label{eqn:qftLecture9:660}
\begin{aligned}
\hatU(\Balpha)
&= e^{i \Balpha \cdot \hat{\BL}} \\
&= \exp\lr{i \alpha_k
\int d^3 x x^i \hat{T}^{0j} \epsilon_{ijk}
}
\end{aligned}
\end{equation}
From \ref{eqn:qftLecture9:560} we can quantize with \( T^{0j} = \partial^0 \phi \partial^j \phi \rightarrow \hat{\pi} \lr{\spacegrad \phihat}_j\), or
\begin{equation}\label{eqn:qftLecture9:900}
\begin{aligned}
\hatU(\Balpha)
&=
\exp\lr{i \alpha_k
\int d^3 x x^i \hat{\pi} (\spacegrad \phihat)_j \epsilon_{ijk}
} \\
&=
\exp\lr{i \Balpha \cdot
\int d^3 x \hat{\pi} \spacegrad \phihat \cross \Bx
}
\end{aligned}
\end{equation}
(up to a sign in the exponent which doesn’t matter)
\begin{equation}\label{eqn:qftLecture9:680}
\begin{aligned}
\phihat(\By) \rightarrow \hatU(\alpha) \phihat(\By) \hatU^\dagger(\alpha)
&\approx
\phihat(\By) +
i \Balpha \cdot
\antisymmetric{
\int d^3 x \hat{\pi}(\Bx) \spacegrad \phihat(\Bx) \cross \Bx
}
{
\phihat(\By)
} \\
&=
\phihat(\By) +
i \Balpha \cdot
\int d^3 x
(-i) \delta^3(\Bx – \By)
\spacegrad \phihat(\Bx) \cross \Bx \\
&=
\phihat(\By) +
\Balpha \cdot \lr{ \spacegrad \phihat(\By ) \cross \By}
\end{aligned}
\end{equation}
Explicitly, in coordinates, this is
\begin{equation}\label{eqn:qftLecture9:700}
\begin{aligned}
\phihat(\By)
&\rightarrow
\phihat(\By) +
\alpha^i
\lr{
\partial^j \phihat(\By) y^k \epsilon_{jki}
} \\
&=
\phihat(\By) –
\epsilon_{ikj} \alpha^i y^k \partial^j \phihat \\
&=
\phihat( y^j – \epsilon^{ikj} \alpha^i y^k ).
\end{aligned}
\end{equation}
This is a rotation. To illustrate, pick \( \Balpha = (0, 0, \alpha) \), so \( y^j \rightarrow y^j – \epsilon^{ikj} \alpha y^k \delta_{i3} = y^j – \epsilon^{3kj} \alpha y^k \), or
\begin{equation}\label{eqn:qftLecture9:n}
\begin{aligned}
y^1 &\rightarrow y^1 – \epsilon^{3k1} \alpha y^k = y^1 + \alpha y^2 \\
y^2 &\rightarrow y^2 – \epsilon^{3k2} \alpha y^k = y^2 – \alpha y^1 \\
y^3 &\rightarrow y^3 – \epsilon^{3k3} \alpha y^k = y^3,
\end{aligned}
\end{equation}
or in matrix form
\begin{equation}\label{eqn:qftLecture9:720}
\begin{bmatrix}
y^1 \\
y^2 \\
y^3 \\
\end{bmatrix}
\rightarrow
\begin{bmatrix}
1 & \alpha & 0 \\
-\alpha & 1 & 0 \\
0 & 0 & 1
\end{bmatrix}
\begin{bmatrix}
y^1 \\
y^2 \\
y^3 \\
\end{bmatrix}.
\end{equation}