## Gauge-covariant extensions of Killing tensors and conservation laws

J.W. van Holten a

## Nikhef

Science Park 105, Amsterdam NL

and

Lorentz Institute, Leiden University

Niels Bohrweg 2, Leiden NL

## Abstract

In classical and quantum mechanical systems on manifolds with gauge-field fluxes, constants of motion are constructed from gauge-covariant extensions of Killing vectors and tensors. This construction can be carried out using a manifestly covariant procedure, in terms of covariant phase space with a covariant generalization of the Poisson brackets, c.q. quantum commutators. Some examples of this construction are presented.

a e-mail: v.holten@nikhef.nl

## 1 Noether's theorem

This paper discusses symmetries and conservation laws in the context of hamiltonian dynamics. The discussion is framed predominantly in the language of classical dynamics, but the use of Poisson brackets and their correspondence with quantum commutators, guarantees that many results also apply to the operator formulation of quantum dynamics. The main difference is the operator ordering to be implemented in quantum theory, the technicalities of which are not relevant to the issues I focus on.

The connection between continuous symmetries and conservation laws is established by Noether's theorem [1]. I briefly review the theorem by considering infinitesimal transformations on phase-space variables ( x, p ) obtained from a generating function G x, p ( ) through the Poisson brackets

$$\delta x = \{ x, G \} = \frac { \partial G } { \partial p }, \quad \delta G = \{ p, G \} = - \frac { \partial G } { \partial x }.$$

Observe, that these variations are defined such that G itself is an invariant:

$$\delta G = \delta x \frac { \partial G } { \partial x } + \delta p \frac { \partial G } { \partial p } = \{ G, G \} = 0.$$

Under such transformations the hamiltonian of the system changes by

$$\delta H = \{ H, G \} = - \frac { d G } { d t },$$

the change of G along the phase-space trajectory ( x t , p ( ) ( )) t generated by the hamiltonian H . It follows immediately, that G is a constant of motion if the hamiltonian is invariant under the transformations (1).

It is also of some interest to consider the variation of the action

$$S = \int _ { 1 } ^ { 2 } d t \left ( p \, \frac { d x } { d t } - H ( x, p ) \right ).$$

Applying the variations (1)

$$\delta S = \int _ { 1 } ^ { 2 } d t \left [ \frac { d } { d t } \left ( p \, \frac { \partial G } { \partial p } - G \right ) - \{ H, G \} \right ] = \left [ p \frac { \partial G } { \partial p } - G \right ] _ { 1 } ^ { 2 }.$$

Thus variations under which the hamiltonian is invariant, leave the action invariant modulo boundary terms. This is sufficient for G to be a constant of motion.

## 2 Isometries of manifolds

On a manifold with (local) co-ordinates x µ and metric g µν ( x ) the geodesics can be obtained as the trajectories of test-particles with proper-time hamiltonian

$$H = \frac { 1 } { 2 } \, g ^ { \mu \nu } p _ { \mu } p _ { \nu }.$$

Indeed, using the overdot notation for proper-time derivatives, the hamilton equations take the form

$$\dot { x } ^ { \mu } = g ^ { \mu \nu } p _ { \nu }, \quad \dot { p } _ { \mu } = - \frac { 1 } { 2 } \, \frac { \partial g ^ { \nu \lambda } } { \partial x ^ { \mu } } \, p _ { \nu } p _ { \lambda } = \frac { 1 } { 2 } \, \frac { \partial g _ { \nu \lambda } } { \partial x ^ { \mu } } \, \dot { x } ^ { \nu } \dot { x } ^ { \lambda }.$$

The last expression is equivalent to the geodesic equation

$$\ddot { x } ^ { \mu } + \Gamma _ { \lambda \nu } ^ { \ \mu } \, \dot { x } ^ { \lambda } \dot { x } ^ { \nu } = 0.$$

In this language, isometries of the manifold are found as constants of motion which are linear in the momentum:

$$J ( x, p ) = J ^ { \mu } ( x ) p _ { \mu }, \quad \left \{ J, H \right \} = ( \nabla _ { \mu } J _ { \nu } ) \, p ^ { \mu } p ^ { \nu } = 0,$$

where (in a somewhat hybrid notation) the contravariant components of the momentum are p µ = ˙ x µ . Hence the covariant coefficient functions J µ form a Killing vector, a solution of the Killing equation

$$\nabla _ { \mu } J _ { \nu } + \nabla _ { \nu } J _ { \mu } = 0 \quad \Leftrightarrow \quad J ^ { \lambda } \frac { \partial g _ { \mu \nu } } { \partial x ^ { \lambda } } + \frac { \partial J ^ { \lambda } } { \partial x ^ { \mu } } \, g _ { \lambda \nu } + \frac { \partial J ^ { \lambda } } { \partial x ^ { \nu } } \, g _ { \mu \lambda } = 0. \quad \ \ ( 1 0 )$$

The second (contravariant) form of the equation states that the Lie-derivative of the metric w.r.t. the vector J µ vanishes, which is the usual definition of an isometry. Also note, that the constants of motion defined by Killing vectors are precisely those, for which

$$p \, \frac { \partial G } { \partial p } = G,$$

and therefore generates transformations under which the action is strictly invariant. This is to be expected for an isometry which by construction leaves the line element ds 2 = g µν dx dx µ ν invariant.

Although the coordinate transformations δx µ generated by Killing vectors thus have an elegant interpretation, this does not hold for the corresponding transformations of the canonical momenta p µ :

$$\delta p _ { \mu } = A _ { \mu } ^ { \nu } ( x ) p _ { \nu }, \quad A _ { \mu } ^ { \nu } = - \frac { \partial J ^ { \nu } } { \partial x ^ { \mu } }.$$

This transformation rule is not general covariant; as A ν µ ( x ) is point-dependent, covariance requires δp to be corrected for the parallel displacement generated by

the translation δx µ . Thus a set of covariant transformations in phase space is defined by

$$\Delta x ^ { \mu } = \delta x ^ { \mu }, \quad \Delta p _ { \mu } = \delta p _ { \mu } - \delta x ^ { \lambda } \Gamma _ { \lambda \mu } ^ { \ \nu } p _ { \nu }.$$

In order for these transformations to respect the Poisson brackets (2) and (3), it is then necessary to introduce a covariant derivative

$$\mathcal { D } _ { \mu } G = \frac { \partial G } { \partial x ^ { \mu } } + \Gamma _ { \mu \nu } ^ { \ \lambda } p _ { \lambda } \, \frac { \partial G } { \partial p _ { \nu } } = - \Delta p _ { \mu },$$

such that

$$\Delta G = \Delta x ^ { \mu } \mathcal { D } _ { \mu } G + \Delta p _ { \mu } \frac { \partial G } { \partial p _ { \mu } } = \{ G, G \} = 0,$$

and

$$\Delta H = \Delta x ^ { \mu } \mathcal { D } _ { \mu } H + \Delta p _ { \mu } \frac { \partial H } { \partial p _ { \mu } } = \{ H, G \} = - \frac { d G } { d \tau }.$$

Thus we have constructed a covariant expression for the Poisson bracket of two arbitrary scalar phase-space functions:

$$\{ G, K \} = \mathcal { D } _ { \mu } G \, \frac { \partial K } { \partial p _ { \mu } } - \frac { \partial G } { \partial p _ { \mu } } \, \mathcal { D } _ { \mu } K.$$

## 3 Killing tensors

The metric postulate guarantees that the geodesic hamiltonian (6) is covariantly constant:

$$\nabla _ { \lambda } g _ { \mu \nu } = 0 \quad \Leftrightarrow \quad \mathcal { D } H = 0.$$

The condition for a constant of geodesic motion then takes the simple form

$$\{ G, H \} = p ^ { \mu } \mathcal { D } _ { \mu } G = 0.$$

Applying this to a general expression of homogeneous rank n in the momenta

$$G ( x, p ) = G ^ { \mu _ { 1 } \dots \mu _ { n } } ( x ) p _ { \mu _ { 1 } } \dots p _ { \mu _ { n } } \quad \Rightarrow \quad p ^ { \mu } \mathcal { D } _ { \mu } G = \left ( \nabla _ { \mu _ { n + 1 } } G _ { \mu _ { 1 } \dots \mu _ { n } } \right ) p ^ { \mu _ { 1 } } \dots p ^ { \mu _ { n + 1 } }, \ \ ( 2 0 )$$

the condition for a constant of motion becomes a generalization of the Killing equation (10):

$$\nabla _ { ( \mu _ { n + 1 } } G _ { \mu _ { 1 } \dots \mu _ { n } ) } = 0.$$

The solutions of these equations are therefore known as Killing tensors [3, 4].

The geometrical interpretation of the transformations generated by constants of motion constructed from Killing tensors of rank 2 or higher, is more complicated than for Killing vectors. They do not generate transformations in the manifold, but in the tangent bundle, the physical phase-space:

$$\Delta x ^ { \mu } = \frac { \partial G } { \partial p _ { \mu } } = n G ^ { \mu \mu _ { 1 } \dots \mu _ { n - 1 } } p _ { \mu _ { 1 } } \dots p _ { \mu _ { n - 1 } },$$

$$\Delta p _ { \mu } = - \mathcal { D } _ { \mu } G = - \left ( \nabla _ { \mu } G ^ { \mu _ { 1 } \dots \mu _ { n } } \right ) p _ { \mu _ { 1 } } \dots p _ { \mu _ { n } }.$$

Note, that for transformations with n ≥ 2 the action actually transforms by boundary terms:

$$p _ { \mu } \, \frac { \partial G } { \partial p _ { \mu } } = n G \quad \Rightarrow \quad \Delta S = ( n - 1 ) \left [ G ( x _ { 2 }, p _ { 2 } ) - G ( x _ { 1 }, p _ { 1 } ) \right ]. \quad \quad ( 2 3 )$$

Obviously, as G is a constant of motion the action is still invariant along geodesics.

## Example: Kerr geometry [2, 3]

The Kerr geometry is a Ricci-flat geometry in 4-dimensional space-time: R µν = 0, implying it is a solution of the Einstein equations of General Relativity in empty space-time. This solutions describes a rotating black hole with constant mass M and angular momentum J . The geodesic hamiltonian reads

$$H = \frac { 1 } { 2 \rho ^ { 2 } } \left [ \Delta ^ { 2 } p _ { r } ^ { 2 } + p _ { \theta } ^ { 2 } + \left ( a \sin \theta \, p _ { t } + \frac { p _ { \varphi } } { \sin \theta } \right ) ^ { 2 } - \frac { 1 } { \Delta ^ { 2 } } \left ( ( r ^ { 2 } + a ^ { 2 } ) p _ { t } + a p _ { \varphi } \right ) ^ { 2 } \right ], \ \ ( 2 4 )$$

where a = J/M is the angular momentum per unit of mass, and the notation follows standard conventions:

$$\Delta ^ { 2 } = r ^ { 2 } - 2 M r + a ^ { 2 }, \quad \rho ^ { 2 } = r ^ { 2 } + a ^ { 2 } \cos ^ { 2 } \theta.$$

This space-time geometry admits a rank-2 Killing tensor, coding the well-known Carter constant of motion:

$$\text{canvas convvex} & \text{in move}. \\ K \ = \ \frac { 1 } { 2 \rho ^ { 2 } } \left [ - \Delta ^ { 2 } a ^ { 2 } \cos ^ { 2 } \theta \, p _ { r } ^ { 2 } + r ^ { 2 } p _ { \theta } ^ { 2 } + r ^ { 2 } \left ( a p _ { t } + \frac { p _ { \varphi } } { \sin ^ { 2 } \theta } \right ) ^ { 2 } \\ & \quad + \frac { a ^ { 2 } \cos ^ { 2 } \theta } { \Delta ^ { 2 } } \left ( ( r ^ { 2 } + a ^ { 2 } ) p _ { t } + a p _ { \varphi } \right ) ^ { 2 } \right ]. \\ - - \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot$$

This constant of motion can be generalized to neutral and charged, spinless and spinning particles in Kerr-Newman space-time of charged black holes [3, 5, 6, 7].

## 4 Bracket algebra

The generators of hamiltonian symmetry transformations, defining constants of motion, define a Lie-algebra by their bracket relations (17); this follows from the Jacobi identity

$$\left \{ \left \{ G, K \right \}, J \right \} + \left \{ \left \{ K, J \right \}, G \right \} + \left \{ \left \{ J, G \right \}, K \right \} = 0.$$

Applying the identity to the special case where one of the functions is the hamiltonian: J = H , implies that constants of motion ( G,K ) satisfy

$$\left \{ \left \{ G, K \right \}, H \right \} = 0.$$

Hence the bracket of two constants of motion produces another constant of motion, and the set of such constants is closed under the bracket operation. Here we summarize some properties of this algebra [4].

The generators J = J µ p µ linear in momentum define a Lie subalgebra of the full algebra:

$$\{ J _ { 1 }, J _ { 2 } \} = J _ { 3 }, \quad J _ { 3 } ^ { \mu } = J _ { 2 } ^ { \nu } \, \nabla _ { \nu } J _ { 1 } ^ { \mu } - J _ { 1 } ^ { \nu } \, \nabla _ { \nu } J _ { 2 } ^ { \mu }.$$

The constants of motion of rank n ≥ 2 then define representations of this subalgebra , characterized by the transformation rule 1

$$\{ G _ { 1 }, J \} = G _ { 2 }, \quad G _ { 2 } ^ { \mu _ { 1 } \dots \mu _ { n } } = J ^ { \lambda } \nabla _ { \lambda } G ^ { \mu _ { 1 } \dots \mu _ { n } } - n G ^ { \lambda ( \mu _ { 1 } \dots \mu _ { n - 1 } } \nabla _ { \lambda } J ^ { \mu _ { n } ) }. \quad ( 2 9 )$$

The Jacobi identity

$$\left \{ \left \{ G, J _ { 1 } \right \}, J _ { 2 } \right \} - \left \{ \left \{ G, J _ { 2 } \right \}, J _ { 1 } \right \} = \left \{ G, \left \{ J _ { 1 }, J _ { 2 } \right \} \right \},$$

then guarantees that the transformations generated by the linear J x, p ( ) on the rankn tensors G x, p ( ) obey the composition rules of the subalgebra spanned by the linear J x, p ( ) themselves.

As for the brackets of a generator G ( n ) of rank n , and a generator G ( m ) of rank m , with both ( n, m ) ≥ 2, it is easily recognized that any non-vanishng brackets of such quantities must be a generator G ( n + m -1) , of rank n + m -1 exceeding both n and m :

$$\left \{ G ^ { ( n ) }, G ^ { ( m ) } \right \} \sim G ^ { ( n + m - 1 ) }.$$

Therefore non-vanishing brackets of higher-rank generators will potentially lead to an infinite set of generators of arbitrary high rank. Well-known examples of such algebras are the Virasoro and Kac-Moody algebras of 2-D conformal field theories.

## 5 Abelian gauge interactions

Geodesic motion on a manifold applies to the motion of a pure mass point subject only to geometrical forces. However, Noether's theorem is very general and can be applied also in presence of external force fields [8, 9]. In this section I consider abelian gauge interactions transmitted by a vector field A µ ( x ), acting on a point mass m with charge q . For such a particle the canonical hamiltonian reads

$$H = \frac { 1 } { 2 m } \, g ^ { \mu \nu } \left ( p _ { \mu } - q A _ { \mu } \right ) \left ( p _ { \nu } - q A _ { \nu } \right ).$$

The hamiltonian equations of motion (7), (8) are generalized to

$$p _ { \mu } = m g _ { \mu \nu } \dot { x } ^ { \nu } + q A _ { \mu }, \quad g _ { \mu \nu } \left ( \ddot { x } ^ { \nu } + \Gamma _ { \kappa \lambda } ^ { \ \nu } \dot { x } ^ { \kappa } \dot { x } ^ { \lambda } \right ) = \frac { q } { m } \, F _ { \mu \nu } \, \dot { x } ^ { \nu },$$

1 The parenthesis around the indices ( µ ...µ 1 n ) denote complete symmetrization with unit weight.

which is the Lorentz force law on curved manifolds.

A drawback of this hamiltonian formulation is, that the canonical momentum is not gauge invariant; under gauge transformations

$$A ^ { \prime } _ { \mu } = A _ { \mu } + \nabla _ { \mu } \Lambda, \quad p ^ { \prime } _ { \mu } = p _ { \mu } + q \nabla _ { \mu } \Lambda.$$

Therefore it is preferable to work with the covariant momentum [7, 10]

$$\pi _ { \mu } = p _ { \mu } - q A _ { \mu } = m g _ { \mu \nu } \dot { x } ^ { \nu },$$

in terms of which the hamiltonian takes the simple form

$$H = \frac { 1 } { 2 m } \, g ^ { \mu \nu } \pi _ { \mu } \pi _ { \nu }.$$

If the covariant derivative D now is generalized to

$$\mathcal { D } _ { \mu } G = \frac { \partial G } { \partial x ^ { \mu } } + \Gamma _ { \mu \nu } ^ { \ \lambda } \pi _ { \lambda } \, \frac { \partial G } { \partial \pi _ { \nu } },$$

then as before the metric postulate can be used to show that

$$\mathcal { D } _ { \mu } H = 0.$$

Moreover, the correct dynamics is reproduced as usual by

$$\frac { d G } { d \tau } = \left \{ G, H \right \},$$

when supplemented by the covariant brackets

$$\{ G, K \} = \mathcal { D } _ { \mu } G \, \frac { \partial K } { \partial \pi _ { \mu } } - \frac { \partial G } { \partial \pi _ { \mu } } \, \mathcal { D } _ { \mu } K + q F _ { \mu \nu } \, \frac { \partial G } { \partial \pi _ { \mu } } \frac { \partial K } { \partial \pi _ { \nu } }.$$

In particular, this bracket reproduces the Ricci identity in the form

$$\{ \pi _ { \mu }, \pi _ { \nu } \} = q F _ { \mu \nu }.$$

In view of the identity (38) and the anti-symmetry of the field strength tensor F µν , the condition for a scalar quantity G x, π ( ) to be constant of motion becomes [10]

$$\{ G, H \} = 0 \quad \Rightarrow \quad \pi ^ { \mu } \mathcal { D } _ { \mu } G = q \pi ^ { \mu } F _ { \mu \nu } \, \frac { \partial G } { \partial \pi _ { \nu } }.$$

In particular, if G can be expressed as a power series

$$G ( x, \pi ) = \sum _ { n } \frac { 1 } { n! } \, G ^ { ( n ) \, \mu _ { 1 } \dots \mu _ { n } } ( x ) \, \pi _ { \mu _ { 1 } } \dots \pi _ { \mu _ { n } },$$

then eq. (42) takes the form of a first-order p.d.e.

$$\nabla _ { ( \mu _ { 1 } } G ^ { ( n ) } _ { \mu _ { 2 } \dots \mu _ { n + 1 } ) } = q F ^ { \ \nu } _ { ( \mu _ { 1 } } G ^ { ( n + 1 ) } _ { \mu _ { 2 } \dots \mu _ { n + 1 } ) \nu }.$$

This is a generalization of equation (21) for Killing tensors, forming a hierarchy of equation connecting tensors of different rank. Nevertheless, the existence of a Killing tensor of rank n still provides a constant of motion: it allows the series expansion (43) to be truncated, with all higher-rank components vanishing:

$$\nabla _ { ( \mu _ { 1 } } G ^ { ( n ) } _ { \mu _ { 2 } \dots \mu _ { n + 1 } ) } = 0 \ \Rightarrow \ \ G ^ { ( n + k ) } _ { \mu _ { 1 } \dots \mu _ { n + k } } = 0, \ \forall k \geq 0,$$

while all lower-rank components are obtained by solving eq. (44) with the allready known higher-rank ones as inhomegenous source term on the right-hand side [10].

Finally note, that it is straightforward to include a scalar potential Φ( x ) as well:

$$H = \frac { 1 } { 2 m } \, g ^ { \mu \nu } \pi _ { \mu } \pi _ { \nu } + \Phi, \ \Rightarrow \ \mathcal { D } _ { \mu } H = \frac { \partial \Phi } { \partial x ^ { \mu } }.$$

This modifies the generalization of the Killing equation to include a gradient term [11]

$$\pi ^ { \mu } \mathcal { D } _ { \mu } G = q \pi ^ { \mu } F _ { \mu \nu } \frac { \partial G } { \partial \pi _ { \nu } } + \frac { \partial \Phi } { \partial x ^ { \mu } } \frac { \partial G } { \partial \pi _ { \mu } }.$$

## Example: a quantum-dot model

A model of a quantum-dot, consisting of 2 electrons moving as a non-relativistic bound pair with Coulomb interaction in a confining harmonic potential and a magnetic field, was studied in ref. [12]. Ignoring the center-of-mass motion, the pair can be described in axial co-ordinates ( ρ, z, ϕ ) as a single particle with hamiltonian

$$H _ { C M } = \frac { 1 } { 2 } \, g ^ { i j } \pi _ { i } \pi _ { j } + \Phi,$$

where

$$\text{where} \quad g _ { i j } = \text{diag} ( 1, 1, \rho ^ { 2 } ), \quad \Phi = \frac { 1 } { 2 } \left ( \omega _ { 0 } ^ { 2 } \rho ^ { 2 } + \omega _ { z } ^ { 2 } z ^ { 2 } \right ) - \frac { \kappa } { \sqrt { \rho ^ { 2 } + z ^ { 2 } } }. \quad \quad ( 4 9 )$$

The definition of the dynamics is completed by the brackets

$$\{ G, K \} = \mathcal { D } _ { i } G \frac { \partial K } { \partial \pi _ { i } } - \frac { \partial G } { \partial \pi _ { i } } \mathcal { D } _ { i } K - 2 \omega _ { L } \rho \left ( \frac { \partial G } { \partial \pi _ { \rho } } \frac { \partial K } { \partial \pi _ { \varphi } } - \frac { \partial G } { \partial \pi _ { \varphi } } \frac { \partial K } { \partial \pi _ { \rho } } \right ),$$

where ω L = eB/ 2 is the Larmor frequency.

To guarantee that a Killing tensor can be completed to a full constant of motion, solving eq. (47), it may be necessary to tune the parameters in the scalar potential. In the case of the quantum dot, there is a rank-4 Killing tensor which can be completed to a full constant of motion provided the magnetic field is tuned

to a value such that the Larmor frequency is a fixed combination of the harmonic frequencies:

$$\omega _ { L } ^ { 2 } + \omega _ { 0 } ^ { 2 } = 4 \omega _ { z } ^ { 2 }.$$

Then the following expression is a solution of eq. (47) for this system:

$$G = \rho ^ { 2 } \pi _ { z } ^ { 4 } - 2 \rho z \pi _ { \rho } \pi _ { z } ^ { 3 } + z ^ { 2 } \pi _ { \rho } ^ { 2 } \pi _ { z } ^ { 2 } + \frac { 1 } { \rho ^ { 2 } } \pi _ { \varphi } ^ { 4 } + \pi _ { \rho } ^ { 2 } \pi _ { \varphi } ^ { 2 } + \left ( 2 + \frac { z ^ { 2 } } { \rho ^ { 2 } } \right ) \pi _ { z } ^ { 2 } \pi _ { \varphi } ^ { 2 }$$

$$G & = \rho ^ { 2 } \pi _ { z } ^ { 4 } - 2 \rho z \pi _ { \rho } ^ { 3 } z ^ { 2 } \pi _ { \rho } ^ { 2 } z ^ { 2 } + \frac { 1 } { \rho ^ { 2 } } \pi _ { \rho } ^ { 4 } + \pi _ { \rho } ^ { 2 } \pi _ { \rho } ^ { 2 } + \left ( 2 + \frac { z ^ { z } } { \rho ^ { 2 } } \right ) \pi _ { z } ^ { 2 } \pi _ { \rho } ^ { 2 } \\ & + 2 \omega _ { L } \pi _ { \varphi } \left ( \rho ^ { 2 } \pi _ { \rho } ^ { 2 } + ( 2 \rho ^ { 2 } + z ^ { 2 } ) \pi _ { z } ^ { 2 } \right ) + \left [ 2 \omega _ { z } ^ { 2 } z ^ { 3 } \rho + \frac { 2 \kappa z \rho } { \sqrt { \rho ^ { 2 } } + z ^ { 2 } } \right ] \pi _ { \rho } \pi _ { z } \\ & + \left [ ( 2 \omega _ { z } ^ { 2 } - \omega _ { 0 } ^ { 2 } ) z ^ { 2 } \rho ^ { 2 } + 2 \omega _ { L } ^ { 2 } \rho ^ { 4 } - \frac { 2 \kappa \rho ^ { 2 } } { \sqrt { \rho ^ { 2 } } + z ^ { 2 } } \right ] \pi _ { z } ^ { 2 } + \omega _ { L } ^ { 2 } \rho ^ { 4 } \pi _ { \rho } ^ { 2 } \\ & + \left [ 2 \omega _ { z } ^ { 2 } z ^ { 2 } + ( \omega _ { 0 } ^ { 2 } - 5 \omega _ { L } ^ { 2 } ) \rho ^ { 2 } - \frac { 2 \kappa } { \sqrt { \rho ^ { 2 } } + z ^ { 2 } } \right ] \pi _ { \varphi } ^ { 2 } \\ & - 2 \omega _ { L } \pi _ { \varphi } \left [ ( 3 \omega _ { L } ^ { 2 } - \omega _ { 0 } ^ { 2 } ) \rho ^ { 4 } - 2 \omega _ { z } ^ { 2 } \rho ^ { 2 } z ^ { 2 } + \frac { 2 \kappa \rho ^ { 2 } } { \sqrt { \rho ^ { 2 } } + z ^ { 2 } } \right ] + \omega _ { z } ^ { 4 } \rho ^ { 2 } z ^ { 4 } + 2 \omega _ { z } ^ { 2 } \omega _ { L } ^ { 2 } \rho ^ { 4 } z ^ { 2 } \\ & - \omega _ { L } ^ { 2 } ( 3 \omega _ { L } ^ { 2 } - 4 \omega _ { z } ^ { 2 } ) \rho ^ { 6 } + \frac { 2 \kappa } { \sqrt { \rho ^ { 2 } } + z ^ { 2 } } \left ( \omega _ { z } ^ { 2 } \rho ^ { 2 } z ^ { 2 } - \omega _ { L } ^ { 2 } \rho ^ { 4 } \right ) + \frac { \kappa ^ { 2 } } { 2 } \, \rho ^ { 2 } + z ^ { 2 }.$$

$$- \omega _ { L } ^ { 2 } ( 3 \omega _ { L } ^ { 2 } - 4 \omega _ { z } ^ { 2 } ) \rho ^ { 6 } + \frac { 2 \kappa } { \sqrt { \rho ^ { 2 } + z ^ { 2 } } } \left ( \omega _ { z } ^ { 2 } \rho ^ { 2 } z ^ { 2 } - \omega _ { L } ^ { 2 } \rho ^ { 4 } \right ) + \frac { \kappa ^ { 2 } } { 2 } \, \frac { \rho ^ { 2 } - z ^ { 2 } } { \rho ^ { 2 } + z ^ { 2 } }.$$

## 6 Non-abelian gauge interactions

The dynamics of a particle with non-abelian gauge interactions is described by Wong's generalization of the Lorentz force law [13]

$$g _ { \mu \nu } \left ( \ddot { x } ^ { \nu } + \Gamma _ { \kappa \lambda } ^ { \ \nu } \dot { x } ^ { \kappa } \dot { x } ^ { \lambda } \right ) = \frac { g } { m } \, t _ { a } F _ { \mu \nu } ^ { a } \dot { x } ^ { \nu }, \quad \dot { t } _ { a } + g f _ { a b } ^ { \ \ c } t _ { c } A _ { \mu } ^ { \ \ b } \dot { x } ^ { \mu }.$$

Here g is the coupling constant, f c ab are the structure constants of the Lie algebra of gauge charges t a , and the non-abelian field-strength tensor is

$$F _ { \mu \nu } ^ { a } = \nabla _ { \mu } A _ { \nu } ^ { a } - \nabla _ { \nu } A _ { \mu } ^ { a } + g f _ { b c } ^ { \ a } A _ { \mu } ^ { b } A _ { \nu } ^ { c }.$$

The theory can be cast in hamiltonian form [10] by introducing covariant momenta π µ and a canonical hamiltonian

$$H = \frac { 1 } { 2 m } \, g ^ { \mu \nu } \pi _ { \mu } \pi _ { \nu },$$

supplemented with brackets

$$& \{ G, K \} = \mathcal { D } _ { \mu } G \frac { \partial K } { \partial \pi _ { \mu } } - \frac { \partial G } { \partial \pi _ { \mu } } \mathcal { D } _ { \mu } K + g t _ { a } F _ { \mu \nu } ^ { a } \frac { \partial G } { \partial \pi _ { \mu } } \frac { \partial K } { \partial \pi _ { \nu } } + f _ { a b } ^ { c } t _ { c } \frac { \partial G } { \partial t _ { a } } \frac { \partial K } { \partial t _ { b } }, \\ & \mathcal { D } _ { \mu } G \equiv \frac { \partial G } { \partial x ^ { \mu } } + \Gamma _ { \mu \nu } ^ { \lambda } \pi _ { \lambda } \frac { \partial G } { \partial \pi _ { \nu } } + g f _ { a b } ^ { c } t _ { c } A _ { \mu } ^ { a } \frac { \partial G } { \partial t _ { b } }. \\ \psi & \quad. \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{$$

In particular these brackets guarantee the Ricci identity and the Lie-algebra of gauge charges:

$$\bar { \Phi } _ { 0 } \{ \pi _ { \mu }, \pi _ { \nu } \} = g t _ { a } \, F _ { \mu \nu } ^ { a }, \quad \{ t _ { a }, t _ { b } \} = f _ { a b } ^ { \ c } t _ { c }.$$

It is now straightforward to derive the condition for the existence of constants of motion:

$$\cdot \\ \{ G, H \} = 0 \quad \Rightarrow \quad \pi ^ { \mu } \mathcal { D } _ { \mu } G = g t _ { a } F ^ { a } _ { \mu \nu } \, \pi ^ { \mu } \frac { \partial G } { \partial \pi _ { \nu } }.$$

In components this gives us a hierarchy of equations similar to (44) for the abelian case:

$$\mathcal { D } _ { \mu } G ^ { ( 0 ) } & = g t _ { a } F _ { \mu } ^ { a \nu } G _ { \nu } ^ { ( 1 ) }, \\ \mathcal { D } _ { \mu } G _ { \nu } ^ { ( 1 ) } & + \mathcal { D } _ { \nu } G _ { \mu } ^ { ( 1 ) } = g t _ { a } \left ( F _ { \mu } ^ { a \lambda } G _ { \lambda \nu } ^ { ( 2 ) } + F _ { \nu } ^ { a \lambda } G _ { \lambda \mu } ^ { ( 2 ) } \right ), \\ \mathcal { D } _ { \mu } G _ { \nu \lambda } ^ { ( 2 ) } & + \mathcal { D } _ { \nu } G _ { \lambda \mu } ^ { ( 2 ) } + \mathcal { D } _ { \lambda } G _ { \mu \nu } ^ { ( 2 ) } = g t _ { a } \left ( F _ { \mu } ^ { a \kappa } G _ { \kappa \lambda \nu } ^ { ( 3 ) } + F _ { \nu } ^ { a \kappa } G _ { \kappa \lambda \mu } ^ { ( 2 ) } + F _ { \lambda } ^ { a \kappa } G _ { \kappa \mu \nu } ^ { ( 3 ) } \right ),$$

...

Example: 2-D SU (2) Yang-Mills point charge

As a simple example consider the dynamics of a non-abelian point charge in a static magnetic SU (2)-field in the euclidean plane. In the 2-dimensional plane such a magnetic field can be written as

$$F _ { i j } ^ { a } = \varepsilon _ { i j } B ^ { a },$$

where a represents an adjoint vector component in 3-dimensional internal SU (2)space. The free Yang-Mills equation then implies

$$\nabla _ { i } B ^ { a } + g \varepsilon ^ { a b c } A _ { i } ^ { b } B ^ { c } = 0 \quad \Rightarrow \quad \nabla _ { i } B ^ { a \, 2 } = 0,$$

and the modulus of B a is constant. Now the direction of B a can be rotated point-wise in the plane by a local gauge transformation, and this freedom can be used to make B a constant. Such a constant B a -field is derived from a vector potential

$$A _ { i } ^ { a } = - \frac { 1 } { 2 } \, \varepsilon _ { i j } x ^ { j } B ^ { a }.$$

Now the 2-dimensional euclidean plane is invariant under translations and rotations, guaranteeing conservation of momentum and angular momentum. Our construction implies, that the angular momentum will obtain a field dependent addition:

$$J = \varepsilon _ { i j } x _ { i } \pi _ { j } + \frac { g } { 2 } \, t _ { a } B ^ { a } x _ { i } ^ { 2 }.$$

In addition, there is also a pair of of rank-2 Killing tensors, generating a conserved Runge-Lenz vector:

$$K _ { i } = x _ { i } \pi _ { j } ^ { 2 } - \pi _ { i } \pi _ { j } x _ { j } + g B ^ { a } t _ { a } \left ( \frac { 1 } { 2 } \, \varepsilon _ { i j } \pi _ { j } x _ { k } ^ { 2 } + x _ { i } x _ { j } \varepsilon _ { j k } \pi _ { k } \right ) + \frac { 1 } { 2 } \, \left ( g B ^ { a } t _ { a } \right ) ^ { 2 } x _ { i } \, x _ { j } ^ { 2 }. \ \ ( 6 4 )$$

## Acknowledgement

The work described in this paper is supported by the Foundation for Fundamental Research of Matter (FOM), as part of the research programme Theoretical Particle Physics in the Era of the LHC .

## References

- [1] E. Noether, Invariante Variationsprobleme , Nachr. D. Knig. Gesellsch. D. Wiss. Zu G¨ttingen, Math-phys. Klasse 1918: o 235
- [2] B. Carter, Global Structure of the Kerr Family of Gravitational Fields , Phys. Rev. 174 (1968), 1559
- [3] B. Carter, Killing tensor quantum numbers and conserved currents in curved space Phys. Rev. D16 (1977), 3395
- [4] J.W. van Holten and R.H. Rietdijk, Symmetries and motions in manifolds , J. Geom. and Physics 11 (1993), 559; arXiv:hep-th/9205074
- [5] R. Penrose, Ann. N. Y. Acad. Sci. 224 (1973), 125;
- [6] R. Floyd, The Dynamics of Kerr Fields PhD. Thesis (London 1973)
- [7] G.W. Gibbons, R.H. Rietdijk and J.W. van Holten, SUSY in the sky , Nucl. Phys. B404 (1993), 42; arXiv:hep-th/9303112
- [8]
- R. Jackiw and N.S. Manton, Symmetries and conservation laws in gauge theories . Ann. Phys. (NY) 127 (1980), 257

- [9] C. Duval and P. Horvathy, Particles with internal structure: the geometry of classical motions and conservation laws , Ann. Phys. (NY) 142 (1982), 10
- [10] J.W. van Holten, Covariant hamiltonian dynamics Phys. Rev. D75 (2007) 025027; arXiv:hep-th/0612216
- [11]
- M. Visinescu, Higher order first integrals of motion in a gauge covariant Hamiltonian framework Mod. Phys. Lett. A25 (2010), 341; arXiv:0910.3474v2 [hep-th]
- M. Cariglia, G.W. Gibbons, J.-W. van Holten, P.A. Horvathy, P. Kos´ ınski and P.-M. Zhang, Killing tensors and canonical geometry in classical dynamics
- [12] Class. Quantum Grav. (in press); arXiv:1401.8195 [hep-th]
- [13] S.K. Wong, Field and particle equations for the classical Yang-Mills field and particles with isotopic spin , Nuov. Cim. 65A (1970), 689
J. W. van Holten
Gauge-covariant extensions of Killing tensors and conservation laws
In classical and quantum mechanical systems on manifolds with gauge-field
fluxes, constants of motion are constructed from gauge-covariant extensions of
Killing vectors and tensors. This construction can be carried out using a
manifestly covariant procedure, in terms of covariant phase space with a
covariant generalization of the Poisson brackets, c.q. quantum commutators.
Some examples of this construction are presented.