{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "ff8465b6",
   "metadata": {},
   "source": [
    "| [**Lattice QCD I**](https://homepages.uni-regensburg.de/~lec17310/teaching/wise2122/lqft.html) | - | [**Christoph Lehner**](https://go.ur.de/lehner) |\n",
    "|:-|--|-:|\n",
    "\n",
    "# Chapter 5: symmetries of fundamental field theories\n",
    "\n",
    "In this chapter we explore how symmetries constrain the construction of an action for a theory of interest.  We start by discussing the spacetime symmetries."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1c216fb8",
   "metadata": {},
   "source": [
    "## The Lorentz group\n",
    "\n",
    "### Minkowski space\n",
    "Let us consider a photon moving with the speed of light $c$.  It shall\n",
    "propagate for a distance $d\\: x$ in an infinitesimal time $dt$, i.e.,\n",
    "$$\n",
    "\\begin{align}\n",
    "  c^2 dt^2 - d\\: x^2 = 0\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "If we consider a transformation of space and time coordinates $(t,x)$\n",
    "to $(t',x')$ the statement that the speed of light $c$ is the same in\n",
    "the new coordinate system is equivalent to the statement that also\n",
    "$$\n",
    "\\begin{align}\n",
    "  c^2 dt'^2 - d\\: x'^2 = 0\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "This property can now be expressed in a convenient mathematical\n",
    "representation by introducing vectors in a four-dimensional\n",
    "pseudo-Euclidean vector space with metric\n",
    "$$\n",
    "\\begin{align}\n",
    "  (g_{\\mu\\nu}) =\n",
    "  \\begin{pmatrix}\n",
    "    1 &    &    &    \\\\\n",
    "      & -1 &    &    \\\\\n",
    "      &    & -1 &    \\\\\n",
    "      &    &    & -1\n",
    "  \\end{pmatrix}.\n",
    "\\end{align}\n",
    "$$\n",
    "This vector space is called **Minkowski space**.  The infinitesimal\n",
    "time $dt$ and the corresponding vector $d\\: x$ are combined to a four\n",
    "vector $(dx^\\mu)=(c\\: dt,d\\: x)$ and the above equation can\n",
    "be written as\n",
    "$$\n",
    "\\begin{align}\n",
    "  ds^2 = dx_\\mu dx^\\mu = g_{\\mu\\nu} dx^\\mu dx^\\nu = 0\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "In this framework the transformations of coordinates that leave the\n",
    "speed of light invariant are just the isometries that leave the inner\n",
    "products of infinitesimal difference vectors invariant.\n",
    "\n",
    "For convenience we adopt natural units in the remainder of this chapter\n",
    "and set $c = 1$ (as we previously also had set $\\hbar = 1$)."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d55561f5",
   "metadata": {},
   "source": [
    "## Poincaré group\n",
    "The group of isometries of the Minkowski space is the **Poincaré\n",
    "  group** consisting of all transformations of the affine form\n",
    "  $$\n",
    "\\begin{align}\n",
    "  x'^\\mu = {\\Lambda^\\mu}_\\nu x^\\nu + T^\\mu\\,,\n",
    "\\end{align}\n",
    "$$\n",
    "with\n",
    "$$\n",
    "\\begin{align}\n",
    " dx'^\\mu dx'_\\mu = dx^\\mu dx_\\mu\n",
    "\\end{align}\n",
    "$$\n",
    "where $dx^\\mu$ is an infinitesimal difference vector in Minkowski\n",
    "spacetime, ${\\Lambda^\\mu}_\\nu$ is a real four-by-four matrix and\n",
    "$T^\\mu$ is a four vector describing translations.\n",
    "The above equation implies that\n",
    "$$\n",
    "\\begin{align}\n",
    "  dx'^\\mu dx'_\\mu = g_{\\mu\\nu} {\\Lambda^\\mu}_\\alpha{\\Lambda^\\nu}_\\beta\n",
    "  dx^\\alpha dx^\\beta = g_{\\alpha\\beta} dx^\\alpha dx^\\beta\n",
    "\\end{align}\n",
    "$$\n",
    "or\n",
    "$$\n",
    "\\begin{align}\n",
    "  g_{\\mu\\nu} {\\Lambda^\\mu}_\\alpha{\\Lambda^\\nu}_\\beta = g_{\\alpha\\beta}\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "This defining condition can be written in matrix form as\n",
    "$$\n",
    "\\begin{align}\\label{eqn:deflt}\n",
    "  \\Lambda^T g \\Lambda = g\n",
    "\\end{align}\n",
    "$$\n",
    "and thus\n",
    "$$\n",
    "\\begin{align}\\label{eqn:defltdet}\n",
    "  \\det \\Lambda = \\pm 1\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "The subgroup defined by $T^\\mu=0$, i.e., the subgroup of all linear\n",
    "transformations, is the **Lorentz group** and its elements are\n",
    "called **Lorentz transformations**."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "416c8245",
   "metadata": {},
   "source": [
    "## Restricted Lorentz group\n",
    "\n",
    "Let us first consider Lorentz transformations that are continuously\n",
    "connected to the identity transformation $\\Lambda = \\mathbb{1}$.  Lorentz\n",
    "transformations with this property live in a subgroup called the **restricted Lorentz group**.  Since a continuous transformation cannot\n",
    "change the sign of $\\det \\Lambda$, restricted Lorentz\n",
    "transformations have $\\det \\Lambda = 1$.\n",
    "\n",
    "A well-known subgroup of the restricted Lorentz group is the group of\n",
    "rotations with\n",
    "$$\n",
    "\\begin{align}\n",
    "  \\Lambda =\n",
    "  \\begin{pmatrix}\n",
    "    1 & 0 & 0 & 0 \\\\\n",
    "    0 &   &   &   \\\\\n",
    "    0 &   & R &   \\\\\n",
    "    0 &   &   &  \n",
    "  \\end{pmatrix},\n",
    "\\end{align}\n",
    "$$\n",
    "where $R^T R=\\mathbb{1}$, $\\det R = 1$.  Now by first rotating the spatial\n",
    "components appropriately we can restrict the remaining discussion to\n",
    "the two-dimensional subspace of vectors $(dx^\\mu)=(dt,dx,0,0)$.  The\n",
    "relevant Lorentz transformations are then of the form\n",
    "$$\n",
    "\\begin{align}\n",
    "  \\Lambda =\n",
    "  \\begin{pmatrix}\n",
    "    {\\Lambda^0}_0 & {\\Lambda^0}_1 & 0 & 0 \\\\\n",
    "    {\\Lambda^1}_0 & {\\Lambda^1}_1 & 0 & 0 \\\\\n",
    "    0 & 0 & 1 & 0 \\\\\n",
    "    0 & 0 & 0 & 1\n",
    "  \\end{pmatrix}.\n",
    "\\end{align}\n",
    "$$\n",
    "Thus for infinitesimal transformations $\\Lambda = \\mathbb{1} + G$ the defining\n",
    "condition of $\\det \\Lambda = 1 + \\mathrm{Tr}\\: G = \\pm 1$ yields $G^T g + g G = 0$, and\n",
    "therefore\n",
    "$$\n",
    "\\begin{align}\\notag\n",
    "  0 &=\n",
    "  \\begin{pmatrix}\n",
    "    {G^0}_0 & {G^1}_0 \\\\\n",
    "    {G^0}_1 & {G^1}_1\n",
    "  \\end{pmatrix}\n",
    "  \\begin{pmatrix}\n",
    "    1 & 0 \\\\\n",
    "    0 & -1\n",
    "  \\end{pmatrix}  \n",
    "  +\n",
    "  \\begin{pmatrix}\n",
    "    1 & 0 \\\\\n",
    "    0 & -1\n",
    "  \\end{pmatrix}  \n",
    "  \\begin{pmatrix}\n",
    "    {G^0}_0 & {G^0}_1 \\\\\n",
    "    {G^1}_0 & {G^1}_1\n",
    "  \\end{pmatrix} \\\\\n",
    "  &=\n",
    "  \\begin{pmatrix}\n",
    "    {G^0}_0 & -{G^1}_0 \\\\\n",
    "    {G^0}_1 & -{G^1}_1\n",
    "  \\end{pmatrix}\n",
    "  +\n",
    "  \\begin{pmatrix}\n",
    "    {G^0}_0 & {G^0}_1 \\\\\n",
    "    -{G^1}_0 & -{G^1}_1\n",
    "  \\end{pmatrix},\n",
    "\\end{align}\n",
    "$$\n",
    "or ${G^0}_0={G^1}_1=0$ and ${G^0}_1={G^1}_0$.  A finite transformation\n",
    "is thus given by"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "06f2445d",
   "metadata": {},
   "source": [
    "$$\n",
    "\\begin{align}\n",
    "  \\Lambda = \\exp\n",
    "  \\begin{pmatrix}\n",
    "    0 & s \\\\\n",
    "    s & 0\n",
    "  \\end{pmatrix}\n",
    "  =\\begin{pmatrix}\n",
    "    \\cosh s & \\sinh s \\\\\n",
    "    \\sinh s & \\cosh s\n",
    "  \\end{pmatrix}\n",
    "\\end{align}\n",
    "$$\n",
    "with arbitrary $s \\in \\mathbb{R}$.  Let us try to understand what the\n",
    "parameter $s$ means.  Consider an infinitesimal vector $(dt,dx)$ that\n",
    "transforms to\n",
    "$$\n",
    "\\begin{align}\n",
    "  \\begin{pmatrix}\n",
    "    dt^\\prime \\\\ \n",
    "    dx^\\prime\n",
    "  \\end{pmatrix}\n",
    "  &=\n",
    "  \\Lambda \n",
    "  \\begin{pmatrix}\n",
    "    dt \\\\ \n",
    "    dx\n",
    "  \\end{pmatrix} \\\\\n",
    "  &=\n",
    "  \\begin{pmatrix}\n",
    "    dt \\cosh s+dx \\sinh s \\\\ \n",
    "    dt \\sinh s+dx \\cosh s\n",
    "  \\end{pmatrix}\n",
    "  .\n",
    "\\end{align}\n",
    "$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8e780ed8",
   "metadata": {},
   "source": [
    "Now we define a transformed velocity \n",
    "$$\n",
    "\\begin{align}\n",
    "  v'=\\frac{dx'}{dt'}=\\frac{v \\cosh s+\\sinh s}{v \\sinh s+\\cosh s}\n",
    "\\end{align}\n",
    "$$\n",
    "with $v=dx/dt$.  If we have $v=0$ in the untransformed system we have\n",
    "$v'=\\tanh s$ in the transformed system.  Therefore transformations of\n",
    "this type describe a change of coordinates to a frame of reference\n",
    "that moves with a constant velocity of $\\tanh s$ relative to the\n",
    "original frame of reference.  These are the **boosts** in the\n",
    "special theory of relativity with **rapidity** $s$.\n",
    "\n",
    "Let us define $\\beta = \\tanh s$.  Since $\\cosh^2 s - \\sinh^2 s = 1$,\n",
    "we can show that\n",
    "$$\n",
    "\\begin{align}\n",
    "\\cosh s=\\frac1{\\sqrt{1-\\tanh^2 s}} = \\frac1{\\sqrt{1-\\beta^2}} = \\gamma\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "Therefore we can express the transformation also by the matrix\n",
    "$$\n",
    "\\begin{align}\n",
    "  \\Lambda(\\beta) = \\begin{pmatrix} \\gamma(\\beta) & \\gamma(\\beta)\n",
    "    \\beta \\\\ \\gamma(\\beta) \\beta & \\gamma(\\beta)\n",
    "  \\end{pmatrix}.\n",
    "\\end{align}\n",
    "$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f6a5574c",
   "metadata": {},
   "source": [
    "## Discrete Lorentz transformations\n",
    "Consider the vector $(x^\\mu)=(t,0)$ which is invariant under rotations\n",
    "and transforms to\n",
    "$$\n",
    "\\begin{align}\n",
    "  (x'^\\mu) =\n",
    "  \\begin{pmatrix}\n",
    "    t \\cosh s \\\\ \n",
    "    t \\sinh s\n",
    "  \\end{pmatrix}  \n",
    "\\end{align}\n",
    "$$\n",
    "under a boost with rapidity $s$.  Since $\\cosh s > 0$, we conclude\n",
    "that the sign of $x^0$ is invariant under boosts and thus under the\n",
    "complete restricted Lorentz group.\n",
    "\n",
    "Therefore, in order to obtain all possible Lorentz transformations,\n",
    "the discrete Lorentz transformation\n",
    "$$\n",
    "\\begin{align}\n",
    "  T =\n",
    "  \\begin{pmatrix}\n",
    "    -1 & & & \\\\\n",
    "    & 1 & & \\\\\n",
    "    & & 1 & \\\\\n",
    "    & & & 1\n",
    "  \\end{pmatrix}\n",
    "\\end{align}\n",
    "$$\n",
    "needs to be included in addition to restricted Lorentz\n",
    "transformations.  This is the **time reversal** operator.\n",
    "Furthermore the **space inversion** or **parity** operator\n",
    "$$\n",
    "\\begin{align}\n",
    "  P =\n",
    "  \\begin{pmatrix}\n",
    "    1 & & & \\\\\n",
    "    & -1 & & \\\\\n",
    "    & & -1 & \\\\\n",
    "    & & & -1\n",
    "  \\end{pmatrix}\n",
    "\\end{align}\n",
    "$$\n",
    "is also not a part of the restricted Lorentz group and needs to be\n",
    "included separately.\n",
    "\n",
    "The quotient group of the Lorentz group and the restricted Lorentz\n",
    "group is the discrete group with elements\n",
    "$$\n",
    "\\begin{align}\n",
    "  \\mathbb{1}, P, T, PT\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "In other words, the Lorentz group can be partitioned in four\n",
    "disconnected parts defined by\n",
    "$$\n",
    "\\begin{align}\n",
    "  \\det \\Lambda = \\pm 1\\,,\\qquad  \\textrm{Sign}\\: {\\Lambda^0}_0  = \\pm 1\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "We call transformations with $\\det \\Lambda = 1$ **proper** Lorentz\n",
    "transformations and transformations with $\\textrm{Sign}\\: {\\Lambda^0}_0 = 1$ **orthochronous**\n",
    "Lorentz transformations."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "0ec5e0ca",
   "metadata": {},
   "source": [
    "## Generators of the restricted Lorentz group\n",
    "\n",
    "Recall that infinitesimal restricted Lorentz transformations $\\Lambda\n",
    "= \\mathbb{1} + G$ satisfy\n",
    "$$\n",
    "\\begin{align}\\label{eqn:ltdefalg}\n",
    "  G^T g + g G = 0\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "We write $G$ in block form\n",
    "$$\n",
    "\\begin{align}\n",
    "  G =\n",
    "  \\begin{pmatrix}\n",
    "    G_{00} & G_{01} \\\\\n",
    "    G_{10} & G_{11}\n",
    "  \\end{pmatrix},\n",
    "\\end{align}\n",
    "$$\n",
    "where $G_{00}$ only acts on the temporal component, $G_{11}$ only acts\n",
    "on the spatial components, and $G_{01}$ and $G_{10}$ mix spatial and\n",
    "temporal components.  In this way the above equation can be\n",
    "expressed as\n",
    "$$\n",
    "\\begin{align}\\notag\n",
    "  0 &= \\begin{pmatrix}\n",
    "    1 & 0 \\\\\n",
    "    0 & -\\mathbb{1}_3\n",
    "  \\end{pmatrix}\n",
    "  \\begin{pmatrix}\n",
    "    G_{00} & G_{01} \\\\\n",
    "    G_{10} & G_{11}\n",
    "  \\end{pmatrix}\n",
    "  +\n",
    "  \\begin{pmatrix}\n",
    "    G_{00} & G^T_{10} \\\\\n",
    "    G^T_{01} & G^T_{11}\n",
    "  \\end{pmatrix}\n",
    "  \\begin{pmatrix}\n",
    "    1 & 0 \\\\\n",
    "    0 & -\\mathbb{1}_3\n",
    "  \\end{pmatrix} \\\\ &=\n",
    "  \\begin{pmatrix}\n",
    "    2 G_{00} & G_{01}-G^T_{10} \\\\\n",
    "    G^T_{01}- G_{10} & -G_{11}-G^T_{11}\n",
    "  \\end{pmatrix},\n",
    "\\end{align}\n",
    "$$\n",
    "where $\\mathbb{1}_3$ is the three-dimensional identity matrix.\n",
    "Therefore the defining conditions for generators of the\n",
    "restricted Lorentz group are\n",
    "$$\n",
    "\\begin{align}\n",
    " G_{01} = G^T_{10}\\,,\\qquad G_{00}=0\\,,\\qquad G_{11}^T=-G_{11}\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "This implies the following generators of the restricted Lorentz group.\n",
    "\n",
    "The boosts are generated by\n",
    "$$\n",
    "\\begin{align}\n",
    "  K_i =\n",
    "  \\begin{pmatrix}\n",
    "    0 & e_i^T \\\\\n",
    "    e_i & 0\n",
    "  \\end{pmatrix}\n",
    "\\end{align}\n",
    "$$\n",
    "with $(e_i)_j = \\delta_{ij}$ and $i=1,2,3$.  They satisfy \n",
    "$$\n",
    "\\begin{align}\\notag\n",
    "  [K_i,K_j] &=\n",
    "  \\begin{pmatrix}\n",
    "    e_i^T e_j - e_j^T e_i & 0 \\\\\n",
    "    0 & (e_i e_j^T - e_j e_i^T)_{ab}\n",
    "  \\end{pmatrix} =\n",
    "  \\begin{pmatrix}\n",
    "    0 & 0 \\\\\n",
    "    0 & (\\delta_{ia}\\delta_{jb} - \\delta_{ja}\\delta_{ib})_{ab}\n",
    "  \\end{pmatrix} \\\\&\n",
    "  =-\\varepsilon_{ijk} J_k\n",
    "\\end{align}\n",
    "$$\n",
    "with $J_k$ defined below.\n",
    "\n",
    "The rotations are generated by\n",
    "$$\n",
    "\\begin{align}\n",
    "  J_i =\n",
    "  \\begin{pmatrix}\n",
    "    0 & 0 \\\\\n",
    "    0 & L_i\n",
    "  \\end{pmatrix}\n",
    "\\end{align}\n",
    "$$\n",
    "with $(L_i)_{jk} = -\\varepsilon_{ijk}$ and $i=1,2,3$.  They satisfy\n",
    "$$\n",
    "\\begin{align}\\notag\n",
    "  [J_i,J_j] &= \\begin{pmatrix}\n",
    "    0 & 0 \\\\\n",
    "    0 & ([L_i,L_j])_{ab}\n",
    "  \\end{pmatrix}=\n",
    "  \\begin{pmatrix}\n",
    "    0 & 0 \\\\\n",
    "    0 & (\\varepsilon_{ial}\\varepsilon_{jlb} -\n",
    "    \\varepsilon_{jal}\\varepsilon_{ilb})_{ab}\n",
    "  \\end{pmatrix}\\\\&=\n",
    "  \\begin{pmatrix}\n",
    "    0 & 0 \\\\\n",
    "    0 & (\\delta_{ib}\\delta_{aj} - \\delta_{jb}\\delta_{ia})_{ab}\n",
    "  \\end{pmatrix}\n",
    "  =\\varepsilon_{ijk} J_k\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "Hence boosts do not form a subgroup of the restricted Lorentz group,\n",
    "but rotations do.  "
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c1e56496",
   "metadata": {},
   "source": [
    "Note that\n",
    "$$\n",
    "\\begin{align}\n",
    "  [K_i,J_j] &=  \n",
    "  \\begin{pmatrix}\n",
    "    0 & e_i^T \\\\\n",
    "    e_i & 0\n",
    "  \\end{pmatrix}\n",
    "  \\begin{pmatrix}\n",
    "    0 & 0 \\\\\n",
    "    0 & L_j\n",
    "  \\end{pmatrix}-\n",
    "  \\begin{pmatrix}\n",
    "    0 & 0 \\\\\n",
    "    0 & L_j\n",
    "  \\end{pmatrix}\n",
    "  \\begin{pmatrix}\n",
    "    0 & e_i^T \\\\\n",
    "    e_i & 0\n",
    "  \\end{pmatrix} \\\\\n",
    "  &=\n",
    "  \\begin{pmatrix}\n",
    "    0 & (e^T_i L_j)_a \\\\\n",
    "    (-L_j e_i)_a & 0\n",
    "  \\end{pmatrix} \\\\\n",
    "  &=\n",
    "  \\begin{pmatrix}\n",
    "    0 & \\varepsilon_{ija} \\\\\n",
    "    \\varepsilon_{ija} & 0\n",
    "  \\end{pmatrix}\n",
    "  =\\varepsilon_{ijk} K_k\n",
    "\\end{align}\n",
    "$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ac8ee8eb",
   "metadata": {},
   "source": [
    "The Lie algebra of the restricted Lorentz group is therefore given by\n",
    "$$\n",
    "\\begin{align}\n",
    "  [K_i,K_j] =-\\varepsilon_{ijk} J_k\\,,\\qquad\n",
    "  [J_i,J_j] = \\varepsilon_{ijk} J_k\\,,\\qquad\n",
    "  [K_i,J_j] = \\varepsilon_{ijk} K_k\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "A finite transformation is given by\n",
    "$$\n",
    "\\begin{align}\\label{eqn:ltfinite}\n",
    "  \\Lambda = \\exp[\\: s \\cdot \\: K + \\: \\phi \\cdot \\: J]\\,,\n",
    "\\end{align}\n",
    "$$\n",
    "where $\\: \\phi$ contains the angles of a rotation and $\\: s$ contains\n",
    "the rapidities of a boost.\n",
    "\n",
    "A convenient representation of the generators is given by\n",
    "$$\n",
    "\\begin{align}\n",
    "  S^\\pm_i = \\frac12(\\pm K_i + i J_i)\n",
    "\\end{align}\n",
    "$$\n",
    "with $(S^\\pm_i)^\\dagger = S_i^\\pm$ and $i=1,2,3$.  We find\n",
    "$$\n",
    "\\begin{align}\\notag\n",
    "  [S^a_i,S^b_j] &= (ab[K_i,K_j] + i b[J_i,K_j]+ i a[K_i,J_j]\n",
    "  -[J_i,J_j])/4 \\\\& = i\\varepsilon_{ijk}[i[(1+ab)/4] J_k + [(a+b)/4]\n",
    "  K_k] = \\delta_{ab} i\\varepsilon_{ijk} S^a_k\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "Therefore the group algebra factorizes in a direct product of two\n",
    "$\\textrm{SU}(2)$ algebras (this is of course not true in terms of groups).  We\n",
    "can express $J_i$ and $K_i$ in terms of $S_i^\\pm$ as\n",
    "$$\n",
    "\\begin{align}\n",
    "i J_i = S_i^+ + S_i^- \\,,\\qquad K_i = S_i^+ - S_i^-\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "Therefore $\\Lambda = \\exp[\\: s \\cdot \\: K + \\: \\phi \\cdot \\: J]$ can be written as\n",
    "$$\n",
    "\\begin{align}\\label{eqn:lts}\n",
    "  \\Lambda &= \\exp[s_i (S_i^+ - S_i^-) -i \\phi_i(S_i^+ + S_i^-)] =\n",
    "  \\exp[-i x_i S_i^+] \\exp[-i x_i^* S_i^-]\n",
    "\\end{align}\n",
    "$$\n",
    "with $x_i = \\phi_i + i s_i$."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5ac12023",
   "metadata": {},
   "source": [
    "## Translations in space and time\n",
    "The Casimir operators of $S^+$ and $S^-$ can now be used to classify\n",
    "the representations of the restricted Lorentz group.  These Casimir\n",
    "operators are, however, no invariants of representations of the\n",
    "complete Poincaré group since they do not commute with all\n",
    "translations of space and time.  In this section we show that the spin\n",
    "of a massive particle is, nevertheless, a well-defined quantity.\n",
    "\n",
    "We extend the Minkowski space by a fifth dimension so that we can\n",
    "express a general transformation of the Poincaré group conveniently as\n",
    "$$\n",
    "\\begin{align}\n",
    " x' = \\Gamma(\\Lambda,T) x\n",
    "\\end{align}\n",
    "$$\n",
    "with Lorentz transformation $\\Lambda$, a four-dimensional\n",
    "translation vector $(T_\\mu)$, $(x_\\mu)=(x_0,x_1,x_2,x_3,1)$, and\n",
    "$$\n",
    "\\begin{align}\n",
    "  \\Gamma(\\Lambda,T) =\n",
    "  \\begin{pmatrix}\n",
    "    \\Lambda & (T_\\mu) \\\\\n",
    "    0 & 1\n",
    "  \\end{pmatrix}\n",
    "\\end{align}\n",
    "$$\n",
    "in block notation.  The generators of translations in space and time\n",
    "$P_\\mu$ are therefore given by the matrices\n",
    "$$\n",
    "\\begin{align}\n",
    "  P_\\mu =\n",
    "  \\begin{pmatrix}\n",
    "    0 & (\\delta_{\\mu\\nu}) \\\\\n",
    "    0 & 0\n",
    "  \\end{pmatrix}\n",
    "\\end{align}\n",
    "$$\n",
    "in block notation.  A finite translation is given by\n",
    "$$\n",
    "\\begin{align}\n",
    "  \\Gamma(\\mathbb{1},T) = \\exp\\left[ \\:\\sum_{\\mu=0}^3 T_\\mu P_\\mu\\: \\right].\n",
    "\\end{align}\n",
    "$$\n",
    "We can now determine the algebra of the complete Poincaré group,\n",
    "$$\n",
    "\\begin{align}\\notag\n",
    "  [P_\\mu,P_\\nu] &= 0\\,, & [P_0,J_i] &= 0\\,, & [P_0,K_i] &= -P_i\\,, \\\\\n",
    "  [P_i,J_j] &= \\varepsilon_{ijk}P_k\\,, &  [P_i,K_j] &= -\\delta_{ij} P_0\\,, &  [K_i,K_j] &=-\\varepsilon_{ijk} J_k\\,, \\notag\\\\\n",
    "  [J_i,J_j] &= \\varepsilon_{ijk} J_k\\,, &  [K_i,J_j] &= \\varepsilon_{ijk} K_k\\,. &\n",
    "\\end{align}\n",
    "$$\n",
    "The Poincaré algebra has two Casimir operators.  The first one is\n",
    "given by\n",
    "$$\n",
    "\\begin{align}\n",
    "  C_1 = P_\\mu P^\\mu = P_0^2 - P_i^2\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "Note that for the matrix $P_\\mu$ given above this always vanishes, however,\n",
    "we can also consider different representations satisfying the same algebra with\n",
    "different $C_1$.\n",
    "We check explicitly that\n",
    "$$\n",
    "\\begin{align}\n",
    "  [P_\\mu,C_1] &= 0 \\,,\\\\\n",
    "  [J_i, C_1] &= [J_i, P_0^2] - [J_i,P_j^2] = -[J_i,P_j]P_j-P_j[J_i,P_j] \\notag\\\\\n",
    "&=2\\varepsilon_{ijk} P_k P_j = -2\\varepsilon_{ijk} P_k P_j = 0\\,,\\\\\n",
    "[K_i, C_1] &= [K_i,P_0]P_0 + P_0[K_i,P_0] -[K_i,P_j]P_j-P_j[K_i,P_j] \\notag\\\\\n",
    "&= 2 P_i P_0 -2P_i P_0 = 0\n",
    "\\end{align}\n",
    "$$\n",
    "for arbitrary $i$ and $\\mu$.  Let us pause at this point and ask what\n",
    "this means for a theory of a free particle with energy $E$ and\n",
    "momentum $\\: p$.  In quantum mechanics the generator of the\n",
    "translations in space, $P_i$, measures the $i$th component of the\n",
    "momentum, and the generator of the translations in time, $P_0$, measures\n",
    "the energy.  Therefore if we let $C_1$ act on a free particle state\n",
    "$\\vert{E,\\: p}\\rangle$ with energy $E$ and momentum $\\: p$ we find\n",
    "$$\n",
    "\\begin{align}\n",
    "  C_1 \\vert {E,\\: p} \\rangle = (E^2-\\: p^2) \\vert{E,\\: p} \\rangle = m^2 \\vert{E,\\: p} \\rangle\\,,\n",
    "\\end{align}\n",
    "$$\n",
    "where $m$ is the mass of the particle.  We can conclude that the mass\n",
    "of a particle is invariant under the Poincaré group and can be\n",
    "considered a well-defined property of a particle."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b6e96d7e",
   "metadata": {},
   "source": [
    "The second Casimir operator $C_2$ can be conveniently defined in terms\n",
    "of the **Pauli-Lubanski vector** $W^\\mu$ with\n",
    "$$\n",
    "\\begin{align}\n",
    "  W_0 = J_j P_j\\,,\\qquad\n",
    "  W_i = P_0 J_i - \\varepsilon_{ijk} K_j P_k\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "It is given by\n",
    "$$\n",
    "\\begin{align}\n",
    "  C_2 = W^\\mu W_\\mu = (W_0)^2 - (W_i)^2\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "In order to prove that $C_2$ is indeed a Casimir operator we first\n",
    "show that $W_\\mu$ commutes with translations, i.e.,\n",
    "$$\n",
    "\\begin{align}\n",
    "  [P_\\mu,W_0] &= [P_\\mu,J_j P_j] = [P_\\mu,J_j] P_j = (1-\\delta_{\\mu0}) \\varepsilon_{\\mu j k} P_k P_j = 0\\,,\\\\\n",
    "  [P_j,W_i] &= P_0 [P_j,J_i] - \\varepsilon_{ilk} [P_j,K_l] P_k = P_0 P_k(\\varepsilon_{jik} + \\varepsilon_{ijk}) = 0\\,,\\\\\n",
    "  [P_0,W_i] &= - \\varepsilon_{ilk} [P_0,K_l] P_k =\\varepsilon_{ilk} P_l P_k = 0\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "Next we discuss the commutators of $W_\\mu$ with boosts and calculate\n",
    "$$\n",
    "\\begin{align}\n",
    "  [K_j,W_0] &= [K_j,J_i P_i] = [K_j,J_i] P_i + J_i[K_j,P_i] \\notag\\\\&= \\varepsilon_{jik} K_k P_i + J_j P_0 = W_j\n",
    "\\end{align}\n",
    "$$\n",
    "and\n",
    "$$\n",
    "\\begin{align}\n",
    "  [K_j,W_i] &= [K_j,P_0 J_i] - \\varepsilon_{ilk} [K_j,K_l P_k]\n",
    "  \\notag\\\\&= P_0[K_j,J_i]+[K_j,P_0]J_i - \\varepsilon_{ilk} K_l[K_j,\n",
    "  P_k] - \\varepsilon_{ilk} [K_j,K_l] P_k \\notag\\\\&=\n",
    "  \\varepsilon_{jik}P_0 K_k+P_j J_i - \\varepsilon_{ilj} K_l P_0 +\n",
    "  \\varepsilon_{ilk} \\varepsilon_{jlr} J_r P_k \\notag\\\\&=\n",
    "  \\varepsilon_{jik}[P_0, K_k]+P_j J_i +\n",
    "  (\\delta_{ij}\\delta_{kr}-\\delta_{ir}\\delta_{kj}) J_r P_k \\notag\\\\&=\n",
    "  -\\varepsilon_{jik}P_k +[P_j, J_i] + \\delta_{ij}J_k P_k \\notag=\n",
    "  -\\varepsilon_{jik}P_k +\\varepsilon_{jik}P_k + \\delta_{ij}J_k P_k \\notag\\\\&=\n",
    "  \\delta_{ij} W_0\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "We finally calculate the commutators of $W_\\mu$ with rotations and find\n",
    "$$\n",
    "\\begin{align}\n",
    "  [J_j,W_0] &= [J_j,J_i P_i] = [J_j,J_i] P_i + J_i[J_j,P_i] =\n",
    "\\varepsilon_{jik} J_k P_i - \\varepsilon_{ijk} J_i P_k \\notag\\\\&=\n",
    "\\varepsilon_{jik} J_k P_i - \\varepsilon_{kji} J_k P_i = 0\\,,\\\\\n",
    "  [J_j,W_i] &=[J_j,P_0 J_i] - \\varepsilon_{ilk} [J_j,K_l P_k] \\notag\\\\&=\n",
    "  \\varepsilon_{jik} P_0 J_k - \\varepsilon_{ilk} [J_j,K_l]P_k - \\varepsilon_{ilk} K_l[J_j,P_k] \\notag\\\\&=\n",
    "  \\varepsilon_{jik} P_0 J_k + \\varepsilon_{ilk} \\varepsilon_{ljr}K_r P_k + \\varepsilon_{ilk} \\varepsilon_{kjr} K_l P_r \\notag\\\\&=\n",
    "  \\varepsilon_{jik} P_0 J_k + (\\varepsilon_{ilk} \\varepsilon_{ljr}+\\varepsilon_{ljk}\\varepsilon_{irl})K_r P_k \\notag\\\\&=\n",
    "  \\varepsilon_{jik} P_0 J_k + (\\delta_{kj}\\delta_{ir}-\\delta_{jr}\\delta_{ik})K_r P_k \\notag\\\\&=\n",
    "  \\varepsilon_{jik} P_0 J_k + \\varepsilon_{lji}\\varepsilon_{lkr} K_r P_k \\notag\\\\&=\n",
    "  \\varepsilon_{jik} (P_0 J_k + \\varepsilon_{klr} K_r P_l) = \\varepsilon_{jik} W_k\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "We observe that $W_\\mu$ has the same commutation relations with the\n",
    "other parts of the algebra as $P_\\mu$, and therefore $C_2$ is also a\n",
    "Casimir operator.\n",
    "\n",
    "For a massive particle we can calculate the action of $C_2$ in its\n",
    "rest frame, i.e.,\n",
    "$$\n",
    "\\begin{align}\n",
    "  C_2 \\vert{m,0}\\rangle = -m^2 J_i^2 \\vert{m,0}\\rangle\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "Therefore $\\vert{m,0}\\rangle$ must also be an eigenstate of $J_i^2$ and the\n",
    "corresponding eigenvalues $s(s+1)$ correspond to the **spin** or\n",
    "intrinsic rotation of the point-like particle.  In other words,\n",
    "massive particles can be classified according to their spin as defined\n",
    "by their behavior under the rotation group.\n",
    "\n",
    "For a massless particle there is no rest frame and thus the situation\n",
    "is more complicated.  It turns out that for massless particles the\n",
    "projection of the spin to the momentum,\n",
    "$$\n",
    "\\begin{align}\n",
    "  \\lambda = \\: J \\cdot \\hat P\\,,\n",
    "\\end{align}\n",
    "$$\n",
    "is a well-defined property and assumes the role of the spin of massive\n",
    "particles.  This property is called **helicity**.\n",
    "\n",
    "For a detailed discussion of the representation theory of the complete\n",
    "Poincaré group, see, e.g., books of Weinberg or Ryder."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "2fb82a71",
   "metadata": {},
   "source": [
    "In the following discussion, we construct a Lagrangian of massive spin $0$ and spin $1/2$\n",
    "particles that are invariant under orthochronous Lorentz\n",
    "transformations.\n",
    "\n",
    "## Spinor representations\n",
    "Note that the sub-sectors $+$ and $-$ of the restricted Lorentz group\n",
    "both transform identically under rotations with\n",
    "$$\n",
    "\\begin{align}\n",
    "  \\Lambda &= \\exp[-i\\phi_i S_i^\\pm]\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "Since $\\phi_i$ are the angles of a rotation in space and the $S_i$\n",
    "span the algebra of $\\textrm{SU}(2)$ the different representations of $S$\n",
    "correspond to different spin states.  Possible representations of\n",
    "$S^+ \\oplus S^-$ are\n",
    "$$\n",
    "\\begin{align}\n",
    "  0 \\oplus 0\\,,\\qquad \\frac12 \\oplus 0\\,,\\qquad 0 \\oplus \\frac12\\,,\\qquad\n",
    "  \\frac12 \\oplus \\frac12\\,,\\qquad \\ldots\\,.\n",
    "\\end{align}\n",
    "$$\n",
    "The\n",
    "$$\n",
    "0 \\oplus 0\n",
    "$$\n",
    "representation is the singlet representation, i.e., fields $\\phi$ transforming in this representation do not change under a restricted Lorentz transformation.  Such fields therefore correspond to spin $0$ particles.\n",
    "In the next lecture, we complete the construction of a minimal Lagrangian for a spin $0$ particle and extend the discussion to spin $1/2$ particles as well."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "02bb6965",
   "metadata": {},
   "source": [
    "## The Lagrangian of spin 0 fields\n",
    "\n",
    "We first consider the $0 \\oplus 0$ representation of $S^+ \\oplus S^-$ for which fields\n",
    "$\\phi$ transform as\n",
    "$$\n",
    "\\begin{align}\n",
    "\\phi \\to \\phi' = \\phi\n",
    "\\end{align}\n",
    "$$\n",
    "under the restricted Lorentz group.  Additionally, the derivatives\n",
    "transform as\n",
    "$$\n",
    "\\begin{align}\n",
    "\\partial_\\mu \\to \\partial_\\mu^\\prime  &= {\\Lambda_{\\mu}}^\\nu \\partial_\\nu  \\\\\n",
    "\\partial^\\mu \\to {\\partial^\\mu}^\\prime &= {\\Lambda^{\\mu}}_\\nu \\partial^\\nu\n",
    "\\end{align}\n",
    "$$\n",
    "such that\n",
    "$$\n",
    "\\begin{align}\n",
    "\\partial_\\mu \\partial^\\mu \\to \\partial_\\alpha \\partial^\\beta {\\Lambda_{\\mu}}^\\alpha {\\Lambda^{\\mu}}_\\beta =\\partial_\\alpha \\partial^\\beta {g^{\\alpha}}_\\beta = \\partial_\\mu \\partial^\\mu \\,.\n",
    "\\end{align}\n",
    "$$\n",
    "We also remind the reader that since $\\vert \\det \\Lambda \\vert = 1$, the integral measure $d^4x$ is invariant under the restricted Lorentz group.\n",
    "\n",
    "The action of chapter 4\n",
    "$$\n",
    "\\begin{align}\n",
    " S = \\int d^4x {\\cal L}(\\phi(x),\\partial_\\mu \\phi(x)) = \\int d^4x \\frac12 \\phi(x) (-\\partial_\\mu \\partial^\\mu - m^2) \\phi(x)\n",
    "\\end{align}\n",
    "$$\n",
    "is therefore invariant under the restricted Lorentz group.  The factor of $1/2$ and for that matter any factor of the kinetic term (derivative term) can be absorbed in a re-definition of the fields, over which we will integrate.  The only relevant quantity here is the relative coefficient of the $\\phi^2$ mass term and the kinetic term.\n",
    "\n",
    "### Linear terms in the field are redundant\n",
    "Consider shifting the fields by a constant, i.e., $\\phi(x) \\to \\phi_0 + \\phi(x)$ then\n",
    "$$\n",
    "\\begin{align}\n",
    " S \\to -\\int d^4x \\frac12 \\phi_0^2 m^2 - \\int d^4x \\frac12 \\phi(x) (2m^2 \\phi_0) - \\int d^4x \\frac12 \\phi_0 \\partial_\\mu \\partial^\\mu \\phi(x) + \\int d^4x \\frac12 \\phi(x) (-\\partial_\\mu \\partial^\\mu - m^2) \\phi(x)\n",
    "\\end{align}\n",
    "$$\n",
    "so this transformation generates a linear term in $\\phi(x)$, a constant term, and a term proportional to\n",
    "$$\n",
    "\\begin{align}\n",
    "\\int d^4x \\partial_\\mu \\partial^\\mu \\phi(x) \\,.\n",
    "\\end{align}\n",
    "$$\n",
    "Such a term, however, vanishes as long as we assume that $\\phi$ vanishes at sufficiently long distances or we consider a path integral in which the fields live on a torus.  Therefore linear terms in $\\phi$ in the action can be absorbed by a constant shift in fields in the path integral and are therefore redundant.\n",
    "\n",
    "### Action dimensionality and higher-dimensional contributions\n",
    "Finally, we note that the action must have mass dimension 0.  This can be quickly verified for the 1d free particle case\n",
    "$$\n",
    "\\begin{align}\n",
    "[S] = [ \\int dt p^2/2/m ] = [t] + 2[p] - [m] = -1 + 2 - 1 = 0 \\,.\n",
    "\\end{align}\n",
    "$$\n",
    "Because of this, any additional terms in the action with higher mass dimensions must also be accompanied with additional dimensionful constant prefactors.  We will investigate the role of such terms when we discuss the continuum limit of field theories.\n",
    "\n",
    "### Klein-Gordon equation\n",
    "\n",
    "The equations of motion\n",
    "$$\n",
    "\\begin{align}\n",
    "\\frac{\\partial {\\cal L}}{\\partial \\phi} - \\partial_\\mu \\frac{\\partial {\\cal L}}{\\partial (\\partial_\\mu \\phi)}\n",
    "&=-(m^2 + \\partial_\\mu \\partial^\\mu) \\phi = 0 \\,.\n",
    "\\end{align}\n",
    "$$\n",
    "yield the Klein-Gordon equations of relativistic scalar particles\n",
    "$$\n",
    "\\begin{align}\n",
    "(m^2 + \\partial_\\mu \\partial^\\mu) \\phi = 0 \\,.\n",
    "\\end{align}\n",
    "$$"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "572d0c1d",
   "metadata": {},
   "source": [
    "## The Lagrangian of spin $1/2$ fields\n",
    "\n",
    "We now try to construct a theory of spin $1/2$ particles such as electrons, neutrons or protons.  To this end we\n",
    "first consider fields that transform in the $(1/2) \\oplus 0$\n",
    "representations of $S^+ \\oplus S^-$.  We set\n",
    "\\begin{align}\\label{eqn:spinansatz}\n",
    "  S^+_i = \\frac12 \\sigma_i\\,,\\qquad S^-_i = 0\n",
    "\\end{align}\n",
    "with Pauli matrices $\\sigma_i$ and consider two-dimensional spinors\n",
    "$\\psi^+$ which transform as\n",
    "\\begin{align}\n",
    "  \\psi'^+ = \\exp[(s_i -i \\phi_i)\\sigma_i/2]\\psi^+\n",
    "\\end{align}\n",
    "under the restricted Lorentz group.\n",
    "\n",
    "Let us try to construct a Lagrangian with fields $\\psi^+$.  Each term\n",
    "in the Lagrangian has to satisfy the following properties: (i) Due to\n",
    "relativity each term has to be a Lorentz scalar.  (ii) The Lagrangian\n",
    "has to be real.  (iii) Each term has to have mass dimension of $4$\n",
    "(the action has to be dimensionless).\n",
    "\n",
    "### The mass term\n",
    "One may be tempted to write down a simple mass term of the form\n",
    "\\begin{align}\\label{eqn:simplemass}\n",
    " {\\cal L}_\\text{mass} = m (\\psi^+)^\\dagger \\psi^+\\,.\n",
    "\\end{align}\n",
    "Unfortunately, such a term does not satisfy (i) and is therefore not\n",
    "allowed in the Lagrangian.  We discuss how a proper mass term can be\n",
    "constructed if we consider the representation $(1/2) \\oplus (1/2)$ of\n",
    "$S^+ \\oplus S^-$ below.  This is the mass term relevant for QCD.  It\n",
    "is, however, instructive to consider another way to construct an\n",
    "invariant mass term that involves only $(\\psi^+)^T$ and $\\psi^+$, the\n",
    "**Majorana** mass term $(\\psi^+)^T \\sigma_2 \\psi^+$.  First note\n",
    "that $\\sigma_i^T = \\sigma_i (-1)^{\\delta_{i2}}$ with anticommutator\n",
    "$\\{\\sigma_i,\\sigma_j\\}=2\\delta_{ij}\\mathbb{1}$, and therefore\n",
    "\\begin{align}\\label{eqn:sigmatwo}\n",
    "  \\sigma_i^T \\sigma_2 \\sigma_i =\n",
    "  (-1)^{\\delta_{i2}}\\sigma_i\\sigma_2\\sigma_i =\n",
    "  (-1)^{\\delta_{i2}}(-\\sigma_2+2\\delta_{i2}\\sigma_i)=-\\sigma_2\\,,\n",
    "\\end{align}\n",
    "where no sum over $i$ is implied.  Thus $\\sigma_i^T \\sigma_2=-\\sigma_2\n",
    "\\sigma_i$, and for infinitesimal transformations with coordinates $x_i\n",
    "\\ll 1$ we find\n",
    "\\begin{align}\\notag\n",
    "  (\\psi^+)^T \\sigma_2 \\psi^+ \\to &(\\psi^+)^T (\\mathbb{1} -i x_i \\sigma_i^T/2)\\sigma_2\n",
    "  (\\mathbb{1} -i x_i \\sigma_i/2)\\psi^+ \\\\&= (\\psi^+)^T \\sigma^2 \\psi^+ -(i/2) x_i\n",
    "  (\\psi^+)^T(\\sigma_i^T\\sigma^2+\\sigma^2\\sigma_i)\\psi^+ \\notag\\\\&= (\\psi^+)^T\n",
    "  \\sigma^2 \\psi^+\\,.\n",
    "\\end{align}\n",
    "In order to make this term real we need to also include its complex\n",
    "conjugate.  Since $\\sigma_2$ is purely imaginary we write\n",
    "\\begin{align}\n",
    "  {\\cal L}_{\\text{Majorana mass}} = im ((\\psi^+)^T \\sigma^2 \\psi^+ -\n",
    "  (\\psi^+)^\\dagger \\sigma^2 (\\psi^+)^*)\\,.\n",
    "\\end{align}\n",
    "Note that for a two-component field $(\\psi^+)^T=(a,b)$ we find\n",
    "$(\\psi^+)^T \\sigma^2 \\psi^+=i(b a - a b)$.  Therefore if we consider\n",
    "$a$ and $b$ to be ordinary numbers, the Majorana mass term would\n",
    "vanish identically.  However, in a quantized theory we will later show that $a$ and $b$\n",
    "anticommute since they correspond to fermions, and the Majorana mass\n",
    "term is nonzero."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "728eaec0",
   "metadata": {},
   "source": [
    "### The kinetic term\n",
    "In this subsection we consider terms of the form\n",
    "\\begin{align}\n",
    "  (\\psi^+)^\\dagger R \\psi^+\\,,\n",
    "\\end{align}\n",
    "where $R$ contains objects that transform non-trivially under the\n",
    "restricted Lorentz group.  We use the first non-trivial ansatz including Lorentz vectors\n",
    "\\begin{align}\n",
    "  R = M_\\mu v^\\mu\\,,\n",
    "\\end{align}\n",
    "where $v^\\mu$ is a contravariant vector, $M_\\mu$ is a matrix in the\n",
    "two-dimensional spin space and the sum over $\\mu$ is implied.  Note\n",
    "that $M_\\mu$ is not a Lorentz vector.  Therefore under Lorentz\n",
    "transformations we find\n",
    "\\begin{align}\n",
    "  R' = M_\\mu v'^\\mu = M_\\mu {\\Lambda^\\mu}_\\nu v^\\nu\\,.\n",
    "\\end{align}\n",
    "In order to construct an invariant term we need\n",
    "\\begin{align}\\label{eqn:definvarkin}\n",
    "  (\\psi^+)^\\dagger R \\psi^+ = (\\psi^+)^\\dagger\\exp[(s_i +i \\phi_i)\\sigma_i/2] R' \\exp[(s_i -i \\phi_i)\\sigma_i/2]\\psi^+\\,.\n",
    "\\end{align}\n",
    "Let us first consider a infinitesimal boost in $r$ direction, i.e.,\n",
    "$\\: \\phi=0$, $s_i=\\delta_{ir} s$ with $s \\ll 1$, and\n",
    "\\begin{align}\n",
    "  v'^\\mu = v^\\mu + s {K_r^\\mu}_\\nu v^\\nu\\,.\n",
    "\\end{align}\n",
    "Now Eq. \\eqref{eqn:definvarkin} gives\n",
    "\\begin{align}\n",
    "  M_\\mu v^\\mu&\\stackrel!= [\\mathbb{1} + s \\sigma_r/2] M_\\mu v'^\\mu [\\mathbb{1} + s\n",
    "  \\sigma_r/2] \\notag\\\\&=[\\mathbb{1} + s \\sigma_r/2] M_\\mu [v^\\mu + s {K_r^\\mu}_\\nu\n",
    "  v^\\nu] [\\mathbb{1} + s \\sigma_r/2] \\notag\\\\& =M_\\mu v^\\mu + s (M_\\mu\n",
    "  {K_r^\\mu}_\\nu v^\\nu + \\sigma_r M_\\mu v^\\mu/2 + v^\\mu M_\\mu\n",
    "  \\sigma_r/2)\\,.\n",
    "\\end{align}\n",
    "This has to hold for all $v^\\mu$ so that we need\n",
    "\\begin{align}\n",
    "  0= M_\\nu {K_r^\\nu}_\\mu + \\{\\sigma_r, M_\\mu\\} /2 = M_0\n",
    "  \\delta_{r\\mu} + \\delta_{\\mu0} M_r + \\{\\sigma_r, M_\\mu\\} /2\\,.\n",
    "\\end{align}\n",
    "Now this means that\n",
    "\\begin{align}\\label{eqn:defclt}\n",
    "M_r = -\\{\\sigma_r, M_0\\} /2  \\,,\\qquad M_0 \\delta_{ri} =- \\{\\sigma_r, M_i\\} /2\\,.\n",
    "\\end{align}\n",
    "Next, we consider a rotation about the $r$ axis, i.e.,\n",
    "$\\: s=0$, $\\phi_i=\\delta_{ir} \\phi$ with $\\phi \\ll 1$, and\n",
    "\\begin{align}\n",
    "  v'^\\mu = v^\\mu + \\phi {J_r^\\mu}_\\nu v^\\nu\\,.\n",
    "\\end{align}\n",
    "Now Eq. \\eqref{eqn:definvarkin} gives\n",
    "\\begin{align}\\notag\n",
    "  M_\\mu v^\\mu&\\stackrel!= M_\\mu v'^\\mu =[\\mathbb{1} + i \\phi \\sigma_r/2]\n",
    "  M_\\mu [v^\\mu + \\phi {J_r^\\mu}_\\nu v^\\nu] [\\mathbb{1} - i\\phi \\sigma_r/2]\n",
    "  \\\\& =M_\\mu v^\\mu + \\phi (M_\\nu {J_r^\\nu}_\\mu v^\\mu + i\n",
    "  [\\sigma_r/2, M_\\mu] v^\\mu)\\,.\n",
    "\\end{align}\n",
    "This has to hold for all $v^\\mu$ so that we need\n",
    "\\begin{align}\n",
    "  0 = M_\\nu {J_r^\\nu}_\\mu + i [\\sigma_r/2, M_\\mu] =\n",
    "  -\\varepsilon_{\\mu ri} (1-\\delta_{\\mu0})M_i + i [\\sigma_r/2,\n",
    "  M_\\mu]\\,,\n",
    "\\end{align}\n",
    "and thus\n",
    "\\begin{align}\n",
    "  [\\sigma_r, M_\\mu] = i2\\varepsilon_{r\\mu i} (1-\\delta_{\\mu0})M_i\\,.\n",
    "\\end{align}\n",
    "For $\\mu=0$ this means that $[\\sigma_r, M_0] = 0$ for arbitrary\n",
    "$r$.  This is only satisfied for \n",
    "\\begin{align}\n",
    "  M_0 = c \\mathbb{1}\\,.\n",
    "\\end{align}\n",
    "For $\\mu=j$ with $j=1,2,3$ this means that \n",
    "\\begin{align}\n",
    "  [\\sigma_r, M_j] = i2\\varepsilon_{rji}M_i\\,.\n",
    "\\end{align}\n",
    "We know that this is satisfied by the Pauli matrices \n",
    "\\begin{align}\n",
    "  M_j = \\sigma_j\\,.\n",
    "\\end{align}\n",
    "We determine $c$ from Eq. \\eqref{eqn:defclt} and $\\{\\sigma_r,\n",
    "\\sigma_i\\} = 2\\mathbb{1} \\delta_{ri}$ and find $c=-1$.  It is easy to check\n",
    "that if we would have considered the sector $-$ instead of $+$ the\n",
    "solution would be $c=1$.  We define\n",
    "$(M_\\mu)=(\\sigma^+_\\mu)=(-\\mathbb{1},\\sigma_1,\\sigma_2,\\sigma_3)$ \n",
    "and $(\\sigma_+^\\mu)=(-\\mathbb{1},-\\sigma_1,-\\sigma_2,-\\sigma_3)$ so that\n",
    "\\begin{align}\n",
    "  (\\psi^+)^\\dagger \\sigma^+_\\mu v^\\mu \\psi^+ = (\\psi^+)^\\dagger \\sigma^+_\\nu g^{\\mu\\nu} v_\\mu \\psi^+ = (\\psi^+)^\\dagger \\sigma_+^\\mu v_\\mu \\psi^+\n",
    "\\end{align}\n",
    "is invariant under the restricted Lorentz group.  While $\\sigma_+^\\mu$\n",
    "does not transform as a Lorentz vector, we can conclude that\n",
    "$$\n",
    "\\begin{align}\\label{eqn:lorentzvector}\n",
    "  (\\psi^+)^\\dagger \\sigma_+^\\mu \\psi^+ \n",
    "\\end{align}\n",
    "$$\n",
    "does transform as a Lorentz vector.  Note that the relevant matrices\n",
    "for the $-$ sector are\n",
    "$(\\sigma^-_\\mu)=(\\mathbb{1},\\sigma_1,\\sigma_2,\\sigma_3)$ and\n",
    "$(\\sigma_-^\\mu)=(\\mathbb{1},-\\sigma_1,-\\sigma_2,-\\sigma_3)$.\n",
    "\n",
    "By substituting $v_\\mu = \\partial_\\mu$ we can thus construct an\n",
    "invariant kinetic term that only involves $+$ fields."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5b1e797d",
   "metadata": {},
   "source": [
    "### Chirality\n",
    "Let us consider all orthochronous Lorentz transformations, i.e, let us\n",
    "include the parity operator in addition to the restricted Lorentz\n",
    "transformations.  The action of parity is defined by\n",
    "\\begin{align}\\label{eqn:defparityop}\n",
    "  \\Lambda(s,\\phi) P = P^2 \\Lambda(s,\\phi) P = P \\Lambda(-s,\\phi)\n",
    "\\end{align}\n",
    "due to $P^2=1$, $P K_i P = - K_i$, and $P J_i P = J_i$.\n",
    "Equation \\eqref{eqn:defparityop} has to hold for all representations,\n",
    "and therefore the action $D(\\Lambda)$ of Lorentz transformations\n",
    "$\\Lambda$ on $\\psi_+$ yields\n",
    "$$\n",
    "\\begin{align}\\label{eqn:parityact}\n",
    "  (D(P) \\psi'_+) &= D(P) D(\\Lambda(s,\\phi))\\psi_+  =\n",
    "  D(\\Lambda(-s,\\phi)) (D(P) \\psi_+)\n",
    "\\end{align}\n",
    "$$\n",
    "with $\\psi'_+ = D(\\Lambda(s,\\phi))\\psi_+$.  We observe that the field\n",
    "$D(P)\\psi^+$ transforms according to the $1/2$ representation of\n",
    "$S^-$.  Therefore if we want to construct a theory that is also\n",
    "invariant under parity, we need to include a spin $1/2$ representation\n",
    "of $S^-$ as well.  The twofold structure that emerges from the $(1/2)\n",
    "\\oplus (1/2)$ representation of $S^+ \\oplus S^-$ is called **chirality**.\n",
    "\n",
    "We consider a spinor\n",
    "\\begin{align}\n",
    "  \\psi =\n",
    "  \\begin{pmatrix}\n",
    "    \\psi_- \\\\\n",
    "    \\psi_+\n",
    "  \\end{pmatrix},\n",
    "\\end{align}\n",
    "where $\\psi_\\pm$ transform according to the $1/2$ representation of\n",
    "$S^\\pm$.  The action of parity shall be given by\n",
    "\\begin{align}\n",
    "  D(P)\\psi =\n",
    "  \\begin{pmatrix}\n",
    "    \\psi_+ \\\\\n",
    "    \\psi_-\n",
    "  \\end{pmatrix},\n",
    "\\end{align}\n",
    "in accordance with Eq. \\eqref{eqn:parityact}.  We can write down a\n",
    "mass term\n",
    "\\begin{align}\n",
    "  {\\cal L}_{mass} &= m \\bar\\psi \\psi\n",
    "\\end{align}\n",
    "with\n",
    "\\begin{align}\n",
    "  \\bar\\psi = \\begin{pmatrix}\n",
    "    \\psi_+^\\dagger & \\psi_-^\\dagger\n",
    "  \\end{pmatrix}\n",
    "\\end{align}\n",
    "that is invariant under orthochronous Lorentz transformations, see Eq. \\eqref{eqn:lts}, repeated here:\n",
    "$$\n",
    "\\begin{align}\n",
    "  \\Lambda &=\n",
    "  \\exp[-i x_i S_i^+] \\exp[-i x_i^* S_i^-] \\,.\n",
    "\\end{align}\n",
    "$$\n",
    "\n",
    "We already know that\n",
    "\\begin{align}\n",
    "  (\\psi^+)^\\dagger \\sigma_+^\\mu \\partial_\\mu \\psi^+\n",
    "\\end{align}\n",
    "and\n",
    "\\begin{align}\n",
    "  (\\psi^-)^\\dagger \\sigma_-^\\mu \\partial_\\mu \\psi^-\n",
    "\\end{align}\n",
    "are both invariant under the restricted Lorentz group.  Under\n",
    "parity we have $\\psi^+ \\leftrightarrow \\psi^-$ and $\\partial_i \\to\n",
    "-\\partial_i$ for $i=1,2,3$ so that\n",
    "\\begin{align}\n",
    "  \\sigma_+^\\mu \\partial_\\mu \\leftrightarrow -\\sigma_-^\\mu \\partial_\\mu\\,.\n",
    "\\end{align}\n",
    "Therefore we can construct a real and Lorentz invariant kinetic term\n",
    "\\begin{align}\n",
    "  {\\cal L}_{kinetic} &= i[(\\psi^-)^\\dagger \\sigma_-^\\mu \\partial_\\mu\n",
    "  \\psi^- - (\\psi^+)^\\dagger \\sigma_+^\\mu \\partial_\\mu\n",
    "  \\psi^+] \\notag\\\\\n",
    "  &=\\bar\\psi i \\gamma^\\mu \\partial_\\mu \\psi\n",
    "\\end{align}\n",
    "with\n",
    "\\begin{align}\n",
    "  \\gamma^\\mu =\n",
    "  \\begin{pmatrix}\n",
    "    0 & -\\sigma^\\mu_+ \\\\ \\sigma^\\mu_- & 0\n",
    "  \\end{pmatrix}.\n",
    "\\end{align}\n",
    "The factor $i$ is needed since $\\partial_\\mu$ is anti-Hermitian, i.e.,\n",
    "\\begin{align}\\label{eqn:antihermitianpartial}\n",
    "  \\langle{\\psi'}\\vert \\partial_x \\vert{\\psi}\\rangle &= \\int dx \\:\n",
    "  \\psi'^*(x) \\partial_x \\psi(x) =- \\int dx\n",
    "  (\\partial_x\\psi'^*(x)) \\psi(x) \\notag\\\\&\n",
    "  =-\\langle{\\psi}\\vert \\partial_\\mu \\vert{\\psi'}\\rangle^*\n",
    "\\end{align}\n",
    "for arbitrary fields $\\psi$ and $\\psi'$ with vanishing spacetime\n",
    "boundary contributions.\n",
    "\n",
    "We write out the gamma matrices as\n",
    "\\begin{align}\\label{eqn:defgamma}\n",
    "  \\gamma^0 = \\begin{pmatrix}\n",
    "    0 & \\mathbb{1} \\\\ \\mathbb{1} & 0\n",
    "  \\end{pmatrix},\\qquad \\gamma^i = \\begin{pmatrix} 0 & \\sigma_i \\\\\n",
    "    -\\sigma_i & 0\n",
    "  \\end{pmatrix}\n",
    "\\end{align}\n",
    "with $i=1,2,3$ and note that\n",
    "\\begin{align}\n",
    "  \\bar\\psi = \\psi^\\dagger \\gamma^0\\,.\n",
    "\\end{align}\n",
    "The total Lagrangian of a noninteracting, massive spin $1/2$ particle\n",
    "of mass $m$ is thus given by\n",
    "\\begin{align}\\label{eqn:origlagrspinoh}\n",
    "  {\\cal L} &= \\bar\\psi (i \\gamma^\\mu \\partial_\\mu - m) \\psi\\,.\n",
    "\\end{align}\n",
    "It is apparent that this Lagrangian is also invariant under\n",
    "translations of space and time.  The corresponding equation of motion\n",
    "is the **Dirac equation** of a free spin $1/2$ field\n",
    "\\begin{align}\n",
    "  (i\\gamma^\\mu \\partial_\\mu - m)\\psi = (i D - m) \\psi = 0\n",
    "\\end{align}\n",
    "with **Dirac operator** $D = \\gamma^\\mu \\partial_\\mu$.\n",
    "\n",
    "Note that we do not have to consider $\\psi_-$ and $\\psi_+$ as\n",
    "independent fields.  If we identify\n",
    "\\begin{align}\n",
    "  \\psi_- = i \\sigma^2\\psi_+^*\n",
    "\\end{align}\n",
    "it follows from $\\sigma_2 \\sigma_i^* \\sigma_2 =-\\sigma_i$, see\n",
    "Eq. \\eqref{eqn:sigmatwo}, that under restricted Lorentz transformations\n",
    "\\begin{align}\n",
    "  \\psi'_- &= i\\sigma^2 [\\exp[(s_i-i\\phi_i)\\sigma_i/2]\\psi_+]^* =\n",
    "  \\exp[(s_i+i\\phi_i)\\sigma^2\\sigma_i^*\\sigma^2/2] i\\sigma^2\\psi_+^*\n",
    "  \\notag\\\\&= \\exp[(-s_i-i\\phi_i)\\sigma_i/2] \\psi_-\\,,\n",
    "\\end{align}\n",
    "in accordance with Eq. \\eqref{eqn:parityact}.  The mass terms then\n",
    "become Majorana mass terms, and it can be shown that the fields\n",
    "$\\psi_+$ become their own antiparticles.  This, however, implies that\n",
    "they are not allowed to carry a nonzero charge and therefore excludes\n",
    "this scenario for the quarks of QCD."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "6a3e8159",
   "metadata": {},
   "source": [
    "## Gamma matrices and Lorentz structure\n",
    "Before we continue with the discussion of gauge symmetries a few notes\n",
    "about the algebra of gamma matrices are in order.  The gamma matrices\n",
    "satisfy the **Clifford-algebra** relation\n",
    "\\begin{align}\n",
    "  \\{ \\gamma^\\mu, \\gamma^\\nu \\} = \\gamma^\\mu \\gamma^\\nu + \\gamma^\\nu\n",
    "  \\gamma^\\mu = 2 g^{\\mu\\nu}\\,.\n",
    "\\end{align}\n",
    "The parity operator can be written in terms of $\\gamma^0$ as\n",
    "\\begin{align}\n",
    "  D(P)\\psi = \\gamma^0 \\psi\\,.\n",
    "\\end{align}\n",
    "Furthermore, it is convenient to define\n",
    "\\begin{align}\n",
    "  \\gamma^5=i\\gamma^0 \\gamma^1 \\gamma^2 \\gamma^3 =\n",
    "  \\begin{pmatrix}\n",
    "    -\\mathbb{1} & \\\\ & \\mathbb{1}\n",
    "  \\end{pmatrix}\n",
    "\\end{align}\n",
    "which allows to project on the $-$ and $+$ sectors by\n",
    "\\begin{align}\n",
    "  P_\\pm = \\frac{\\mathbb{1} \\pm \\gamma^5}2\\,,\n",
    "\\end{align}\n",
    "where $\\mathbb{1}$ is the identity matrix in the respective space.  The\n",
    "matrix $\\gamma^5$ anticommutes with all other gamma matrices,\n",
    "\\begin{align}\n",
    "  \\{ \\gamma^5, \\gamma^\\mu \\} = 0\n",
    "\\end{align}\n",
    "with $\\mu=0,1,2,3$.  \n",
    "\n",
    "Note that the gamma matrices can be used to construct field bilinears\n",
    "that transform in a well-defined way under the orthochronous Lorentz\n",
    "group.  Under restricted Lorentz transformations $\\Lambda$ we find\n",
    "\\begin{align}\n",
    "  v^\\mu &= \\bar\\psi \\gamma^\\mu \\psi &\\to && {\\Lambda^\\mu}_\\nu v^\\nu\\,,\\notag\\\\\n",
    "  a^\\mu &= \\bar\\psi \\gamma^\\mu \\gamma^5 \\psi &\\to&& {\\Lambda^\\mu}_\\nu a^\\nu\\,,\\notag\\\\\n",
    "  s &= \\bar\\psi \\psi &\\to&& s\\,,\\notag\\\\\n",
    "  p &= \\bar\\psi \\gamma^5 \\psi &\\to&& p\\,,\n",
    "\\end{align}\n",
    "see Eq. \\eqref{eqn:lorentzvector}.  The action of parity $P$ on $v$,\n",
    "$a$, $s$ and $p$ is given by\n",
    "\\begin{align}\n",
    "  v^\\mu &&\\to&& \n",
    "  -v^\\mu + 2 g^{0\\mu} v^0\\,,\\notag\\\\\n",
    "  a^\\mu &&\\to&& a^\\mu -2 g^{0\\mu} a^0\\,,\\notag\\\\\n",
    "  s &&\\to&& s\\,,\\notag\\\\\n",
    "  p &&\\to&& -p\\,.\n",
    "\\end{align}\n",
    "Therefore $v^\\mu$ transforms as a vector, $a^\\mu$ transforms as an\n",
    "axial vector, $s$ transforms as a scalar and $p$ transforms as a\n",
    "pseudoscalar."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "6923e861",
   "metadata": {},
   "source": [
    "## Gauge symmetry\n",
    "In the last section we have constructed a relativistically invariant\n",
    "Lagrangian of a massive spin $1/2$ field.  Up to now the particles\n",
    "represented by the field do not interact with each other.  In the\n",
    "following we add a local internal symmetry (or **gauge symmetry**)\n",
    "to the Lagrangian and show that such a modification introduces an\n",
    "interaction between the spin $1/2$ particles that is mediated by\n",
    "massless spin $1$ particles."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d9747815",
   "metadata": {},
   "source": [
    "### Internal symmetries\n",
    "Consider the Lagrangian of Eq. \\eqref{eqn:origlagrspinoh}, i.e.,\n",
    "\\begin{align}\\tag{\\ref{eqn:origlagrspinoh}}\n",
    "  {\\cal L} &= \\bar\\psi (i \\gamma^\\mu \\partial_\\mu - m) \\psi\n",
    "\\end{align}\n",
    "with fields $\\psi$ in spinor space.  The operation of the matrices\n",
    "$\\gamma^\\mu$ on $\\psi$ is given by the matrix-vector multiplication in\n",
    "this space.  The most trivial way to add an additional symmetry\n",
    "$S_i$ is to choose a new symmetry group $S$ that is a direct product of\n",
    "the Poincaré symmetry group $S_p$ and $S_i$,\n",
    "\\begin{align}\n",
    "  S = S_p \\otimes S_i\\,.\n",
    "\\end{align}\n",
    "In such a modification we call $S_i$ an **internal symmetry** of the\n",
    "Lagrangian.  The fields $\\psi$ must transform in representations\n",
    "of the bigger symmetry group $S$ and therefore live in a product space\n",
    "of the spinor space and the vector space of the internal symmetry."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "a3509f13",
   "metadata": {},
   "source": [
    "### Local symmetries\n",
    "Let us choose $S_i$ to consist of spacetime-dependent transformations\n",
    "of $\\psi(x)$ with infinitesimal transformations $G(x)$ defined by\n",
    "\\begin{align}\\label{eqn:loctr}\n",
    "  \\psi(x) \\to \\psi(x) + i G(x) \\psi(x)\\,,\n",
    "\\end{align}\n",
    "where the action of $G(x)$ on $\\psi(x)$ is the matrix-vector\n",
    "multiplication in the internal symmetry space.  We ignore terms of\n",
    "order $G^2$ throughout the remainder of this section.  The mass\n",
    "term of Eq. \\eqref{eqn:origlagrspinoh},\n",
    "\\begin{align}\n",
    "  {\\cal L}_\\text{mass} &= m \\bar\\psi \\psi\\,,\n",
    "\\end{align}\n",
    "is symmetric under Eq. \\eqref{eqn:loctr} if\n",
    "\\begin{align}\n",
    "  G(x)^\\dagger = G(x)\\,,\n",
    "\\end{align}\n",
    "i.e., if $G(x)$ generates unitary transformations.  The kinetic term\n",
    "\\begin{align}\n",
    "  {\\cal L}_\\text{kinetic} &= \\bar\\psi (i \\gamma^\\mu \\partial_\\mu)\n",
    "  \\psi\\,,\n",
    "\\end{align}\n",
    "however, transforms to\n",
    "\\begin{align}\n",
    "  {\\cal L}'_\\text{kinetic} &= {\\cal L}_\\text{kinetic} -\\bar\\psi (\\partial_\\mu G(x)) \\gamma^\\mu \\psi\n",
    "\\end{align}\n",
    "under unitary transformations.  An invariant term can only be\n",
    "constructed if we replace\n",
    "\\begin{align}\n",
    "  \\partial_\\mu \\to D_\\mu\n",
    "\\end{align}\n",
    "with\n",
    "\\begin{align}\\label{eqn:defcovder}\n",
    "  D_\\mu \\to [\\mathbb{1}+i G(x)] D_\\mu [\\mathbb{1}-i G(x)^\\dagger]\n",
    "\\end{align}\n",
    "under Eq. \\eqref{eqn:loctr}, where $\\mathbb{1}$ is the identity matrix.  We\n",
    "call $D_\\mu$ a **covariant derivative**.  The covariant derivative\n",
    "has to generate a kinetic term for the spin $1/2$ fields, and therefore\n",
    "we use the ansatz\n",
    "\\begin{align}\n",
    "  D_\\mu = \\partial_\\mu + i A_\\mu\\,,\n",
    "\\end{align}\n",
    "where $A_\\mu$ has to transform under $S_i$ in a way that satisfies\n",
    "Eq. \\eqref{eqn:defcovder}.  Note that $A_\\mu$ can act non-trivially on\n",
    "the internal symmetry space.  Since $\\partial_\\mu$ is an\n",
    "anti-Hermitian operator, we require $A=A^\\dagger$ so that the\n",
    "Lagrangian is real.  In accordance with Eq. \\eqref{eqn:defcovder} we\n",
    "request that the transformed $A'_\\mu$ satisfies\n",
    "\\begin{align}\n",
    "  \\partial_\\mu + i A'_\\mu &= [\\mathbb{1}+i G(x)] (\\partial_\\mu + i A_\\mu) [\\mathbb{1}-i G(x)^\\dagger] \\notag\\\\\n",
    "  &= \\partial_\\mu + i A_\\mu -i (\\partial_\\mu G(x)) - [G(x), A_\\mu]\\,,\n",
    "\\end{align}\n",
    "where we used the Hermiticity of $G(x)$.  Thus we can construct an\n",
    "invariant kinetic term if $A_\\mu$ transforms as\n",
    "\\begin{align}\n",
    "  A'_\\mu = A_\\mu- (\\partial_\\mu G(x)) + i [G(x), A_\\mu]\\,.\n",
    "\\end{align}\n",
    "We conclude that we can construct a Lagrangian\n",
    "\\begin{align}\n",
    "  {\\cal L} &= \\bar\\psi (i \\gamma^\\mu D_\\mu - m) \\psi = \\bar\\psi (i\n",
    "  \\gamma^\\mu \\partial_\\mu - m) \\psi - \\bar\\psi \\gamma^\\mu A_\\mu \\psi\n",
    "\\end{align}\n",
    "that is invariant under the symmetry group $S$ with internal symmetry\n",
    "$S_i$ defined by the infinitesimal transformation\n",
    "\\begin{align}\\label{eqn:deflocsymnew}\n",
    "  \\psi(x) &\\to \\psi(x) + i G(x) \\psi(x)\\,,\\notag\\\\\n",
    "  A_\\mu & \\to A_\\mu- (\\partial_\\mu G(x)) + i [G(x), A_\\mu]\\,.\n",
    "\\end{align}\n",
    "Note that the $A_\\mu$ also transform in the fundamental representation\n",
    "of the restricted Lorentz symmetry group,\n",
    "\\begin{align}\n",
    "  A_\\mu \\to {\\Lambda_\\mu}^\\nu A_\\nu\n",
    "\\end{align}\n",
    "under Lorentz transformation $\\Lambda$.  The spin operator of the\n",
    "fundamental representation of the restricted Lorentz group is given by\n",
    "$S_j = i J_j$ with $S^2 = s(s+1)$ and $s=1$.  Therefore we have\n",
    "introduced fields $A_\\mu$ of spin $1$ that interact with the spin\n",
    "$1/2$ fields due to the term\n",
    "\\begin{align}\n",
    " {\\cal L}_\\text{interaction} &=  -\\bar\\psi \\gamma^\\mu A_\\mu \\psi\n",
    "\\end{align}\n",
    "in the Lagrangian.  Since the Lagrangian has to be invariant under\n",
    "Eq. \\eqref{eqn:deflocsymnew} the fields $A_\\mu$ are not allowed to\n",
    "have a quadratic mass term and must therefore correspond to massless\n",
    "particles.  They can, however, have a kinetic term that allows them to\n",
    "propagate in spacetime.  To second order in $\\partial_\\mu A_\\nu$ the\n",
    "only term that is invariant under $S_i$ and $S_p$ is proportional to\n",
    "\\begin{align}\n",
    " {\\cal L}_\\text{YM} \\propto  \\textrm{Tr}\\: [D_\\mu,D_\\nu][D^\\mu,D^\\nu]\\,,\n",
    "\\end{align}\n",
    "where the trace $\\textrm{Tr}$ acts on the internal symmetry space.  This is\n",
    "the **Yang-Mills** term.  The invariance under $S_i$ is\n",
    "due to the covariance of\n",
    "\\begin{align}\n",
    "  [D_\\mu,D_\\nu] \\to [\\mathbb{1}+i G(x)] [D_\\mu,D_\\nu] [\\mathbb{1}-i G(x)^\\dagger]\n",
    "\\end{align}\n",
    "under Eq. \\eqref{eqn:deflocsymnew}.  We define the **field-strength tensor**\n",
    "\\begin{align}\n",
    "  F_{\\mu\\nu} &= -i[D_\\mu,D_\\nu] =\n",
    "  (\\partial_\\mu A_\\nu) - (\\partial_\\nu A_\\mu) + i[A_\\mu,A_\\nu]\n",
    "\\end{align}\n",
    "and express the total Lagrangian conveniently as\n",
    "\\begin{align}\n",
    "  {\\cal L} &= \\bar\\psi (i \\gamma^\\mu D_\\mu - m) \\psi + \\frac12\\textrm{Tr}\\:\n",
    "  F_{\\mu\\nu}F^{\\mu\\nu}\\,.\n",
    "\\end{align}\n",
    "\n",
    "Note that if the local symmetry group is not abelian, the term\n",
    "$[A_\\mu,A_\\nu]$ introduces a self-interaction between the massless\n",
    "spin $1$ particles.\n",
    "\n",
    "If we choose the first non-trivial unitary symmetry group $\\textrm{U}(1)$, we\n",
    "recover the theory of electrodynamics coupled to a spin $1/2$ field.\n",
    "The photons are now given by the spin $1$ fields $A_\\mu$.  The\n",
    "equations of motion of the fields $A_\\mu$ can readily be identified\n",
    "with Maxwell's equations of electrodynamics."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "96d72447",
   "metadata": {},
   "source": [
    "## The Lagrangian of QCD\n",
    "The internal symmetry group of QCD is given by $\\textrm{SU}(3)$.  If we choose\n",
    "$A_\\mu$ to live in the group algebra of $\\textrm{SU}(3)$, we can write\n",
    "\\begin{align}\n",
    "  A_\\mu = A_\\mu^a \\lambda_a\\,,\n",
    "\\end{align}\n",
    "where the matrices $\\lambda_a$ ($a=1,\\ldots,8$) span the algebra of\n",
    "$\\textrm{SU}(3)$.  The eight fields $A_\\mu^a$ now correspond to eight\n",
    "independent gluons.  The quarks live in the internal symmetry space of\n",
    "$\\textrm{SU}(3)$.  Its fundamental representation is three-dimensional and\n",
    "therefore there are three different quark fields, or three different\n",
    "**colors** of quarks (The name quantum chromodynamics is due\n",
    "  to this interpretation of the three different quark fields as\n",
    "  different colors of quarks.).  The bound states of quarks and\n",
    "anti-quarks, called **hadrons**, must transform as singlets of $S_i$\n",
    "and are therefore **color neutral**.  There are two types of hadrons:\n",
    "**mesons** and **baryons**.  Mesons, such as the pion, are bosonic\n",
    "hadrons that consist of a quark and an anti-quark.  Baryons, such as\n",
    "the proton or neutron, are fermionic hadrons that consist of three\n",
    "quarks.\n",
    "\n",
    "Note that since $\\textrm{SU}(3)$ is non-abelian, gluons are self-interacting.\n",
    "This property can be shown to lead to the **asymptotic freedom** of QCD, i.e., for high energies the strength of the\n",
    "interaction becomes weaker.\n",
    "\n",
    "We rescale the fields $A_\\mu \\to g A_\\mu$ and change the prefactor of\n",
    "the kinetic term of gluons so that we can adjust the strength of the\n",
    "interaction of quarks and gluons explicitly.  The Lagrangian of a\n",
    "quark coupled to the gluons then reads\n",
    "\\begin{align}\n",
    "  {\\cal L} &= \\bar\\psi (i \\gamma^\\mu \\partial_\\mu - m) \\psi + \\frac12\\textrm{Tr}\\:\n",
    "  F_{\\mu\\nu}F^{\\mu\\nu} - g A_\\mu^a \\bar\\psi \\gamma^\\mu \\lambda_a \\psi\n",
    "\\end{align}\n",
    "with\n",
    "\\begin{align}\n",
    "  F_{\\mu\\nu} &= (\\partial_\\mu A_\\nu) - (\\partial_\\nu A_\\mu) +i g\n",
    "  [A_\\mu,A_\\nu]\\,.\n",
    "\\end{align}\n",
    "It turns out that in nature there are more than one kind of quarks\n",
    "which differ by their mass and electromagnetic charge.  One currently\n",
    "has experimental evidence for $6$ different types of quarks, called\n",
    "different quark **flavors**, of which three have a fractional\n",
    "electromagnetic charge of $+2/3$ and three have a fractional\n",
    "electromagnetic charge of $-1/3$.  Two quarks are very light and thus\n",
    "play an important role in the low-energy physics of QCD discussed in\n",
    "the remainder of this thesis.  They are called **up** and **down**\n",
    "quarks (corresponding to their respective fractional electromagnetic\n",
    "charges $+2/3$ and $-1/3$).  The next heavier quark is called **strange**\n",
    "quark and has a fractional electromagnetic charge of\n",
    "$-1/3$.  Their masses are related approximately by\n",
    "\\begin{align}\n",
    "  \\frac{m_s}{m_d} \\approx 20\\,,\\qquad \\frac{m_u}{m_d} \\approx \\frac12\\,,\n",
    "\\end{align}\n",
    "where $m_u$, $m_d$, $m_s$ are the masses of up, down and strange quark.\n",
    "Note that these relations are only\n",
    "order-of-magnitude estimates.  The total Lagrangian of QCD thus reads\n",
    "\\begin{align}\n",
    "  {\\cal L}_\\text{QCD} &= \\sum_{f=1}^6\\bar\\psi_f (i\n",
    "  \\gamma^\\mu \\partial_\\mu - m_f) \\psi_f + \\frac12\\textrm{Tr}\\: F_{\\mu\\nu}F^{\\mu\\nu}\n",
    "  \\notag\\\\&\\quad - g \\sum_{f=1}^6 A_\\mu^a \\bar\\psi_f \\gamma^\\mu\n",
    "  \\lambda_a \\psi_f\\,.\n",
    "\\end{align}\n",
    "We observe that, depending on the quark masses, the total Lagrangian\n",
    "has an additional symmetry in flavor space.  This symmetry will force some of the hadrons to be particularly light.  More on this in a later chapter."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
