The idea of gauge symmetries was introduced in Maxwell's theory. Maxwell's equations are

∇
⋅
E
→
=
ρ
ϵ
0
∇
×
B
→
−
ϵ
0
μ
0
∂
E
→
∂
t
=
μ
0
J
→
∇
×
E
→
+
∂
B
→
∂
t
=
0
∇
⋅
B
→
=
0

where
ρ
is the charge density and
J
→
the current density. The last two equations can be solved by writing fields in terms of a scalar potential,
ϕ
, and a vector potential,
A
→
:

E
→
=
−
∇
ϕ
−
∂
A
→
∂
t
B
→
=
∇
×
A
→
.

The potentials uniquely determine the fields, but the fields do not uniquely determine the potentials - we can make the changes:

ϕ
′
=
ϕ
+
∂
Λ
∂
t
A
→
′
=
A
→
−
∇
Λ

without effecting the electric and magnetic fields, where
Λ
(
x
→
,
t
)
is an arbitrary function of space-time . These are called gauge transformations. There is an elegant relativistic notation: the gauge field is

A
μ
=
(
ϕ
,
A
→
)

and the above gauge transformations read,

A
μ
′
=
A
μ
+
∂
μ
Λ
.

The so-called field strength tensor is introduced,

F
μ
ν
=
∂
μ
A
ν
−
∂
ν
A
μ

which is easily shown to be invariant under gauge transformations. In components,

F
0
i
=
E
i
,
ϵ
i
j
k
F
j
k
=
B
i
.

Maxwell's source-free action is given by:

S
=
−
1
2
∫
d
4
x
(
F
μ
ν
F
μ
ν
)
.

The ability to varying the gauge potential at different points in space and time without changing the physics is called a local invariance. Electromagnetic theory possess the simplest kind of local gauge symmetry called
U
(
1
)
. A theory that displays local gauge invariance is called a gauge theory. In order to formulate other gauge theories we turn the above reasoning inside out. This is the subject of the next section.

We know from quantum mechanics that if we replace the wave-function,
ψ
(
x
)
, describing the electron field by

ψ
′
(
x
)
=
exp
(
i
θ
)
ψ
(
x
)

that it leaves physical predictions unchanged. We consider the imposition of local invariance on the phase of the electron field,

ψ
′
(
x
)
=
Ω
ψ
(
x
)
=
exp
(
i
θ
(
x
)
)
ψ
(
x
)

The problem is that derivatives of
ψ
(
x
)
are not covariant under this transformation:

∂
μ
(
exp
(
i
θ
(
x
)
)
ψ
(
x
)
)
=
Ω
∂
μ
ψ
(
x
)
+
∂
μ
Ω
ψ
(
x
)
.

In order to cancel out the second unwanted term, one introduces a new derivative operator
D
μ
that is covariant. To construct
D
μ
, one introduces a new field, the connection
A
μ
:

D
μ
=
∂
μ
+
i
g
A
μ
(
x
)
.

Then

(
D
μ
ψ
)
′
=
∂
μ
ψ
′
+
i
g
A
μ
′
ψ
′
=
Ω
∂
μ
ψ
+
(
∂
Ω
)
ψ
+
i
g
A
μ
′
Ω
ψ

The term
∂
μ
Ω
is precisely cancelled out by requiring the connection field transforms as

A
μ
′
(
x
)
=
A
μ
(
x
)
+
i
g
[
∂
μ
Ω
(
x
)
]
Ω
−
1
(
x
)
E
q
1.
.

We then have that

(
D
μ
ψ
)
′
=
Ω
D
μ
ψ
.

Note that
E
q
1
is equivalent to

A
μ
′
(
x
)
=
A
μ
(
x
)
+
1
g
∂
μ
θ
(
x
)

which looks the same as a gauge transformation of the gauge potential of Maxwell's theory. It is possible to construct an invariant action for the connection field itself. We want an action that only has two derivatives (since actions with higher derivatives are not unitary). Define the quantity:

F
μ
ν
=
−
i
g
[
D
μ
,
D
μ
]
=
−
i
g
[
∂
μ
+
i
g
A
μ
(
x
)
,
∂
ν
+
i
g
A
ν
(
x
)
]

=
−
i
g
(
[
∂
μ
,
∂
ν
]
+
i
g
(
∂
μ
A
ν
−
∂
ν
A
μ
)
−
g
2
[
A
μ
,
A
ν
]
)

=
∂
μ
A
ν
−
∂
ν
A
ν
.

The unique action with only two derivatives is given by:

S
=
−
1
2
∫
d
4
x
(
F
μ
ν
F
μ
ν
)
.

Therefore, one can derive electromagnetic theory from arguments based solely on symmetry.

We now generalize the above reasoning to general gauge groups. One begins with the generators of some Lie algebra:

[
T
i
,
T
j
]
=
i
f
i
j
k
T
k
Let there be a fermion field that transforms as

Ψ
′
↦
Ω
^
(
x
)
Ψ
(
x
)
=
exp
(
i
θ
i
(
x
)
T
i
)
Ψ
(
x
)
Again the derivatives of
Ψ
(
x
)
are not covariant under this transformation. We introduce a covariant derivative

D
μ
=
I
∂
μ
+
i
g
A
μ
(
x
)
with connection field given by

A
μ
(
x
)
=
A
μ
i
(
x
)
T
i
We require that
A
μ
(
x
)
transforms as:

A
μ
′
(
x
)
=
Ω
^
A
μ
(
x
)
Ω
^
−
1
+
i
g
Ω
^
(
∂
μ
Ω
^
−
1
)
.

We define the field strength operator

F
μ
ν
=
−
i
g
[
D
μ
,
D
ν
]
=
∂
μ
A
ν
−
∂
ν
A
μ
+
i
g
[
A
μ
,
A
ν
]
=
(
∂
μ
A
ν
i
−
∂
ν
A
μ
i
+
g
f
i
j
k
A
μ
j
A
ν
k
)
T
i
.

As
D
μ
is covariant, this means that the
F
μ
ν
i
tensor is also covariant:

F
μ
ν
↦
F
μ
ν
′
=
Ω
^
F
μ
ν
Ω
^
−
1
Note that
F
μ
ν
is only invariant under gauge transformations if
Ω
^
is a scalar, that is, only in the case of electromagnetism.

We can now construct an invariant action out of this tensor. Again we want an action that only has two derivatives. The simplest choice is the trace of the commutator:

Tr
(
Ω
^
F
μ
ν
Ω
^
−
1
Ω
^
F
μ
ν
Ω
^
−
1
)
=
Tr
(
F
μ
ν
F
μ
ν
)
The unique action with only two derivatives is given by:

S
=
−
1
2
∫
d
4
x
T
r
(
F
μ
ν
F
μ
ν
)
=
−
1
2
∫
d
4
x
Tr
(
F
μ
ν
i
T
j
F
j
μ
ν
T
j
)
This is the action for Yang-mills theory.

We consider a change of representation in the quantum Maxwell gauge theory. The idea is to introduce a basis of states labeled by loops
∣
γ
⟩
whose inner product with the connection states is given by

⟨
A
∣
γ
⟩
=
W
(
γ
)
=
exp
[
i
e
∫
γ
d
y
α
A
α
(
y
)
]
The loop functional
W
(
γ
)
is the Wilson loop for the abelian
U
(
1
)
case.

We consider for simplicity (and because later we will see this is the relevant gauge group in LQG) an
S
U
(
2
)
Yang–Mills theory in four dimensions. The field variable of the continuous theory is an
S
U
(
2
)
connection (or gauge potential)
A
μ
i
(
x
)
, where
i
is an index in the Lie algebra of
S
U
(
2
)
. We can write for this field

A
μ
(
x
)
=
A
μ
i
(
x
)
τ
i
where
τ
i
are the
s
u
(
2
)
generators, that is the Pauli matrices multiplied by
i
/
2
. note that unlike with Maxwell's theory, the connections
A
μ
(
x
)
are matrix-valued and don't commute, that is they are non-Abelian gauge theories. We must take this into account when defining the corresponding version of the holonomy for
S
U
(
2
)
Yang–Mills theory.

We first describe the quantum theory in terms of connection variable.

In the connection representation the configuration variable is
A
a
i
and its conjugate momentum is the (densitized) triad
E
~
i
a
. It is most natural to consider wavefunctions
Ψ
(
A
a
i
)
. This is known as the connection representation. The canonical variables get promoted to quantum operators:

A
^
a
i
Ψ
[
A
]
=
A
a
i
Ψ
[
A
]
(analogous to the position representation
q
^
ψ
(
q
)
=
q
ψ
(
q
)
) and the triads are functional derivatives,

E
~
^
a
i
Ψ
[
A
]
=
−
i
δ
Ψ
[
A
]
δ
A
a
i
(analogous to
p
^
ψ
(
q
)
=
−
i
d
ψ
(
q
)
d
q
)

Let us return to the classical Yang–Mills theory. It is possible to encode the gauge invariant information of the theory in terms of `loop-like' variables.

We need the notion of a holonomy. A holonomy is a measure of how much the initial and final values of a spinor or vector differ after parallel transport around a closed loop; it is denoted

h
γ
[
A
]
Knowledge of the holonomies is equivalent to knowledge of the connection, up to gauge equivalence. Holonomies can also be associated with an edge; under a Gauss Law these transform as

(
h
e
′
)
α
β
=
U
α
γ
−
1
(
x
)
(
h
e
)
γ
σ
U
σ
β
(
y
)
.
For a closed loop
x
=
y
if we take the trace of this, that is, putting
α
=
β
and summing we obtain

(
h
e
′
)
α
α
=
U
α
γ
−
1
(
x
)
(
h
e
)
γ
σ
U
σ
α
(
x
)
=
[
U
σ
α
(
x
)
U
α
γ
−
1
(
x
)
]
(
h
e
)
γ
σ
=
δ
σ
γ
(
h
e
)
γ
σ
=
(
h
e
)
γ
γ
or

Tr
h
γ
′
=
Tr
h
γ
.
Thus the trace of an holonomy around a closed loop is gauge invariant. It is denoted

W
γ
[
A
]
and is called a Wilson loop. The explicit form of the holonomy is

h
γ
[
A
]
=
P
exp
{
−
∫
γ
0
γ
1
d
s
γ
˙
a
A
a
i
(
γ
(
s
)
)
T
i
}
where
γ
is the curve along which the holonomy is evaluated, and
s
is a parameter along the curve,
P
denotes path ordering meaning factors for smaller values of
s
appear to the left, and
T
i
are matrices that satisfy the
s
u
(
2
)
algebra

[
T
i
,
T
j
]
=
2
i
ϵ
i
j
k
T
k
.
The Pauli matrices satisfy the above relation. It turns out that there are infinitely many more examples of sets of matrices that satisfy these relations, where each set comprises
(
N
+
1
)
×
(
N
+
1
)
matrices with
N
=
1
,
2
,
3
,
…
, and where none of these can be thought to `decompose' into two or more examples of lower dimension. They are called different irreducible representations of the
s
u
(
2
)
algebra. The most fundamental representation being the Pauli matrices. The holonomy is labelled by a half integer
N
/
2
according to the irreducible representation used.

An important theorem about Yang–Mills gauge theories is Giles' theorem, according to which if one gives the trace of the holonomy of a connection for all possible loops on a manifold one can, in principle, reconstruct all the gauge invariant information of the connection. That is, Wilson loops constitute a basis of gauge invariant functions of the connection. This key result is the basis for the loop representation for gauge theories and gravity.

The use of Wilson loops explicitly solves the Gauss gauge constraint. As Wilson loops form a basis we can formally expand any Gauss gauge invariant function as,

Ψ
[
A
]
=
∑
γ
Ψ
[
γ
]
W
γ
[
A
]
.

This is called the loop transform. We can see the analogy with going to the momentum representation in quantum mechanics. There one has a basis of states
exp
(
i
k
x
)
labelled by a number
k
and one expands

ψ
[
x
]
=
∫
d
k
ψ
(
k
)
exp
(
i
k
x
)
.
and works with the coefficients of the expansion
ψ
(
k
)
.

The inverse loop transform is defined by

Ψ
[
γ
]
=
∫
[
d
A
]
Ψ
[
A
]
W
γ
[
A
]
.
This defines the loop representation. Given an operator
O
^
in the connection representation,

Φ
[
A
]
=
O
^
Ψ
[
A
]
,
Eq 1
one should define the corresponding operator
O
^
′
on
Ψ
[
γ
]
in the loop representation via,

Φ
[
γ
]
=
O
^
′
Ψ
[
γ
]
,
Eq 2
where
Φ
[
γ
]
is defined by the usual inverse loop transform,

Φ
[
γ
]
=
∫
[
d
A
]
Φ
[
A
]
W
γ
[
A
]
.
Eq 3

A transformation formula giving the action of the operator
O
^
′
on
Ψ
[
γ
]
in terms of the action of the operator
O
^
on
Ψ
[
A
]
is then obtained by equating the R.H.S. of
E
q
2
with the R.H.S. of
E
q
3
with
E
q
1
substituted into
E
q
3
, namely

O
^
′
Ψ
[
γ
]
=
∫
[
d
A
]
W
γ
[
A
]
O
^
Ψ
[
A
]
,
or

O
^
′
Ψ
[
γ
]
=
∫
[
d
A
]
(
O
^
†
W
γ
[
A
]
)
Ψ
[
A
]
,
where by
O
^
†
we mean the operator
O
^
but with the reverse factor ordering (remember from simple quantum mechanics where the product of operators is reversed under conjugation). We evaluate the action of this operator on the Wilson loop as a calculation in the connection representation and rearranging the result as a manipulation purely in terms of loops (one should remember that when considering the action on the Wilson loop one should choose the operator one wishes to transform with the opposite factor ordering to the one chosen for its action on wavefunctions
Ψ
[
A
]
).

The introduction of Ashtekar variables cast general relativity in the same language as gauge theories. It was in particular the inability to have good control over the space of solutions to the Gauss' law and spatial diffeomorphism constraints that led Rovelli and Smolin to consider a new representation – the loop representation.

To handle the spatial diffeomorphism constraint we need to go over to the loop representation. The above reasoning gives the physical meaning of the operator
O
^
′
. For example, if
O
^
†
corresponded to a spatial diffeomorphism, then this can be thought of as keeping the connection field
A
of
W
γ
[
A
]
where it is while performing a spatial diffeomorphism on
γ
instead. Therefore, the meaning of
O
^
′
is a spatial diffeomorphism on
γ
, the argument of
Ψ
[
γ
]
.

In the loop representation we can then solve the spatial diffeomorphism constraint by considering functions of loops
Ψ
[
γ
]
that are invariant under spatial diffeomorphisms of the loop
γ
. That is, we construct what mathematicians call knot invariants. This opened up an unexpected connection between knot theory and quantum gravity.

The easiest geometric quantity is the area. Let us choose coordinates so that the surface
Σ
is characterized by
x
3
=
0
. The area of small parallelogram of the surface
Σ
is the product of length of each side times
sin
θ
where
θ
is the angle between the sides. Say one edge is given by the vector
u
→
and the other by
v
→
then,

A
=
∥
u
→
∥
∥
v
→
∥
sin
θ
=
∥
u
→
∥
2
∥
v
→
∥
2
(
1
−
cos
2
θ
)
=
∥
u
→
∥
2
∥
v
→
∥
2
−
(
u
→
⋅
v
→
)
2
From this we get the area of the surface
Σ
to be given by

A
Σ
=
∫
Σ
d
x
1
d
x
2
det
q
(
2
)
where
det
q
(
2
)
=
q
11
q
22
−
q
12
2
and is the determinant of the metric induced on
Σ
. This can be rewritten as

det
q
(
2
)
=
ϵ
3
a
b
ϵ
3
c
d
q
a
c
q
b
c
2
.
The standard formula for an inverse matrix is

q
a
b
=
ϵ
b
c
d
ϵ
a
e
f
q
c
e
q
d
f
2
!
det
(
q
)
Note the similarity between this and the expression for
det
q
(
2
)
. But in Ashtekar variables we have
E
~
i
a
E
~
b
i
=
det
(
q
)
q
a
b
. Therefore,

A
Σ
=
∫
Σ
d
x
1
d
x
2
E
~
i
3
E
~
3
i
.
According to the rules of canonical quantization we should promote the triads
E
~
i
3
to quantum operators,

E
~
^
i
3
∼
δ
δ
A
3
i
.
It turns out that the area
A
Σ
can be promoted to a well defined quantum operator despite the fact that we are dealing with product of two functional derivatives and worse we have a square-root to contend with as well. Putting
N
=
2
J
, we talk of being in the *J*-th representation. We note that
∑
i
T
i
T
i
=
J
(
J
+
1
)
1
. This quantity is important in the final formula for the area spectrum. We simply state the result below,

A
^
Σ
W
γ
[
A
]
=
8
π
ℓ
P
l
a
n
c
k
2
β
∑
I
j
I
(
j
I
+
1
)
W
γ
[
A
]
where the sum is over all edges
I
of the Wilson loop that pierce the surface
Σ
.

The formula for the volume of a region
R
is given by

V
=
∫
R
d
3
x
det
(
q
)
=
1
6
∫
R
d
x
3
ϵ
a
b
c
ϵ
i
j
k
E
~
i
a
E
~
j
b
E
~
k
c
.
The quantization of the volume proceeds the same way as with the area. As we take the derivative, and each time we do so we bring down the tangent vector
γ
˙
a
, when the volume operator acts on non-intersecting Wilson loops the result vanishes. Quantum states with non-zero volume must therefore involve intersections. Given that the anti-symmetric summation is taken over in the formula for the volume we would need at least intersections with three non-coplanar lines. Actually it turns out that one needs at least four-valent vertices for the volume operator to be non-vanishing.

We now consider Wilson loops with intersections. We assume the real representation where the gauge group is
S
U
(
2
)
. Wilson loops are an over complete basis as there are identities relating different Wilson loops. These come about from the fact that Wilson loops are based on matrices (the holonomy) and these matrices satisfy identities, the so-called Mandelstam identities (see Mandelstam variables). Given any two
S
U
(
2
)
matrices
A
and
B
it is easy to check that,

Tr
(
A
)
Tr
(
B
)
=
Tr
(
A
B
)
+
Tr
(
A
B
−
1
)
.
This implies that given two loops
γ
and
η
that intersect, we will have,

W
γ
[
A
]
W
η
[
A
]
=
W
γ
∘
η
[
A
]
+
W
γ
∘
η
−
1
[
A
]
where by
η
−
1
we mean the loop
η
traversed in the opposite direction and
γ
∘
η
means the loop obtained by going around the loop
γ
and then along
η
. See figure below. This is called a Mandelstam identity of the second kind. There is the Mandelstam identity of the first kind
W
(
γ
1
∘
γ
2
)
=
W
(
γ
2
∘
γ
1
)
. Spin networks are certain linear combinations of intersecting Wilson loops designed to address the over-completeness introduced by the Mandelstam identities.

In fact spin networks constitute a basis for all gauge invariant functions which minimize the degree of over-completeness of the loop basis, and for trivalent intersections eliminate it entirely.

As mentioned above the holonomy tells you how to propagate test spin half particles. A spin network state assigns an amplitude to a set of spin half particles tracing out a path in space, merging and splitting. These are described by spin networks
γ
: the edges are labelled by spins together with `intertwiners' at the vertices which are prescription for how to sum over different ways the spins are rerouted. The sum over rerouting are chosen as such to make the form of the intertwiner invariant under Gauss gauge transformations.

Theorems establishing the uniqueness of the loop representation as defined by Ashtekar et al. (i.e. a certain concrete realization of a Hilbert space and associated operators reproducing the correct loop algebra – the realization that everybody was using) have been given by two groups (Lewandowski, Okolow, Sahlmann and Thiemann) and (Christian Fleischhack). Before this result was established it was not known whether there could be other examples of Hilbert spaces with operators invoking the same loop algebra, other realizations, not equivalent to the one that had been used so far.

A common method of describing a knot (or link, which are knots of several components entangled with each other) is to consider its projected image onto a plane called a knot diagram. Any given knot (or link) can be drawn in many different ways using a knot diagram. Therefore, a fundamental problem in knot theory is determining when two descriptions represent the same knot. Given a knot diagram, one tries to find a way to assign a knot invariant to it, sometimes a polynomial – called a knot polynomial. Two knot diagrams with different polynomials generated by the same procedure necessarily correspond to different knots. However, if the polynomials are the same, it may not mean that they correspond to the same knot. The better a polynomial is at distinguishing knots the more powerful it is.

In 1984, Jones announced the discovery of a new link invariant, which soon led to a bewildering profusion of generalizations. He had found a new knot polynomial, the Jones polynomial. Specifically, it is an invariant of an oriented knot or link which assigns to each oriented knot or link a polynomial with integer coefficients.

In the late 1980s, Witten coined the term topological quantum field theory for a certain type of physical theory in which the expectation values of observable quantities are invariant under diffeomorphisms.

Witten gave a heuristic derivation of the Jones polynomial and its generalizations from Chern–Simons theory. The basic idea is simply that the vacuum expectation values of Wilson loops in Chern–Simons theory are link invariants because of the diffeomorphism-invariance of the theory. To calculate these expectation values, however, Witten needed to use the relation between Chern–Simons theory and a conformal field theory known as the Wess–Zumino–Witten model (or the WZW model).