Next Article in Journal
Numerical Equilibrium Configurations and Quadrupole Moments of Post-Merger Differentially Rotating Relativistic Stars
Next Article in Special Issue
Weak Gravitation in the 4+1 Formalism
Previous Article in Journal
The NuGrid AGB Evolution and Nucleosynthesis Data Set
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Singular Lagrangians, Constrained Hamiltonian Systems and Gauge Invariance: An Example of the Dirac–Bergmann Algorithm

Department of Physics, North Carolina State University, Raleigh, NC 27695, USA
Universe 2022, 8(3), 171; https://doi.org/10.3390/universe8030171
Submission received: 16 February 2022 / Revised: 1 March 2022 / Accepted: 5 March 2022 / Published: 9 March 2022
(This article belongs to the Special Issue Numerical Relativity)

Abstract

:
The Dirac–Bergmann algorithm is a recipe for converting a theory with a singular Lagrangian into a constrained Hamiltonian system. Constrained Hamiltonian systems include gauge theories—general relativity, electromagnetism, Yang–Mills, string theory, etc. The Dirac–Bergmann algorithm is elegant but at the same time rather complicated. It consists of a large number of logical steps linked together by a subtle chain of reasoning. Examples of the Dirac–Bergmann algorithm found in the literature are designed to isolate and illustrate just one or two of those logical steps. In this paper, I analyze a finite-dimensional system that exhibits all of the major steps in the algorithm. The system includes primary and secondary constraints, first and second class constraints, restrictions on Lagrange multipliers, and both physical and gauge degrees of freedom. This relatively simple system provides a platform for discussing the Dirac conjecture, constructing Dirac brackets, and applying gauge conditions.

1. Introduction

In the early 1950’s, Dirac and Bergmann independently developed the Hamiltonian formalism for systems with singular Lagrangians [1,2,3,4,5,6,7,8,9,10]. These systems, often called “constrained Hamiltonian systems”, include gauge theories. Gauge freedom is more clearly and more completely displayed in the Hamiltonian setting, with the generators of gauge transformations expressed as functions on phase space. Historically, the main motivation for casting gauge theories in Hamiltonian form was to facilitate their canonical quantization. Dirac and Bergmann were primarily motivated by the prospect of developing a quantum theory of gravity based on a Hamiltonian formulation of general relativity.
Textbook treatments of Lagrangian and Hamiltonian mechanics invariably assume that the Lagrangian L ( q , q ˙ ) is nonsingular; that is, that the matrix of second derivatives of L ( q , q ˙ ) with respect to the velocities is invertible. In classical mechanics, the nonsingular case appears to be sufficient to cover problems of physical interest. However, one might argue that textbooks avoid certain physically interesting problems simply because their Lagrangians are singular.
In field theory, the issue of singular Lagrangians and gauge freedom cannot be avoided. Nearly every field theory of physical interest—electrodynamics, Yang–Mills theory, general relativity, relativistic string theory—has gauge freedom.
The Dirac–Bergmann algorithm transforms a singular Lagrangian system into a Hamiltonian system. The formalism is elegant but at the same time rather complex. It consists of a large number of logical steps, linked together by a chain of reasoning that can be difficult to keep straight. Of course, there are many examples in the literature in which the Dirac–Bergmann algorithm is applied, converting a singular Lagrangian into Hamiltonian form. However, to my knowledge, all of these examples are designed to illustrate just one or two of the logical steps in the algorithm. The student of the subject is faced with the task of linking these examples together to create a complete picture of the algorithm.
For those who learn by example, what is needed is a single example that illustrates all of the major logical steps in the Dirac–Bergmann algorithm and shows how these steps are linked together. Such a “complete” example is not easy to identify because there is no obvious way to predict, starting with a particular Lagrangian, which of the steps in the algorithm will be needed.
The system analyzed in this paper is defined by the Lagrangian
L ( q , q ˙ ) = 1 2 ( q 1 + q ˙ 2 + q ˙ 3 ) 2 + ( q ˙ 4 q ˙ 2 ) 2 + ( q 1 + 2 q 2 ) ( q 1 + 2 q 4 )
where the dot denotes a time derivative. The matrix of second derivatives with respect to the velocities is
2 L q ˙ i q ˙ j = 0 0 0 0 0 2 1 1 0 1 1 0 0 1 0 1 .
This matrix is singular; it has rank 2.
As we will see, the system defined by the Lagrangian (1) is relatively complete.1 It contains both primary and secondary constraints, both first and second class constraints, and restrictions on the Lagrange multipliers. The first class constraints for this system are not all primary; this allows us to address the Dirac conjecture. The second class constraints can be eliminated by introducing Dirac brackets. Finally, this system contains both physical and gauge degrees of freedom. The gauge freedom can be eliminated with suitable gauge conditions.
One characteristic of any complete example such as (1) is that the configuration space, the space of q’s, must be at least four-dimensional. Here is why: The number of physical degrees of freedom is equal to the dimension of the configuration space, minus the number of first class constraints, minus half the number of second class constraints. If the example is to have at least one physical degree of freedom, at least two first class constraints (one primary and one secondary), and at least two second class constraints (the number of second class constraints must be even), then the configuration space must be at least four–dimensional.
The study of constrained Hamiltonian systems predates Dirac and Bergmann with earlier work by Rosenfeld [12,13]. Like Rosenfeld, Bergmann and his collaborators [1,2,4,6,7,14] were focused on field theories such as general relativity that are covariant with respect to general four-dimensional coordinate transformations. (For a historical review, see [15].) Dirac took a more basic approach to the problem by considering a generic singular Lagrangian [3,5,8,9,10]. He developed the algorithm for the case of systems with a finite number of degrees of freedom. His view was that the generalization to field theory, with an infinite number of degrees of freedom, would be “merely a formal matter”.
Although the typical starting point for the Dirac–Bergmann analysis is a singular Lagrangian, not all gauge theories can be expressed in terms of a Lagrangian that depends only on the q’s. General relativity is one such example. The Einstein–Hilbert action is a functional of Lagrange multipliers (the lapse function and shift vector) as well as the configuration variables (the spatial metric). Nevertheless, the Dirac–Bergmann algorithm provides the foundation for our interpretation of general relativity as a constrained Hamiltonian system.
In this paper, I apply the Dirac–Bergmann algorithm to the Lagrangian (1), following closely the general treatment given by Henneaux and Teitelboim [11]. In turn, the account of Henneaux and Teitelboim closely follows Dirac’s 1964 Lectures on Quantum Mechanics [10]. Presentations of the Dirac–Bergmann algorithm can also be found in books by Hanson, Regge and Teitelboim [16], Sundermeyer [17], Rothe and Rothe [18], and Lusanna [19].
Throughout the paper, I attempt to explain the reasoning behind the logical steps of the Dirac–Bergmann algorithm, but avoid general proofs. The reader is referred to references [11,16,17,18,19] for more details.
We begin in Section 2 with a derivation of Lagrange’s equations for the singular Lagrangian (1). The general solution is derived, and in Section 3, we discuss the gauge freedom at the Lagrangian level. We begin construction of the Hamiltonian theory in Section 4 with a derivation of the primary constraints and canonical Hamiltonian. In Section 5, we introduce the primary Hamiltonian and the primary action. Section 6 is devoted to a discussion of the initial value problem and the need to go beyond the primary Hamiltonian. In Section 7, we apply Dirac’s consistency conditions to derive the secondary constraints and restrictions on the Lagrange multipliers. The concept of weak equality is introduced in Section 8, along with a formal analysis of the restrictions on the Lagrange multipliers. The total Hamiltonian is computed in Section 9, and in Section 10 we sort the constraints into first and second class. The first class Hamiltonian and gauge generators are identified in Section 11, where we also introduce the Dirac conjecture. In Section 12 and Section 13, we define the extended Hamiltonian and extended action. Dirac brackets are used in Section 14 to eliminate the second class constraints, which yields a partially reduced Hamiltonian. The corresponding partially reduced action is derived in Section 15, and in Section 16 we eliminate the momenta to obtain a partially reduced Lagrangian. Gauge conditions are introduced in Section 17 and Dirac brackets are used to eliminate the constraints and gauge conditions. This yields a fully reduced Hamiltonian. The fully reduced action is derived in Section 18. Finally, Section 19 contains a summary of the Dirac–Bergmann algorithm and a discussion of the Einstein–Hilbert action for general relativity.

2. Lagrangian Analysis

The action is the integral of the Lagrangian (1):
S [ q ] = 0 T d t L ( q , q ˙ ) .
The notation S [ q ] indicates that S is a functional of the complete set of coordinates, q i = { q 1 , q 2 , q 3 , q 4 } . The equations of motion are obtained by extremizing the action. For this example, we are not concerned with boundary conditions and integrate by parts freely. Lagrange’s equations are
0 = δ S δ q 1 = q ˙ 2 + q ˙ 3 + 2 q 1 + q 2 + q 4 ,
0 = δ S δ q 2 = 2 q ¨ 2 q ¨ 3 + q ¨ 4 q ˙ 1 + q 1 + 2 q 4 ,
0 = δ S δ q 3 = q ¨ 2 q ¨ 3 q ˙ 1 ,
0 = δ S δ q 4 = q ¨ 2 q ¨ 4 + q 1 + 2 q 2 .
We can rewrite these as follows. First, add Equation (4b,d), then subtract Equation (4c). This gives
q 1 + q 2 + q 4 = 0 .
Next, subtract this result from Equation (4a) to obtain
q ˙ 2 + q ˙ 3 + q 1 = 0 .
The time derivative of this equation yields Equation (4c). Finally, we find the result
q ¨ 4 q ¨ 2 = q 2 q 4 .
by solving Equation (5a) for q 1 and using the equation of motion (4d).
Equation (5) are equivalent to Lagrange’s Equation (4). In particular, Equation (4a) is the sum of Equation (5a,b); Equation (4b) is the sum of Equation (5a,c) with the time derivative of (5b) subtracted; Equation (4c) is the negative of the time derivative of (5b); Equation (4d) is obtained by subtracting (5c) from (5a).
The equations of motion for this simple linear system are easily solved. Note that the combination q 4 q 2 is determined by (5c) along with initial or boundary data; thus, we have
q 4 ( t ) q 2 ( t ) = A sin t + B cos t ,
where A and B are constants. Now, Equation (5a) gives
q 1 ( t ) + 2 q 2 ( t ) = A sin t B cos t .
If we knew q 2 ( t ) , we could solve the previous two equations for q 1 ( t ) and q 4 ( t ) , then integrate Equation (5b) to obtain q 3 ( t ) . Clearly, we do not have enough information to fully determine each of the q’s as functions of time. One of the q’s must remain undetermined. For example, let us choose q 2 arbitrarily by setting q 2 ( t ) = Ψ ( t ) / 2 for some function Ψ ( t ) . We can then use the equations above to solve for q 1 , q 3 and q 4 :
q 1 ( t ) = A sin t B cos t + Ψ ( t ) ,
q 2 ( t ) = Ψ ( t ) / 2 ,
q 3 ( t ) = A cos t + B sin t + Ψ ( t ) / 2 0 t d s Ψ ( s ) + C ,
q 4 ( t ) = A sin t + B cos t Ψ ( t ) / 2 ,
where C is an integration constant. This is the general solution of the equations of motion.

3. Gauge Invariance

The undetermined function Ψ ( t ) that appears in the general solution (8) can be freely specified. This is the gauge freedom of the theory. We can express the gauge freedom in another way: the Lagrangian (1) and the equations of motion (either (4) or (5)) are invariant under the replacements
q 1 ( t ) q 1 ( t ) + Ψ ( t ) ,
q 2 ( t ) q 2 ( t ) Ψ ( t ) / 2 ,
q 3 ( t ) q 3 ( t ) + Ψ ( t ) / 2 0 t d s Ψ ( s ) ,
q 4 ( t ) q 4 ( t ) Ψ ( t ) / 2 ,
where Ψ ( t ) is an arbitrary function of time.
Although each configuration of the system (that is, each set of q values) corresponds to a specific physical state of the system, the converse is not true. Because of the gauge freedom, there are many sets of q’s that describe one and the same physical state.
Let us examine the gauge freedom more closely, in anticipation of the Hamiltonian description of evolution. To begin, choose the gauge Ψ ( t ) = 0 and consider the general solution (8). This solution describes the evolution of the system from initial data
q 1 ( 0 ) = B , q ˙ 1 ( 0 ) = A ,
q 2 ( 0 ) = 0 , q ˙ 2 ( 0 ) = 0 ,
q 3 ( 0 ) = C A , q ˙ 3 ( 0 ) = B ,
q 4 ( 0 ) = B , q ˙ 4 ( 0 ) = A .
The configuration at some arbitrary final time t = T is
q 1 ( T ) = A sin T B cos T ,
q 2 ( T ) = 0 ,
q 3 ( T ) = C A cos T + B sin T ,
q 4 ( T ) = A sin T + B cos T .
This configuration corresponds to a particular state of the physical system.
We can choose a different gauge in Equation (8). As long as the new gauge satisfies Ψ ( 0 ) = Ψ ˙ ( 0 ) = 0 , the solution will describe evolution from the same initial data (10). For example, with
Ψ ( t ) = ( π 2 4 ) ϵ 8 cos ( π t / T ) 1 + π ϵ 4 π t / T sin ( π t / T )
where ϵ = const , the configuration at t = T is
q 1 ( T ) = A sin T B cos T + ϵ ,
q 2 ( T ) = ϵ / 2 ,
q 3 ( T ) = C A cos T + B sin T + ϵ / 2 ,
q 4 ( T ) = A sin T + B cos T ϵ / 2 .
The configurations (11) and (13) represent the same physical state of the system, since they evolve from the same initial data.
We can express this result more compactly as
δ q 1 = ϵ ,
δ q 2 = ϵ / 2 ,
δ q 3 = ϵ / 2 ,
δ q 4 = ϵ / 2 .
Here, δ q i denotes the change in q i at the generic time T, due to the change in gauge function Ψ ( t ) .
Here is another example. With
Ψ ( t ) = π 2 ϵ 4 T cos ( π t / T ) 1 + π ϵ 2 T π t / T sin ( π t / T )
we obtain a configuration that differs from Equation (11) by
δ q 1 = 0 ,
δ q 2 = 0 ,
δ q 3 = ϵ ,
δ q 4 = 0 .
This configuration is also evolved from the initial data (10), and represents the same physical state as the configurations (11) and (13).
Although the gauge transformation (9) contains a single arbitrary function of time, the gauge invariance naturally splits into two types. The first consists of variations subject to δ q 2 = δ q 3 = δ q 4 = δ q 1 / 2 . The second consists of arbitrary variations in δ q 3 , with δ q 1 = δ q 2 = δ q 4 = 0 . This apparent “doubling” of the gauge freedom arises because the solution (9c) for q 3 ( t ) (unlike the other variables) includes the integral of Ψ ( t ) . There is enough freedom of choice in Ψ ( t ) to allow variations in q 3 that are independent of the variations among the other variables. Both types of gauge transformations leave the physical state of the system unchanged.
The consequences of gauge invariance are most clearly expressed in the Hamiltonian formalism. The extended Hamiltonian defined in Section 12 includes phase space generators for both types of gauge transformations.

4. Primary Constraints and the Canonical Hamiltonian

We now begin construction of the Hamiltonian description of the system. The conjugate momenta are defined as usual by p i = L / q ˙ i . For the Lagrangian (1), we have
p 1 = 0
p 2 = 2 q ˙ 2 + q ˙ 3 q ˙ 4 + q 1
p 3 = q ˙ 2 + q ˙ 3 + q 1
p 4 = q ˙ 4 q ˙ 2
Because the Lagrangian is singular, the matrix of second derivatives 2 L / q ˙ i q ˙ j is not invertible and we cannot solve Equation (17) for the velocities as functions of the coordinates and momenta. The definitions (17) yield two primary constraints,
ϕ 1 p 1 = 0 ,
ϕ 2 p 2 p 3 + p 4 = 0 ,
that restrict the phase space variables p i , q i . We will denote these constraints collectively by ϕ a , where a = 1 , 2 . Note that in this simple example, the primary constraints are independent of the q’s.
There is freedom in choosing how the constraints are written. For example, we could replace the ϕ ’s above with
ϕ ˜ 1 p 1 ,
ϕ ˜ 2 p 1 + p 2 p 3 + p 4 .
In fact, any choice for the constraints is allowed, as long as they follow from the definitions p i = L / q ˙ i and satisfy the regularity conditions. These conditions state that, roughly speaking, the constraints should have nonzero gradients on the constraint surface. More precisely, the Jacobian matrix formed from the derivatives of the constraints with respect to the p’s and q’s should have maximal rank on the constraint subspace [10,11]. For both choices, (18) and (19), the rank of the Jacobian matrix is 2. On the other hand, the set
ϕ ¯ 1 p 1 2 ,
ϕ ¯ 2 p 2 p 3 + p 4 ,
is not permissible because the gradient of ϕ ¯ 1 vanishes on the constraint surface where p 1 = 0 . Correspondingly, the rank of the Jacobian matrix is less than 2 on the constraint surface.
The next step in constructing the Hamiltonian formalism is to compute the canonical Hamiltonian. The canonical Hamiltonian H C is defined from the usual prescription by writing p i q ˙ i L ( q , q ˙ ) in terms of p’s and q’s. Although we cannot solve for all of the q ˙ ’s in terms of p’s and q’s, it can be shown that the combination p i q ˙ i L ( q , q ˙ ) depends only on the phase space variables [10,11]. For our example problem, the canonical Hamiltonian is
H C = 1 2 p 3 2 + p 4 2 2 p 3 q 1 ( q 1 + 2 q 2 ) ( q 1 + 2 q 4 ) .
Note that H C is ambiguous. For example, we could use the primary constraint (18b) to replace the term 2 p 3 q 1 with 2 ( p 2 + p 4 ) q 1 .

5. Primary Hamiltonian and the Primary Action

The primary Hamiltonian H P is obtained from the canonical Hamiltonian H C by adding the primary constraints with Lagrange multipliers,
H P = H C + λ a ϕ a ,
where a sum over the repeated index a is implied. The primary action is built from the primary Hamiltonian in the usual way: S P [ q , p , λ ] = 0 T d t { p i q ˙ i H P } . Explicitly, we have
S P [ q , p , λ ] = 0 T d t { p i q ˙ i 1 2 p 3 2 + p 4 2 2 p 3 q 1 ( q 1 + 2 q 2 ) ( q 1 + 2 q 4 ) λ 1 p 1 λ 2 ( p 2 p 3 + p 4 ) } .
The primary action is a functional of the complete set of phase space coordinates, q i , p i , as well as the Lagrange multipliers λ 1 and λ 2 .
The equations of motion are obtained by extremizing the primary action S P . Extremization with respect to the momenta p i gives
q ˙ 1 = λ 1 ,
q ˙ 2 = λ 2 ,
q ˙ 3 = p 3 q 1 λ 2 ,
q ˙ 4 = p 4 + λ 2 ,
while extremization with respect to the coordinates q i yields
p ˙ 1 = p 3 + q 1 + q 2 + q 4 ,
p ˙ 2 = q 1 + 2 q 4 ,
p ˙ 3 = 0 ,
p ˙ 4 = q 1 + 2 q 2 .
Extremizing the action S P with respect to the Lagrange multipliers gives the constraints,
ϕ 1 p 1 = 0 ,
ϕ 2 p 2 p 3 + p 4 = 0 .
These equations of motion (24) are equivalent to Lagrange’s Equation (4). To show this, we first solve Equation (24c,d,i,j) for the momenta to obtain
p 1 = 0 ,
p 2 = q ˙ 3 q ˙ 4 + q 1 + 2 λ 2 ,
p 3 = q ˙ 3 + q 1 + λ 2 ,
p 4 = q ˙ 4 λ 2 .
Using these results along with Equation (24a,b) for the Lagrange multipliers, we find that Equation (24e,f,g,h) agree precisely with Lagrange’s Equation (4).

6. Hamilton’s Equations and the Initial Value Problem

At this point one might ask whether the task of expressing the singular system (1) in Hamiltonian form is complete. After all, the primary action (23) provides the correct equations of motion for the phase space variables q i and p i . In fact, we can obtain the time evolution of any phase space function F from F ˙ = [ F , H P ] where H P is the primary Hamiltonian and [ · , · ] denotes Poisson brackets. Hamilton’s equations for the coordinates and momenta, q ˙ i = [ q i , H P ] and p ˙ i = [ p i , H P ] , coincide with Equation (24a–h).
Our task of expressing the singular system in Hamiltonian form is not yet complete because we still need to interpret Hamilton’s Equation (24a–h) as an initial value problem. That is, Hamilton’s equations should determine the future history of the system solely from initial data. In contrast, the primary action (23) defines a boundary value problem in which the configuration variables, the q’s, are specified at initial and final times.
The key difference between the equations of motion δ S P = 0 and Hamilton’s equations F ˙ = [ F , H P ] is that the former include the primary constraints, Equation (24i,j), whereas the latter do not. Thus, the phase space trajectories that extremize the action S P must lie entirely in the primary constraint surface. (The primary constraint “surface” is the subspace of phase space that satisfies the primary constraints.) In contrast, the trajectories obtained from Hamiltonian evolution F ˙ = [ F , H P ] are defined throughout the entire phase space. Note that we cannot simply append the primary constraint equations to Hamilton’s equations, because in that case the complete system would not be in Hamiltonian form.
Of course, the physically allowed phase space trajectories must satisfy the primary constraints. With an initial value interpretation of Hamilton’s equations, we can try to enforce the primary constraints with appropriate choices of initial data and Lagrange multipliers. In particular, we can choose initial data that lie on the primary constraint surface ϕ 1 ( 0 ) p 1 ( 0 ) = 0 and ϕ 2 ( 0 ) p 2 ( 0 ) p 3 ( 0 ) + p 4 ( 0 ) = 0 . However, this is not enough, because the primary constraints are not necessarily satisfied at later times as the system evolves into the future.
We can describe the situation as follows. The trajectories that extremize the action S P , the physical trajectories, do not necessarily fill the entire primary constraint surface. Instead they might span only a subspace of the primary constraint surface. If the initial data lie in the primary constraint surface but outside the subspace of physical trajectories, then the primary constraints will not be preserved as the data are evolved.
How should the initial data and Lagrange multipliers be restricted such that the primary constraints hold throughout the evolution? The primary constraints will hold for all time if they hold initially and their time derivatives (to all orders) also vanish initially. In the general case, this leads to a hierarchy of restrictions on the initial data in the form of secondary, tertiary, etc., constraints.2 It can also lead to restrictions on the Lagrange multipliers.
The higher order (secondary, tertiary, etc.) constraints and restrictions on the Lagrange multipliers are not new—imposing them does not change the content or predictions of the physical theory. This is because the higher-order constraints and restrictions on the Lagrange multipliers are direct consequences of the equations of motion (24) that follow from the primary action (23). They are simply “hidden” in those equations. The process of identifying the higher-order constraints and restrictions on Lagrange multipliers reveals these hidden conditions.

7. Consistency Conditions, Secondary Constraints and Restrictions on the Lagrange Multipliers

We can ensure that the primary constraints hold for all time by applying Dirac’s consistency conditions [10]. Begin by computing the time derivatives of the primary constraints with the primary Hamiltonian, ϕ ˙ a = [ ϕ a , H P ] . Now set these equal to zero:
[ ϕ a , H P ] = 0 .
For each value of the index a, there are three possibilities.3 First, [ ϕ a , H P ] might vanish on the constraint surface ϕ a = 0 , so that the consistency condition (26) reduces to the identity 0 = 0 . Second, [ ϕ a , H P ] could be a (non-constant) phase space function that is independent of the Lagrange multipliers. In this case, Equation (26) is a secondary constraint. Finally, [ ϕ a , H P ] might depend on the Lagrange multipliers. Then, Equation (26) fixes one of the Lagrange multipliers in terms of the phase space variables and the other Lagrange multipliers.
The secondary constraints that arise from this process must themselves satisfy the consistency conditions. This can lead to tertiary constraints and more restrictions on the Lagrange multipliers. In turn, the tertiary constraints can lead to quaternary constraints, and so forth. We must continue to apply the consistency conditions until the process naturally stops.
For our example, the primary constraints are
ϕ 1 p 1 = 0 ,
ϕ 2 p 2 p 3 + p 4 = 0 ,
and their time derivatives are
ϕ ˙ 1 = [ ϕ 1 , H P ] = p 3 + q 1 + q 2 + q 4 ,
ϕ ˙ 2 = [ ϕ 2 , H P ] = 2 ( q 1 + q 2 + q 4 ) .
Thus, we find the secondary constraints
ψ 1 p 3 + q 1 + q 2 + q 4 = 0 ,
ψ 2 2 ( q 1 + q 2 + q 4 ) = 0 .
These will be denoted collectively by ψ a .
Applying the consistency conditions to the secondary constraints gives
ψ ˙ 1 = [ ψ 1 , H P ] = p 4 + λ 1 + 2 λ 2 = 0 ,
ψ ˙ 2 = [ ψ 2 , H P ] = 2 ( p 4 + λ 1 + 2 λ 2 ) = 0 .
These equations restrict the Lagrange multipliers to satisfy
p 4 + λ 1 + 2 λ 2 = 0 .
The process has now terminated. In this example, there are no tertiary or higher-order constraints.
Recall from the previous section that our goal was to restrict the initial data and Lagrange multipliers such that the primary constraints vanish for all times under the Hamiltonian evolution defined by H P . We achieve this by imposing the primary constraints at the initial time,
ϕ 1 ( 0 ) p 1 ( 0 ) = 0 ,
ϕ 2 ( 0 ) p 2 ( 0 ) p 3 ( 0 ) + p 4 ( 0 ) = 0 ,
the secondary constraints at the initial time,
ψ 1 ( 0 ) p 3 ( 0 ) + q 1 ( 0 ) + q 2 ( 0 ) + q 4 ( 0 ) = 0 ,
ψ 2 ( 0 ) 2 [ q 1 ( 0 ) + q 2 ( 0 ) + q 4 ( 0 ) ] = 0 ,
and restricting the Lagrange multipliers to satisfy Equation (31) for all time t.
Let us review the reasoning. From Equation (30), the restriction (31) on the Lagrange multipliers tells us that ψ ˙ a = 0 for all time. By Equation (28), ψ a vanishes initially, so we see that ψ a must vanish for all time. Now we use Equation (32) to conclude that ϕ ˙ a must vanish for all time. Since ϕ a vanishes initially, by Equation (Section 7), it follows that the primary constraints ϕ a = 0 must hold for all time t.

8. Weak Equality and Lagrange Multiplier Analysis

It will be useful to follow the general Dirac–Bergmann algorithm closely and carry out a formal analysis of the restriction (31) on the Lagrange multipliers [10,11]. We begin with the concept of weak equality.
Let C A denote the complete set of (primary, secondary, tertiary, etc.) constraints. For our example,
C A ϕ 1 ϕ 2 ψ 1 ψ 2 = p 1 p 2 p 3 + p 4 p 3 + q 1 + q 2 + q 4 2 ( q 1 + q 2 + q 4 )
where the index A runs from 1 to 4.
Two phase space functions F and G are weakly equal if they are equal when the (primary, secondary, tertiary, etc.) constraints hold. In other words, F and G are weakly equal if they coincide on the constraint surface, the subspace of phase space defined by C A = 0 . Weak equality is written as F G .
Functions F and G are strongly equal if they agree throughout phase space. Strong equality is written as F = G .
Now we turn to the formal analysis of the restriction (31) on the Lagrange multipliers. This restriction can be expressed as the weak equality C ˙ A 0 . From Equations (28) and (30), we have
C ˙ A = [ C , H P ] = p 3 + q 1 + q 2 + q 4 2 ( q 1 + q 2 + q 4 ) p 4 + λ 1 + 2 λ 2 2 ( p 4 + λ 2 + 2 λ 2 ) 0 0 0 0
which simplifies to
0 0 p 4 2 p 4 + 0 0 0 0 1 2 2 4 λ 1 λ 2 0 0 0 0 .
This is a system of inhomogeneous linear equations for the Lagrange multipliers. A particular solution is
λ 1 λ 2 particular = 0 p 4 / 2 ,
and the homogeneous solutions are
λ 1 λ 2 homogeneous = 1 1 / 2 λ ,
where λ is arbitrary. The general solution is the sum of particular and homogeneous solutions:
λ 1 λ 2 general = λ ( λ + p 4 ) / 2 .
Thus, the restriction (31) on the Lagrange multipliers yields λ 1 = λ and λ 2 = ( λ + p 4 ) / 2 , where λ is an arbitrary function of time.

9. Total Hamiltonian

The total Hamiltonian  H T is obtained from the primary Hamiltonian H P by inserting the general solution for the Lagrange multipliers:
H T = H P | λ 1 = λ , λ 2 = ( λ + p 4 ) / 2 = H C + λ ϕ 1 ( λ + p 4 ) ϕ 2 / 2 = 1 2 p 3 2 p 2 p 4 + p 3 p 4 2 p 3 q 1 ( q 1 + 2 q 2 ) ( q 1 + 2 q 4 ) + λ ( p 1 p 2 / 2 + p 3 / 2 p 4 / 2 ) .
Physical phase space trajectories are defined by the total Hamiltonian as the weak equality F ˙ [ F , H T ] , with initial data that satisfy the complete set of constraints, C A = 0 .
Hamilton’s equations for the total Hamiltonian H T are
q ˙ 1 λ ,
q ˙ 2 p 4 / 2 λ / 2 ,
q ˙ 3 p 3 + p 4 / 2 q 1 + λ / 2 ,
q ˙ 4 p 3 / 2 p 2 / 2 λ / 2 ,
and
p ˙ 1 p 3 + q 1 + q 2 + q 4 ,
p ˙ 2 q 1 + 2 q 4 ,
p ˙ 3 0 ,
p ˙ 4 q 1 + 2 q 2 .
Since these are weak equalities, we can use the constraints to simplify the results. Observe that the constraints C A = 0 imply p 1 = p 3 = 0 , p 2 + p 4 = 0 and q 1 + q 2 + q 4 = 0 . Therefore, we can set p 1 and p 3 to zero, replace p 4 with p 2 , and replace q 4 with q 1 q 2 . Then Equation (41d,e,g,h) are either redundant or vacuous, and the remaining equations are
q ˙ 1 λ ,
q ˙ 2 p 2 / 2 λ / 2 ,
q ˙ 3 p 2 / 2 q 1 + λ / 2 ,
p ˙ 2 q 1 2 q 2 .
These equations, along with the constraints C A = 0 , give a complete description of the physical system.
Let us check the results. The Lagrange multiplier λ can be eliminated from Equation (42a,b) to give q ˙ 1 + 2 q ˙ 2 p 2 . Now differentiate this equation and eliminate p ˙ 2 with Equation (42d) to obtain q ¨ 1 + 2 q ¨ 2 q 1 2 q 2 . The constraints allow us to set q 1 q 2 q 4 , which gives
q ¨ 4 q ¨ 2 q 2 q 4 .
This is Equation (5c), which follows directly from Lagrange’s equations. The result (5a) from Lagrange’s equations is simply the secondary constraint C 4 2 ( q 1 + q 2 + q 4 ) = 0 . Finally, the result (5b) is obtained by summing Equation (42b,c).
Recall that Equation (5a–c) are equivalent to Lagrange’s equations. Thus, we have verified that Hamilton’s Equation (41), along with the primary and secondary constraints C A = 0 , are equivalent to Lagrange’s equations.

10. First and Second Class Constraints

A first class function F is a phase space function that has weakly vanishing Poisson brackets with all primary and secondary constraints:
[ F , C A ] 0 F is first class .
It can be shown that the Poisson bracket of any two first class functions is itself a first class function [10,11].
The constraints themselves can be first class; constraints that are not first class are called second class. The constraints are separated into first and second class by examining the matrix of Poisson brackets:
[ C A , C B ] = 0 0 1 2 0 0 2 4 1 2 0 0 2 4 0 0 .
The rank of this 4 × 4 matrix is 2, and its nullity is 4 2 = 2 . It follows that there are 2 independent eigenvectors with eigenvalues equal to zero; for example, u A = ( 1 , 1 / 2 , 0 , 0 ) and v A = ( 0 , 0 , 1 , 1 / 2 ) . Then there are two independent combinations of constraints that are first class, namely u A C A and v A C A . (A sum over the repeated index A is implied.) The first class constraints are
C 1 ( f c ) ϕ 1 ϕ 2 / 2 = p 1 p 2 / 2 + p 3 / 2 p 4 / 2 ,
C 2 ( f c ) ψ 1 ψ 2 / 2 = p 3 .
One can check that the first class conditions [ C 1 ( f c ) , C B ] = 0 and [ C 2 ( f c ) , C B ] = 0 hold. The most general first class constraint is a linear combination of C 1 ( f c ) and C 2 ( f c ) .
There are two remaining linear combinations of constraints, which we take to be
C 1 ( s c ) ( ϕ 1 + ϕ 2 ) / 3 = ( p 1 + p 2 p 3 + p 4 ) / 3 ,
C 2 ( s c ) ( ψ 1 + ψ 2 ) / 3 = p 3 / 3 + q 1 + q 2 + q 4 .
These are the second class constraints. They have nonvanishing Poisson brackets with each other,
[ C 2 ( s c ) , C 1 ( s c ) ] = 1 .
The most general second class constraint is a linear combination of C 1 ( s c ) , C 2 ( s c ) , C 1 ( f c ) and C 2 ( f c ) , with nonzero coefficients on one or both of C 1 ( s c ) and C 2 ( s c ) .
The splitting of constraints into first and second class is independent of the splitting into primary and secondary. In this example, the first class constraints are mixtures of primary and secondary constraints. Likewise, the second class constraints are mixtures of primary and secondary constraints.

11. First Class Hamiltonian, Gauge Generators and the Dirac Conjecture

The total Hamiltonian (40) includes the product of an arbitrary Lagrange multiplier λ with the first class constraint C 1 ( f c ) ϕ 1 ϕ 2 / 2 . We refer to C 1 ( f c ) as a primary first class constraint, since it is constructed entirely from primary constraints.
If we remove the primary first class constraint from the total Hamiltonian, what remains is the first class Hamiltonian  H f c . That is, the total Hamiltonian can be written as
H T = H f c + λ C 1 ( f c ) ,
where
H f c = 1 2 p 3 2 p 2 p 4 + p 3 p 4 2 p 3 q 1 ( q 1 + 2 q 2 ) ( q 1 + 2 q 4 )
is the first class Hamiltonian. A common notation for H f c , the notation used by Dirac [10], is H .
We can check directly that the first class Hamiltonian (50) is a first class function. However, this is not necessary, because we know that the constraints are preserved under the time evolution defined by H T . That is, C ˙ A = [ C A , H T ] 0 . Thus, the total Hamiltonian must be first class, [ H T , C A ] 0 . Of course, the primary first class constraint C 1 ( f c ) is first class. It then follows from the definition (49) that H f c must also be a first class function.
The splitting (49) of the total Hamiltonian into the first class Hamiltonian and the primary first class constraint is not special to our example. This splitting will occur for any constrained Hamiltonian system [10,11]. In general, H T will include the products of every primary first class constraint with an arbitrary multiplier.
Primary first class constraints generate gauge transformations. Consider the change in a phase space function F generated by the primary first class constraint C 1 ( f c ) ,
δ F = δ ϵ [ F , C 1 ( f c ) ] .
This transformation does not change the physical state of the system. We can see this by considering F to be evaluated as a function of the q’s and p’s at some particular time t. At an infinitesimally later time t + δ t , this function becomes F ( t + δ t ) = F ( t ) + [ F , H T ] δ t . In terms of the first class Hamiltonian, we have
F ( t + δ t ) = F ( t ) + [ F , H f c ] + λ [ F , C 1 ( f c ) ] δ t .
The Lagrange multiplier is arbitrary, so we can make a different choice during the time interval from t to t + δ t , say, λ ˜ . Then, the function F at time t + δ t will be
F ˜ ( t + δ t ) = F ( t ) + [ F , H f c ] + λ ˜ [ F , C 1 ( f c ) ] δ t .
The physical state of the system at t + δ t should not depend on our choice of Lagrange multiplier, so F ( t + δ t ) and F ˜ ( t + δ t ) must represent the same physical state. The result (51) is obtained by subtracting Equation (52) from Equation (53) and defining δ F F ˜ F and δ ϵ ( λ ˜ λ ) δ t .
For the phase space coordinates, the gauge transformation generated by the primary first class constraint C 1 ( f c ) is
δ q 1 = δ ϵ ,
δ q 2 = δ ϵ / 2 ,
δ q 3 = δ ϵ / 2 ,
δ q 4 = δ ϵ / 2 .
The transformations of the p’s all vanish. This result agrees with the gauge transformation from Equation (14), with the change of notation ϵ δ ϵ . Here, we denote the gauge parameter by δ ϵ because the transformation is infinitesimal; in Equation (14), we used ϵ because the transformation was finite. It is clear that the infinitesimal transformation (54) can be iterated to obtain the finite transformation (14).
In general, a gauge transformation is defined as a transformation δ F = δ ϵ [ F , G ] that does not alter the physical state of the system. The function G is the gauge generator. We have seen that the primary first class constraints generate gauge transformations. However, not all gauge transformations are generated by primary first class constraints. In fact, it can be shown [10,11] that the Poisson bracket between any primary first class constraint and the first class Hamiltonian is itself a first class constraint that generates a gauge transformation.4
For our example problem, the Poisson bracket of the primary first class constraint C 1 ( f c ) and the first class Hamiltonian H f c is
[ C 1 ( f c ) , H f c ] = p 3 .
This is the secondary first class constraint, C 2 ( f c ) ψ 1 ψ 2 / 2 = p 3 . Thus, we see that in this example both the primary and secondary first class constraints are generators of gauge transformations. Explicitly, the transformation δ F = δ ϵ [ F , C 2 ( f c ) ] is
δ q 1 = 0 ,
δ q 2 = 0 ,
δ q 3 = δ ϵ ,
δ q 4 = 0 ,
with the transformations of the p’s all vanishing. We can iterate this infinitesimal gauge transformation (56) to obtain the finite transformation (16).
The “doubling” of the gauge freedom identified in Section 3 appears quite naturally in the Hamiltonian formalism. The two types of gauge transformation are generated by the two first class constraints, C 1 ( f c ) and C 2 ( f c ) .
The Dirac conjecture [10] says that all first class constraints (whether they are primary, secondary, etc., or a combination of primary, secondary, etc.) generate gauge transformations. This conjecture does not hold as a general theorem—there are known examples in which the transformation generated by a secondary first class constraint does not coincide with any invariance of the original Lagrangian system.5 Nevertheless, the Dirac conjecture is usually taken as an assumption. It appears that in practice, for systems of physical interest, all first class constraints generate gauge transformations.

12. Extended Hamiltonian

The Dirac conjecture tells us that all first class constraints generate gauge transformations and should be treated on an equal footing. The extended Hamiltonian  H E is defined by adding all first class constraints C a ( f c ) with Lagrange multipliers γ a to the first class Hamiltonian:
H E = H f c + γ a C a ( f c ) .
(A sum over the index a is implied.)
For our example, the extended Hamiltonian is
H E = 1 2 p 3 2 p 2 p 4 + p 3 p 4 2 p 3 q 1 ( q 1 + 2 q 2 ) ( q 1 + 2 q 4 ) + γ 1 ( p 1 p 2 / 2 + p 3 / 2 p 4 / 2 ) + γ 2 p 3 .
The equations of motion F ˙ [ F , H E ] are
q ˙ 1 γ 1 ,
q ˙ 2 p 4 / 2 γ 1 / 2 ,
q ˙ 3 p 3 + p 4 / 2 q 1 + γ 1 / 2 + γ 2 ,
q ˙ 4 p 3 / 2 p 2 / 2 γ 1 / 2 ,
and
p ˙ 1 p 3 + q 1 + q 2 + q 4 ,
p ˙ 2 q 1 + 2 q 4 ,
p ˙ 3 0 ,
p ˙ 4 q 1 + 2 q 2 .
Let us compare these results to the equations of motion (41) obtained from the total Hamiltonian H T . There are just two differences. The first is trivial: the Lagrange multiplier λ in Equation (41) has changed names to γ 1 in Equation (59). The second difference is significant: the equation for q ˙ 3 has an extra term γ 2 on the right-hand side. This is a new feature of the extended Hamiltonian. It makes explicit the fact that the gauge freedom allows q 3 to be changed arbitrarily, and independently, from the other variables.
We can check the equations of motion for H E following the same reasoning that was applied to the equations of motion for H T . First, recall that the (first and second class) constraints imply p 1 = p 3 = 0 , p 4 = p 2 and q 4 = q 1 q 2 . Then Equation (59e,g) are vacuous, and Equation (59f,h) are redundant. It also follows that with the constraints imposed, Equation (59d) is a consequence of Equation (59a,b). The remaining equations are
q ˙ 1 γ 1 ,
q ˙ 2 p 2 / 2 γ 1 / 2 ,
q ˙ 3 p 2 / 2 q 1 + γ 1 / 2 + γ 2 ,
p ˙ 2 q 1 2 q 2 .
These agree with Equation (42), apart from the change of notation λ γ 1 and the extra term γ 2 on the right-hand side of the q ˙ 3 equation.
By eliminating γ 1 , the equations of motion generated by the extended Hamiltonian H E become
q ˙ 1 + 2 q ˙ 2 p 2 ,
q ˙ 2 + q ˙ 3 q 1 + γ 2 ,
p ˙ 2 q 1 2 q 2 .
If we differentiate the first equation, combine with the second, and use the constraint q 1 + q 2 + q 4 = 0 , we obtain ( q ¨ 4 q ¨ 2 ) ( q 4 q 2 ) . This is the expected result (5c). In fact, the only difference between Hamilton’s equations F ˙ [ F , H E ] and the results (5) (which are equivalent to Lagrange’s equations) is the extra term γ 2 in Equation (61b) above. That term does not appear in the corresponding Lagrangian Equation (5b).
Note that we can use the first class constraints to simplify the extended Hamiltonian H E . For example, using C 2 ( f c ) = p 3 , we can set p 3 = 0 everywhere in Equation (58), except of course in the term γ 2 p 3 . The extended Hamiltonian becomes
H E = 1 2 p 2 p 4 ( q 1 + 2 q 2 ) ( q 1 + 2 q 4 ) + γ 1 ( p 1 p 2 / 2 p 4 / 2 ) + γ 2 p 3 .
This amounts to replacing the Lagrange multiplier γ 2 in Equation (58) by
γ 2 γ 2 p 3 / 2 p 4 / 2 + q 1 γ 1 / 2 .
This replacement does not change the physical content of the theory, since the Lagrange multiplier γ 2 is arbitrary.

13. Extended Action

The equations of motion for the extended theory can be derived from the action [11]
S E [ q , p , γ , σ ] = 0 T d t p i q ˙ i H E σ 1 C 1 ( s c ) σ 2 C 2 ( s c ) ,
which includes the second class constraints with Lagrange multipliers σ a . Recall that the first class constraints C a ( f c ) are included in the extended Hamiltonian H E with multipliers γ a , so S E includes all four constraints.
We can use either form of the extended Hamiltonian, Equation (58) or (62), in the extended action. Let us use Equation (58). Then, the equations of motion that follow from extremizing S E with respect to the momenta p i are
q ˙ 1 = γ 1 + σ 1 / 3 ,
q ˙ 2 = p 4 / 2 γ 1 / 2 + σ 1 / 3 ,
q ˙ 3 = p 3 + p 4 / 2 q 1 + γ 1 / 2 + γ 2 σ 1 / 3 + σ 2 / 3 ,
q ˙ 4 = p 3 / 2 p 2 / 2 γ 1 / 2 + σ 1 / 3 .
Extremizing S E with respect to the coordinates yields
p ˙ 1 = p 3 + q 1 + q 2 + q 4 σ 2 ,
p ˙ 2 = q 1 + 2 q 4 σ 2 ,
p ˙ 3 = 0 ,
p ˙ 4 = q 1 + 2 q 2 σ 2 ,
and the constraints
C 1 ( f c ) p 1 p 2 / 2 + p 3 / 2 p 4 / 2 = 0 ,
C 2 ( f c ) p 3 = 0 ,
C 1 ( s c ) ( p 1 + p 2 p 3 + p 4 ) / 3 = 0 ,
C 2 ( s c ) p 3 / 3 + q 1 + q 2 + q 4 = 0 ,
follow from extremizing S E with respect to the Lagrange multipliers γ a and σ a .
Let us check these equations of motion. The constraints imply p 1 = p 3 = 0 , p 4 = p 2 and q 1 + q 2 + q 4 = 0 . Then the equation of motion (65e) gives σ 2 0 . The constraints also imply q ˙ 1 + q ˙ 2 + q ˙ 4 = 0 . The sum of Equation (65a,b,d) then yields σ 1 0 . Now, if we set σ 1 and σ 2 to zero, the equations of motion (65a–h) agree precisely with Hamilton’s Equation (59) for the extended Hamiltonian H E . In Section 12, we showed that the equations generated by H E agree with Lagrange’s equations apart from the extra term γ 2 in the equation for q ˙ 3 . This term extends the original Lagrangian theory by making explicit the fact that the gauge freedom allows for independent transformations of q 3 .
Finally, we note that the extended action is invariant under the transformation defined by
δ F = ϵ 1 [ F , C 1 ( f c ) ] + ϵ 2 [ F , C 2 ( f c ) ]
for the phase space variables and
δ γ 1 = ϵ ˙ 1 ,
δ γ 2 = ϵ ˙ 2 + ϵ 1 ,
δ σ 1 = 0 ,
δ σ 2 = 0 ,
for the Lagrange multipliers. Here, the gauge parameters ϵ 1 and ϵ 2 are functions of time. These equations express the gauge invariance at the level of the action S E .

14. Dirac Brackets and the Partially Reduced Hamiltonian

We now return to the evolution defined by the extended Hamiltonian H E of Equation (58), and Hamilton’s Equation (59). To obtain a physically allowed trajectory, we must choose initial data that satisfy the four constraints C a ( f c ) = 0 and C a ( s c ) = 0 . Apart from restricting the initial data, the second class constraints play no role in the formalism. It would be convenient if we could restrict the variables from the outset such that the second class constraints are automatically satisfied. For example, we could use C 1 ( s c ) = 0 and C 2 ( s c ) = 0 from Equation (65k,l) to replace q 1 with q 2 q 4 p 3 / 3 and replace p 2 with p 1 + p 3 p 4 .
We are not allowed to apply the second class constraints in this way. For example, consider the Poisson brackets [ q 1 , p 1 ] = 1 . If we were to replace q 1 with q 2 q 4 p 3 / 3 , we would find a different answer: [ q 2 q 4 p 3 / 3 , p 1 ] = 0 . The second class constraints cannot be imposed before Poisson brackets are computed.
Dirac devised a way to allow the second class constraints to be imposed from the outset by modifying the Poisson brackets [10]. The result is the Dirac brackets.
To construct Dirac brackets, we first compute the matrix of Poisson brackets among the second class constraints:
M a b [ C a ( s c ) , C b ( s c ) ] = 0 1 1 0 .
Let
M a b = 0 1 1 0
denote the inverse of M a b . Then, the Dirac brackets [ F , G ] * of two phase space functions F and G are defined by
[ F , G ] * [ F , G ] [ F , C a ( s c ) ] M a b [ C b ( s c ) , G ] .
Dirac brackets, like Poisson brackets, are antisymmetric and satisfy the Jacobi identity [10,11].
Explicitly, the Dirac brackets among the coordinates are
[ q i , q j ] * = 1 9 0 0 1 0 0 0 1 0 1 1 0 1 0 0 1 0 ,
and the Dirac brackets between the q’s and p’s are
[ q i , p j ] * = 1 3 2 1 0 1 1 2 0 1 1 1 3 1 1 1 0 2
For our example, the Dirac brackets among the momenta all vanish: [ p i , p j ] * = 0 .
There are two key properties that make Dirac brackets relevant. First, the Dirac brackets agree weakly with Poisson brackets if one of the two functions is first class. Since the extended Hamiltonian is first class, we have [ F , H E ] * [ F , H E ] for any F. It follows that we can write the equations of motion as
F ˙ [ F , H E ] * ,
using Dirac brackets.
The second key property of the Dirac brackets is that they weakly vanish if one of the functions is a second class constraint: [ F , C a ( s c ) ] * 0 . This allows us to apply the second class constraints before computing brackets. For example, we can use either q 1 or q 2 q 4 p 3 / 3 to compute Dirac brackets with p 1 :
[ q 1 , p 1 ] * = [ q 2 q 4 p 3 / 3 , p 1 ] * = 2 / 3 .
With Dirac brackets, the second class constraints can be treated as strong equations and imposed before computing the equations of motion.
Let us use the second class constraints (65k,l) to eliminate q 1 and p 2 and write the extended Hamiltonian in terms of the smaller set of variables q 2 , q 3 , q 4 , p 1 , p 3 and p 4 . Setting
q 1 = q 2 q 4 p 3 / 3 ,
p 2 = p 1 + p 3 p 4 ,
we have
H P R = 7 p 3 2 / 9 + 1 2 p 4 2 + p 1 p 4 + 2 p 3 ( q 2 + q 4 ) + ( q 2 q 4 ) 2 + γ 1 ( 3 p 1 / 2 ) + γ 2 p 3 .
This is the partially reduced Hamiltonian, obtained from the extended Hamiltonian by applying the second class constraints.
Of course, the partially reduced Hamiltonian is not unique. We could use the second class constraints to eliminate some other pair of variables instead of q 1 and p 2 .
The equations of motion generated by the partially reduced Hamiltonian, F ˙ [ F , H R ] * , are
q ˙ 2 p 1 / 6 p 4 / 2 γ 1 / 2 ,
q ˙ 3 4 p 3 / 3 + p 1 / 6 + p 4 / 2 + q 2 + q 4 + γ 1 / 2 + γ 2 ,
q ˙ 4 p 1 / 3 + p 4 / 2 γ 1 / 2 ,
p ˙ 1 2 p 3 / 3 ,
p ˙ 3 0 ,
p ˙ 4 p 3 / 3 + q 2 q 4 .
We can also use H P R and the Dirac brackets to compute q ˙ 1 and p ˙ 2 . The results are equivalent to those obtained by differentiating the right-hand sides of Equation (74) and using the equations of motion (76).
Let us check the equations of motion. With the second class constraints applied, the first class constraints imply p 1 = p 3 = 0 and p 4 = p 2 . Thus, Equation (76d,e) are vacuous and the remaining equations become
q ˙ 2 p 2 / 2 γ 1 / 2 ,
q ˙ 3 p 2 / 2 + q 2 + q 4 + γ 1 / 2 + γ 2 ,
q ˙ 4 p 2 / 2 γ 1 / 2 ,
p ˙ 2 q 2 + q 4 ,
Compare these to the independent Equation (60) that follow from the extended Hamiltonian. Equation (60b–d) agree with Equation (77a,b,d) once we use q 1 = q 2 q 4 . The final Equation (60a) is obtained by differentiating q 1 = q 2 q 3 in time and using Equation (77a,b).

15. Partially Reduced Action

The partially reduced equations of motion (76) can be obtained from the extended action S E by eliminating the superfluous variables. Note that the equations of motion obtained by varying S E with respect to p 2 , q 1 , σ 1 and σ 2 are Equation (65b,e,k,l), respectively. We can eliminate these variables by solving these equations and substituting the results into the action.6 The results are:
q 1 = p 3 / 3 q 2 q 4 ,
p 2 = p 1 + p 3 p 4 ,
σ 1 = 3 q ˙ 2 + 3 p 4 / 2 + 3 γ 1 / 2 ,
σ 2 = p ˙ 1 + 2 p 3 / 3 .
Inserting these into the extended action (64), we find
S P R [ q 2 , q 3 , q 4 , p 1 , p 3 , p 4 , γ 1 , γ 2 ] = 0 T d t { p 1 p ˙ 3 / 3 + ( p 3 p 4 2 p 1 ) q ˙ 2 + p 3 q ˙ 3 + ( p 4 p 1 ) q ˙ 4 H P R } .
This is the partially reduced action.
The equations of motion obtained from varying S P R with respect to the phase space variables are
2 p ˙ 1 p ˙ 3 + p ˙ 4 p 3 q 2 + q 4 = 0 ,
p ˙ 3 = 0 ,
p ˙ 1 p ˙ 4 p 3 + q 2 q 4 = 0 ,
2 q ˙ 2 q ˙ 4 p ˙ 3 / 3 p 4 / 2 3 γ 1 / 2 = 0 ,
q ˙ 2 + q ˙ 3 + p ˙ 1 / 3 14 p 3 / 9 q 2 q 4 γ 2 = 0 ,
q ˙ 2 + q ˙ 4 p 1 / 2 p 4 = 0 ,
and the equations obtained by varying with respect to the Lagrange multipliers γ 1 and γ 2 are
3 p 1 / 2 = 0 ,
p 3 = 0 .
These are, of course, the first class constraints, reduced by using the second class constraints to eliminate q 1 and p 2 .
We can now solve Equation (80) for the time derivatives of q 2 , q 3 , q 4 , p 1 , p 3 and p 4 . The result coincides with the equations of motion (76) obtained from the partially reduced Hamiltonian H P R and the Dirac brackets.
The partially reduced action S P R is invariant under the transformation defined by
δ F = ϵ 1 [ F , C 1 ( f c ) ] * + ϵ 2 [ F , C 2 ( f c ) ] *
for the phase space variables and
δ γ 1 = ϵ ˙ 1 ,
δ γ 2 = ϵ ˙ 2 + ϵ 1 ,
for the Lagrange multipliers. These equations express the gauge invariance of the theory at the level of the action principle with the second class constraints eliminated.

16. Partially Reduced Lagrangian

It is not too difficult to find a change of variables that will bring S P R into “canonical form”. For example, let
q 2 = Q 1 P 2 / 9 ,
q 3 = Q 2 + P 2 / 9 ,
q 4 = Q 3 P 2 / 9 ,
p 1 = ( P 1 + P 2 P 3 ) / 3 ,
p 3 = P 2 ,
p 4 = ( P 1 + P 2 + 2 P 3 ) / 3 ,
define a new set of variables Q α , P α for the secondary constraint surface. (The index α ranges over 1, 2 and 3.) The partially reduced action becomes
S P R [ Q , P , γ ] = 0 T d t P α Q ˙ α H P R
with
H P R = 1 18 2 P 1 2 + 12 P 2 2 + 2 P 3 2 4 P 1 P 2 + 5 ( P 2 P 1 ) P 3 + 18 P 2 ( Q 1 + Q 3 ) + 9 ( Q 1 Q 3 ) 2 + γ 1 ( P 1 + P 2 P 3 ) / 2 + γ 2 P 2 .
The equations of motion δ S P R = 0 include Q ˙ α = [ Q α , H P R ] and P ˙ α = [ P α , H P R ] , where [ · , · ] are the usual Poisson brackets. The partially reduced action (84) is invariant under the transformation
δ F = ϵ 1 [ F , C 1 ( f c ) ] + ϵ 2 [ F , C 2 ( f c ) ] ,
δ γ 1 = ϵ ˙ 1 ,
δ γ 2 = ϵ ˙ 2 + ϵ 1 ,
where C 1 ( f c ) = ( P 1 + P 2 P 3 ) / 2 and C 2 ( f c ) = P 2 are the first class constraints.
We can use the equations of motion δ S P R / δ P α = 0 to eliminate the momenta P α from the partially reduced action (84). These equations are
Q ˙ 1 = ( 4 P 1 4 P 2 5 P 3 ) / 18 γ 1 / 2 ,
Q ˙ 2 = ( 4 P 1 24 P 2 5 P 3 ) / 18 + Q 1 + Q 3 + γ 1 / 2 + γ 2 ,
Q ˙ 3 = ( 5 P 1 5 P 2 4 P 3 ) / 18 γ 1 / 2 ,
with solutions
P 1 = ( 71 Q ˙ 1 + 9 Q ˙ 2 100 Q ˙ 3 ) / 10 9 ( Q 1 + Q 3 ) / 10 9 γ 1 9 γ 2 / 10 ,
P 2 = 9 ( Q ˙ 1 + Q ˙ 2 Q 1 Q 3 γ 2 ) / 10 ,
P 3 = 10 Q ˙ 1 8 Q ˙ 3 9 γ 1 .
Inserting these results into the partially reduced action, we find
S P R [ Q , γ ] = 0 T d t { ( 71 Q ˙ 1 2 + 9 Q ˙ 2 2 80 Q ˙ 3 2 ) / 20 + ( 9 Q ˙ 2 100 Q ˙ 3 ) Q ˙ 1 / 10 9 ( Q 1 + Q 3 ) ( Q ˙ 1 + Q ˙ 2 ) / 10 ( Q 1 2 + Q 3 2 38 Q 1 Q 3 ) / 20 9 ( Q ˙ 1 + Q ˙ 3 ) γ 1 9 ( Q ˙ 1 + Q ˙ 2 Q 1 Q 3 ) γ 2 / 10 9 ( 10 γ 1 2 γ 2 2 ) / 20 } .
The partially reduced Lagrangian is the integrand of this action.
We can go one step further and eliminate the Lagrange multipliers using the equations of motion δ S / δ γ 1 = 0 and δ S / δ γ 2 = 0 . These equations have solutions
γ 1 = Q ˙ 1 Q ˙ 3 ,
γ 2 = Q ˙ 1 + Q ˙ 2 Q 1 Q 3 .
Inserting these results back into the action yields
S P R [ Q ] = 0 T d t 1 2 ( Q ˙ 1 Q ˙ 3 ) 2 1 2 ( Q 1 Q 3 ) 2 .
This is the action for a harmonic oscillator in the variable Q 1 Q 3 . Note that the coordinate transformation (83) implies ( Q 1 Q 3 ) = ( q 2 q 4 ) , so once again, we find that the coordinate combination q 2 q 4 describes a simple harmonic oscillator. In addition, we observe that the action (91) leaves the variable Q 2 completely unspecified. This expresses the gauge freedom generated by the first class constraint C 2 ( f c ) = P 2 = p 3 .

17. Gauge Conditions and the Fully Reduced Hamiltonian

Let us return to the theory described by the extended Hamiltonian, prior to the elimination of the second class constraints.
For our example problem, phase space is eight-dimensional. The physical trajectories fill the constraint “surface”, which is the four-dimensional subspace where all first and second class constraints hold. Each point in the constraint surface can be mapped into a physically equivalent state by the gauge generators, namely, the first class constraints C a ( f c ) . Since there are two independent gauge generators, each physical state of the system corresponds to a two-dimensional subspace of the constraint surface. The constraint surface is foliated by these two-dimensional slices, referred to as gauge “orbits”.
We can select a single phase space point on each gauge orbit to represent the physical state. We do this by applying gauge conditions. In particular, we will consider a canonical gauge7 which takes the form G a ( q , p ) 0 with a = 1 , 2 . A good canonical gauge condition must not be gauge invariant, otherwise it would allow more than one point on the gauge orbit to represent the physical state of the system. To be precise, the matrix of Poisson brackets of gauge conditions and gauge generators, [ G a , C b ( f c ) ] , must be nonsingular [11].
As an example, let us choose
G 1 = q 1 q 2 ,
G 2 = q 3 + p 4
as our gauge conditions. This is a good gauge:
det [ G a , C b ( f c ) ] = 3 / 2 0 1 / 2 1 = 3 / 2 .
The matrix [ G a , C b ( f c ) ] is nonsingular, as required.
The gauge conditions G a = 0 , like the first and second class constraints, restrict the phase space variables. The full set of restrictions
C A ( a l l ) = { G 1 , G 2 , C 1 ( f c ) , C 2 ( f c ) , C 1 ( s c ) , C 2 ( s c ) }
reduces the available phase space from eight dimensions to two dimensions. (Here, the index A ranges from 1 to 6.) Taken as a whole, the six conditions C A ( a l l ) are second class. We see this by computing the Poisson brackets
M A B [ C A ( a l l ) , C B ( a l l ) ] = 1 6 0 0 9 0 0 0 0 0 3 6 2 4 9 3 0 0 0 0 0 6 0 0 0 0 0 2 0 0 0 6 0 4 0 0 6 0 .
This matrix has a nonzero determinant, det ( M ) = 9 / 4 , which is the condition for the set of constraints and gauge conditions to be second class.
We can eliminate the constraints and gauge conditions by constructing Dirac brackets. The inverse of M A B is
M A B = 1 3 0 0 2 1 0 0 0 0 0 3 0 0 2 0 0 0 0 0 1 3 0 0 2 1 0 0 0 2 0 3 0 0 0 1 3 0 ,
and the Dirac brackets are defined by
[ F , G ] * = [ F , G ] [ F , C A ( a l l ) ] M A B [ C B ( a l l ) , G ] .
The Dirac brackets among the phase space variables are
q 1 , q 3 * = [ q 1 , p 2 ] * = [ q 1 , p 4 ] * = 1 / 3 ,
q 2 , q 3 * = [ q 2 , p 2 ] * = [ q 2 , p 4 ] * = 1 / 3 ,
q 3 , q 4 * = [ q 4 , p 2 ] * = [ q 4 , p 4 ] * = 2 / 3 ,
with all other brackets vanishing.
The constraints can be solved in various ways and the results can be used freely, either before or after computing Dirac brackets. For example, the constraints imply
q 1 = q 4 / 2 ,
q 2 = q 4 / 2 ,
q 3 = p 2 ,
p 1 = 0 ,
p 3 = 0 ,
p 4 = p 2 .
We can use these to eliminate the variables q 1 , q 2 , q 3 , p 1 , p 3 and p 4 . Then the extended Hamiltonian H E becomes the fully reduced Hamiltonian
H F R = 1 2 p 2 2 + 9 4 q 4 2 ,
which depends only on q 4 and p 2 .
The Dirac brackets of the variables that remain are [ q 4 , p 2 ] * = 2 / 3 . Thus, the equations of motion become
q ˙ 4 = [ q 4 , H F R ] * = 2 p 2 / 3 ,
p ˙ 2 = [ p 4 , H F R ] * = 3 q 4 / 2 .
These are the equations for a simple harmonic oscillator with solution
q 4 ( t ) = α sin t + β cos t ,
p 2 ( t ) = ( 3 α / 2 ) cos t + ( 3 β / 2 ) sin t ,
where α and β are arbitrary constants.
With the gauge fixed, the dynamics take place on the fully reduced phase space, the two-dimensional surface defined by the constraints and gauge conditions C A ( a l l ) = 0 . There are many different choices of coordinates for this surface. Instead of solving the constraints for q 4 and p 2 , we could solve them for q 2 and q 3 . In that case, the fully reduced Hamiltonian is
H F R = 9 2 q 2 2 + 1 2 q 3 2
and Hamilton’s equations are
q ˙ 2 = [ q 2 , H F R ] * = q 3 / 3 ,
q ˙ 3 = [ q 3 , H F R ] * = 3 q 2
Again, this describes the simple harmonic oscillator.
We can choose other coordinates on the fully reduced phase space. For example, let q 4 = Q + q 2 and p 2 = P , then use the constraints to eliminate q 1 , q 2 , q 3 , p 1 , p 3 and p 4 . The fully reduced Hamiltonian becomes
H F R = 1 2 ( Q 2 + P 2 ) .
The nonzero Dirac brackets are [ Q , P ] * = 1 , and the equations of motion are simply Q ˙ = P and P ˙ = Q .
In each case, the fully reduced theory exhibits the single physical degree of freedom that we expect.

18. Fully Reduced Action

The fully reduced equations of motion can be derived from the action that includes all of the constraints and gauge conditions. For lack of a better name, let us denote this action with the subscript “all”:
S a l l = 0 T d t p i q ˙ i H f c γ a C a ( f c ) σ a C a ( s c ) ρ a G a .
Now extremize S a l l with respect to variations in p 1 , p 3 , p 4 , q 1 , q 2 , q 3 and q 4 :
q ˙ 1 = γ 1 + σ 1 / 3 ,
q ˙ 3 = q 1 + p 3 + p 4 / 2 + γ 1 / 2 + γ 2 σ 1 / 3 + σ 2 / 3 ,
q ˙ 4 = p 2 / 2 + p 3 / 2 γ 1 / 2 + σ 1 / 3 + ρ 2 ,
p ˙ 1 = q 1 + q 2 + q 4 + p 3 σ 2 ρ 1 ,
p ˙ 2 = q 1 + 2 q 4 + σ 2 ρ 1 ,
p ˙ 3 = ρ 2 .
In addition, vary S a l l with respect to the Lagrange multipliers to obtain the constraints C A ( a l l ) = 0 . The solution of the full set of equations, (107) and C A ( a l l ) = 0 , is given by Equation (99) along with
γ 1 = p 2 / 3 2 p ˙ 3 / 2 + 2 q ˙ 1 / 3 2 q ˙ 4 / 3 ,
γ 2 = p 2 3 q 4 / 4 + p ˙ 1 / 6 + p ˙ 2 / 6 + p ˙ 3 + q ˙ 3 + q ˙ 4 ,
σ 1 = p 2 + 2 p ˙ 3 + q ˙ 1 + 2 q ˙ 4 ,
σ 2 = 3 q 4 / 4 p ˙ 1 / 2 p ˙ 2 / 2 ,
ρ 1 = 3 q 4 / 4 p ˙ 1 / 2 + p ˙ 2 / 2 ,
ρ 2 = p ˙ 3 .
We now insert these results into S a l l to obtain the fully reduced action
S F R [ q 4 , p 2 ] = 0 T d t 3 2 p 2 q ˙ 4 1 2 p 2 2 + 9 4 q 4 2 ,
which is a functional of q 4 and p 2 . The equations of motion δ S F R = 0 are
0 = δ S F R δ q 4 = 3 2 p ˙ 2 9 4 q 4 ,
0 = δ S F R δ p 2 = 3 2 q ˙ 4 p 2 .
These are equivalent to Hamilton’s Equation (101) for the fully reduced Hamiltonian.
We can place S F R into “canonical form” by defining new variables P = p 2 and Q = 3 q 4 / 2 . Then
S F R [ Q , P ] = 0 T d t P Q ˙ 1 2 [ P 2 + Q 2 ] ,
which is the familiar action for the harmonic oscillator.

19. Summary and Discussion

Here is the Dirac–Bergmann algorithm:
  • Compute the conjugate momenta p i = L / q ˙ i and define the canonical Hamiltonian H C as p i q ˙ i L ( q , q ˙ ) , written in terms of p’s and q’s.
  • Identify the primary constraints. The primary Hamiltonian H P is obtained from H C by adding the primary constraints with Lagrange multipliers.
  • Apply Dirac’s consistency conditions to identify higher-order constraints and restrictions on the Lagrange multipliers.
  • The total Hamiltonian H T is found from H P by incorporating the restrictions on the Lagrange multipliers.
  • Separate the primary, secondary, and higher-order constraints into first and second class.
  • The first class Hamiltonian H f c is the part of H T with the primary first class constraints removed.
  • The extended Hamiltonian H E is obtained from the first class Hamiltonian H f c by adding all of the first class constraints with Lagrange multipliers.
  • The partially reduced Hamiltonian H P R is found from H E by using Dirac brackets to eliminate the second class constraints.
  • Gauge freedom is removed by assigning gauge conditions. The fully reduced Hamiltonian H F R is obtained from H E by using Dirac brackets to impose all constraints and gauge conditions.
The theory defined by the singular Lagrangian (1) provides a relatively complete example of each step in the algorithm.
One reason the Dirac–Bergmann algorithm can be confusing is that typical examples are chosen for simplicity, allowing some of the logical steps to be skipped. This causes the distinction between Hamiltonians to become blurred. For example, if there are no restrictions on the Lagrange multipliers, then the primary Hamiltonian H P and the total Hamiltonian H T coincide. Likewise, if there are no secondary (or higher-order) first class constraints, then the total Hamiltonian H T and the extended Hamiltonian H E coincide. Moreover, for theories with no second class constraints and no gauge conditions imposed, Dirac brackets and the reduction process are not needed.
Another confusing aspect of the Dirac–Bergmann algorithm is that for many important theories, the Lagrangian is given to us in a form that contains Lagrange multipliers. For example, consider the Einstein–Hilbert action of general relativity. The Lagrangian density is the spacetime curvature scalar. A 3 + 1 splitting of the spacetime metric [31,32] allows us to write the Lagrangian density as
L = R + K i j ( g i k g j g i j g k ) K k ,
apart from a total derivative term that integrates to the boundary. Here, R and g i j are the spatial scalar curvature and spatial metric. In addition, K i j is the extrinsic curvature of space, built from the lapse function N, shift vector N i and spatial derivatives of g i j . The Lagrangian density depends on the Lagrange multipliers N and N i as well as the configuration space coordinates g i j .
Because the Einstein–Hilbert Lagrangian (112) depends on Lagrange multipliers, it is not analogous to the singular Lagrangian (1) of our example problem. Rather, it is analogous to the partially reduced Lagrangian that appears in the integrand of the action S P R [ Q , γ ] of Equation (89). Recall that the partially reduced Hamiltonian H P R contains only first class constraints and no restrictions on the Lagrange multipliers. Likewise, the Hamiltonian for general relativity [32] is constructed from first class constraints (the Hamiltonian and momentum constraints), and the Lagrange multipliers (the lapse function and shift vector) are unrestricted.
We can attempt to eliminate the lapse and shift from the Einstein–Hilbert action, just as we eliminated the γ ’s from S P R [ Q , γ ] and obtained the result S P R [ Q ] in Equation (91). It is straightforward to eliminate the lapse N—the result is the Baierlein–Sharp–Wheeler action [33]. However, the shift vector cannot be eliminated algebraically because the equations of motion obtained by varying the action with respect to the lapse and shift depend on spatial derivatives of N i .
So, for general relativity, we do not have a singular Lagrangian analogous to Equation (1), and we cannot expect to apply the Dirac–Bergmann algorithm from beginning to end as laid out by Dirac [10]. Nevertheless, the general Dirac–Bergmann algorithm serves as the foundation for our understanding and interpretation of the Hamiltonian form of the theory.

Funding

This research received no external funding.

Acknowledgments

I would like to thank Claudio Bunster and Marc Henneaux for sparking my interest in this subject many years ago. I also thank Ashiqul Islam Dip for helpful discussions.

Conflicts of Interest

The author declares no conflict of interest.

Notes

1
This example does not cover every contingency. In particular, it does not include redundant constraints [11].
2
Some authors use the term “secondary constraints” to refer to all higher–order constraints—that is, all constraints beyond the primary level.
3
This assumes Lagrange’s equations are self–consistent. Otherwise, the consistency conditions could lead to a contradiction such as 1 = 0 .
4
For systems with more than one primary first class constraint, the Poisson brackets of any two primary first class constraints is also a first class constraint that generates a gauge transformation [10,11].
5
Counterexamples to the Dirac conjecture are discussed in Refs. [11,20,21,22,23,24,25,26] and elsewhere. Proofs of the conjecture have been constructed by adopting various simplifying assumptions [11,27,28]. The status of the conjecture is a subtle issue; see for example Refs. [18,29,30].
6
Any action can be reduced by using the equations of motion obtained by varying with respect to a subset of variables, solving those equations for the same subset of variables, then substituting the results into the action. In general, it is not permissible to reduce an action by using the equations obtained by varying with respect to one subset of variables but solving those equations for a different subset of variables.
7
Canonical gauges restrict the phase space variables. Noncanonical gauges [11] involve the Lagrange multipliers.

References

  1. Bergmann, P. Non–Linear Field Theories. Phys. Rev. 1949, 75, 680. [Google Scholar] [CrossRef]
  2. Bergmann, P.; Brunings, J. Non–Linear Field Theories II. Canonical Equations and Quantization. Rev. Mod. Phys. 1949, 21, 480. [Google Scholar] [CrossRef] [Green Version]
  3. Dirac, P. Generalized Hamiltonian dynamics. Can. J. Math. 1950, 2, 129. [Google Scholar] [CrossRef]
  4. Bergmann, P.; Penfield, R.; Schiller, R.; Zatzkis, H. The Hamiltonian of the General Theory of Relativity with Electromagnetic Field. Phys. Rev. 1950, 80, 81. [Google Scholar] [CrossRef]
  5. Dirac, P. The Hamiltonian form of field dynamics. Can. J. Math. 1951, 3, 1. [Google Scholar] [CrossRef]
  6. Anderson, J.; Bergmann, P. Constraints in Covariant Field Theories. Phys. Rev. 1951, 83, 1018. [Google Scholar] [CrossRef]
  7. Bergmann, P.; Goldberg, I. Dirac bracket transformations in phase space. Phys. Rev. 1955, 98, 531. [Google Scholar] [CrossRef]
  8. Dirac, P. Generalized Hamiltonian dynamics. Proc. R. Soc. (Lond.) 1958, 246, 326. [Google Scholar] [CrossRef]
  9. Dirac, P. The theory of gravitation in Hamiltonian form. Proc. R. Soc. (Lond.) 1958, 246, 333. [Google Scholar]
  10. Dirac, P. Lectures on Quantum Mechanics; Belfer Graduate School of Sciences, Yeshiva University: New York, NY, USA, 1964. [Google Scholar]
  11. Henneaux, M.; Teitelboim, C. Quantization of Gauge Systems; Princeton University Press: Princeton, NJ, USA, 1992. [Google Scholar]
  12. Rosenfeld, L. Zur Quantelung der Wellenfelder. Ann. Phys. 1930, 397, 113. [Google Scholar] [CrossRef]
  13. Rosenfeld, L. On the quantization of wave fields. Eur. Phys. J. H 2017, 42, 63. [Google Scholar] [CrossRef]
  14. Penfeld, R. Hamiltonians without Parametrization. Phys. Rev. 1951, 84, 737. [Google Scholar] [CrossRef]
  15. Salisbury, D.C. Peter Bergmann and the invention of constrained Hamiltonian dynamics. Einstein Stud. 2012, 12, 247–257. [Google Scholar] [CrossRef] [Green Version]
  16. Hanson, A.; Regge, T.; Teitelboim, C. Constrained Hamiltonian Systems; Accademia Nazionale dei Lincei: Roma, Italy, 1976. [Google Scholar]
  17. Sundermeyer, K. Constrained Dynamics; Springer: Berlin/Heidelberg, Germany, 1982. [Google Scholar]
  18. Rothe, H.; Rothe, K. Classical and Quantum Dynamics of Constrained Hamiltonian Systems; World Scientific: Singapore, 2010. [Google Scholar]
  19. Lusanna, L. Non-Inertial Frames and Dirac Observables in Relativity; Cambridge University Press: Cambridge, UK, 2019. [Google Scholar]
  20. Cawley, R. Determination of the Hamiltonian in the Presence of Constraints. Phys. Rev. Lett. 1979, 42, 413. [Google Scholar] [CrossRef]
  21. Li, Z. On the Invalidity of a Conjecture of Dirac. Chin. Phys. Lett. 1993, 10, 68–70. [Google Scholar] [CrossRef]
  22. Wang, Y.L.; Wang, Z.P.; Wang, K. Gauge Symmetries and the Dirac Conjecture. Int. J. Theor. Phys. 2009, 48, 1894–1904. [Google Scholar] [CrossRef]
  23. Frenkel, A. Comment on Cawley’s counterexample to a conjecture of Dirac. Phys. Rev. D 1980, 21, 2986–2987. [Google Scholar] [CrossRef]
  24. Sugano, R.; Kimura, T. Comments on Constraints of Gauge Theories. J. Phys. A 1983, 16, 4417. [Google Scholar] [CrossRef]
  25. Wu, B.C. Dirac’s conjecture. Int. J. Theor. Phys. 1994, 33, 1529–1533. [Google Scholar] [CrossRef]
  26. Miskovic, O.; Zanelli, J. Dynamical structure of irregular constrained systems. J. Math. Phys. 2003, 44, 3876–3887. [Google Scholar] [CrossRef] [Green Version]
  27. Costa, M.E.V.; Girotti, H.O.; Simoes, T.J.M. Dynamics of Gauge Systems and Dirac’s Conjecture. Phys. Rev. D 1985, 32, 405. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Cabo, A.; Louis-Martinez, D. On Dirac’s conjecture for Hamiltonian systems with first and second class constraints. Phys. Rev. D 1990, 42, 2726–2735. [Google Scholar] [CrossRef]
  29. Rothe, H.; Rothe, K. Gauge identities and the Dirac conjecture. Ann. Phys. 2004, 313, 479–496. [Google Scholar] [CrossRef] [Green Version]
  30. Pons, J. On Dirac’s incomplete analysis of gauge transformations. Stud. Hist. Phil. Sci. B 2005, 36, 491–518. [Google Scholar] [CrossRef] [Green Version]
  31. Misner, C.; Thorne, K.; Wheeler, J. Gravitation; W.H. Freeman: San Francisco, CA, USA, 1973. [Google Scholar]
  32. Arnowitt, R.L.; Deser, S.; Misner, C.W. The Dynamics of general relativity. Gen. Rel. Grav. 2008, 40, 1997–2027. [Google Scholar] [CrossRef] [Green Version]
  33. Baierlein, R.; Sharp, D.; Wheeler, J. Three-Dimensional Geometry as Carrier of Information about Time. Phys. Rev. 1962, 126, 1864–1865. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Brown, J.D. Singular Lagrangians, Constrained Hamiltonian Systems and Gauge Invariance: An Example of the Dirac–Bergmann Algorithm. Universe 2022, 8, 171. https://doi.org/10.3390/universe8030171

AMA Style

Brown JD. Singular Lagrangians, Constrained Hamiltonian Systems and Gauge Invariance: An Example of the Dirac–Bergmann Algorithm. Universe. 2022; 8(3):171. https://doi.org/10.3390/universe8030171

Chicago/Turabian Style

Brown, J. David. 2022. "Singular Lagrangians, Constrained Hamiltonian Systems and Gauge Invariance: An Example of the Dirac–Bergmann Algorithm" Universe 8, no. 3: 171. https://doi.org/10.3390/universe8030171

APA Style

Brown, J. D. (2022). Singular Lagrangians, Constrained Hamiltonian Systems and Gauge Invariance: An Example of the Dirac–Bergmann Algorithm. Universe, 8(3), 171. https://doi.org/10.3390/universe8030171

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop