Faddeev-Jackiw Hamiltonian Reduction for Free and Gauged Rarita-Schwinger Theories

We study the Faddeev-Jackiw symplectic Hamiltonian reduction for 3+1-dimensional free and Abelian gauged Rarita-Schwinger theories that comprise Grassmannian fermionic fields. We obtain the relevant fundamental brackets and find that they are in convenient forms for quantization. The brackets are independent of whether the theories contain mass or gauge fields, and the structure of constraints and symplectic potentials largely determine characteristic behaviors of the theories. We also note that, in contrast to the free massive theory, the Dirac field equations for free massless Rarita-Schwinger theory cannot be obtained in a covariant way.


I. INTRODUCTION
In 1941, Rarita and Schwinger constructed a theory of spin- 3 2 vector-spinor fields which has a local fermionic gauge-invariance [1]. However, this symmetry is lost when the vector-spinor field has mass or couples to the other lower spin fields. More precisely, in 1961, Johnson and Sudarshan studied massive Rarita-Schwinger field minimally coupled to an external electromagnetic field, and showed that the equal-time commutators and relativistic covariance of the theory are in conflict, which makes the quantization a rather subtle issue [2]. In 1969, Velo and Zwanziger found that the massive gauged extension of the theory also admits superluminal wave propagation. Thus, the causality principle is also violated in the theory [3]. Despite these persistent problems, the massless theory keeps its importance particularly in two aspects. First, the massless (Majorana) Rarita-Schwinger field plays a central role in the construction of covariantly interacting supergravity theory [4][5][6]. The theory describes a generalization of the Rarita-Schwinger fermionic gaugeinvariance and the vector-spinor fields are fermionic superpartner of gravitons, namely gravitinos of the supergravity. In this concept, Das and Freedman showed that the massless theory is free from the non-causal wave propagation and has a unitary propagator structure [7]. Secondly, the massless Rarita-Schwinger theory is valuable for the cancellation of SU(8) gauge anomalies. Unlike the generic anomaly cancellation mechanisms in which the anomalies are supposed to be canceled withing the lower spin fermionic fields, it was shown by Marcus [8] and later studied by Adler [9], that a complete SU (8) gauge theory can be constructed via Rarita-Schwinger fields. In this set-up, the vector-spinor field acquires a crucial role in canceling anomalies arising in the gauge theory. Thus, it is left to determine whether the gauged Rarita-Schwinger fields describe well-behaved, complete classical or quantum field theories. For this purpose, Adler has recently studied minimally gauged massless Rarita-Schwinger theories at both classical and quantum levels in detail [10]. He showed that, unlike the massive case, the massless gauged Rarita-Schwinger theory provides consistent classical and quantum theories with a generalized fermionic gauge-invariance.
Taking the above mentioned observations as inspiration points and noting the hard task of getting proper brackets of constrained systems providing viable quantization, we study the Faddeev-Jackiw (FJ) symplectic Hamiltonian reduction [11,12] for free and gauged Rarita-Schwinger theories. Unlike Dirac's approach for constrained systems [16], FJ symplectic first-order formalism does not require any classification of constraints 1 . In other words, the method avoids analyzing systems by evaluating all commutation relations among the constraints and classifying them accordingly. Apparently, the FJ approach supplies a rather economical way of quantizing constrained systems. In doing so, we find the fundamental brackets for the free and gauged Rarita-Schwinger theories for both massless and massive versions. Here, the brackets are in admissible structures to be quantized. We also observe that the brackets are identical for all kinds of the theories; the brackets are independent of whether the theory is massive or interacting with external electromagnetic field or not. The differences between the theories arise among the constraints they have. We also notice that, in contrast to the massive case, the Dirac field equations for free massless Rarita-Schwinger theory cannot be obtained in a covariant way.
The layout of the paper is as follows: In Sec. II, we recapitulate the fundamental properties of free massless Rarita-Schwinger theory and apply FJ Hamiltonian reduction to the theory. In Sec. III, we turn our attention to the FJ Hamiltonian reduction for free massive Rarita-Schwinger theory. Sec. IV and Sec. V are devoted to the first-order symplectic analysis for Abelian gauged extensions of massless and massive Rarita-Schwinger theories. In Sec. VI, we conclude our results. In the Appendix A, the derivation of the transverse and traceless decomposition of the fields in the free massless Rarita-Schwinger theory is given as a sample. In the Appendix B, we briefly review the FJ approach for constrained and unconstrained systems. We also give an example of the application of symplectic method to anti-commuting spin-1 2 Dirac theory.

II. FREE MASSLESS RARITA-SCHWINGER THEORY
The 3 + 1-dimensional free massless Rarita-Schwinger theory is described by the Lagrangian where ψ µ andψ µ are vector-spinor fields with spinor indices suppressed. We work in the metric signature (+, −, −, −), γ 5 = iγ 0 γ 1 γ 2 γ 3 and {γ µ , γ ν } = 2η µν . We consider the fermionic fields as independent anti-commuting Grassmannian variables. Recall that, unlike the complex Dirac field, for the Grassmannian variables there is no such relation asψ µ = γ 0 ψ + µ . Instead, ψ µ andψ µ are independent generators in the Grassmann algebra. Thus, one can define the conjugation as follows: Notice that this does not mean that Eq.(2) produces a new element in the Grassmannian algebra. This is merely the conjugation of independent variables. Therefore, with the help of the conjugation of the Grassmannian variables (θ 1 θ 2 ) * = θ * 2 θ * 1 , one can show that the Lagrangian in Eq.(1) is selfadjoint up to a boundary term: such that the total derivative term naturally drops at the action level. Moreover, variations with respect to independent variables respectively yield which are the corresponding field equations. From now on, we will work with the first of Eq.(4). But, by following the same steps, one could easily obtain the similar results for the second equation. Notice that by using the identity one can recast the field equation in Eq.(4) as follows Here / ∂ = γ µ ∂ µ and γ · ψ = γ µ ψ µ . Contracting Eq.(6) with γ λ gives Finally, by plugging this result in Eq.(6), the field equation reduces to To obtain the real propagating degrees of freedom, let us now study gauge transformation and corresponding gauge conditions. For this purpose, let us recall that under the local Rarita-Schwinger fermionic gauge transformation the Lagrangian in Eq.(1) transforms as Here ǫ is an arbitrary four-component spinor field. As is seen in Eq.(10), the free massless Rarita-Schwinger Lagrangian changes by a total derivative under the Rarita-Schwinger gauge transformation, which drops at the action level and thus we have a completely gauge-invariant theory. This means that the theory admits a gauge redundancy. To find the correct physical degrees of freedom of the theory, one needs to fix this gauge-freedom. For this purpose, let us assume the Coulomb-like gauge condition where i = 1, 2, 3. In fact, this is a reasonable gauge choice: Any initial data ψ ′ i (x, t) that does not satisfy Eq. (11) can be tuned to the desired form via 2 (See [7] and [17] for further discussions). For the sake of the self-completeness, one needs to examine the theory further to see whether Eq.(11) imposes any additional conditions or not. For this purpose, note that ψ 0 component does not have a time derivative, so it is a Lagrange multiplier. In other words, as in the electromagnetic case, the zeroth component of the vector-spinor field is a zero mode which is followed with a constraint. More precisely, the λ = 0 component of the field equation in Eq.(8) reads One can also get a secondary constraint by contracting the field equation with ∂ λ . But since our primary aim is not analyzing the system by examining all the existing constraints, we leave it as a comment. As is seen in Eq.(13), gauge fixing condition γ i ψ i = 0 imposes γ i ∂ i ψ 0 = 0. Here, since the operator is not invertible, we are not allowed to get ψ 0 = 0 as a corollary of γ i ψ i = 0; yet we assume an additional condition of ψ 0 = 0. Furthermore, splitting the fully contracted equation in Eq.(7) into its space and time components yields In Eq. (14), one should notice that the gauge fixing condition γ i ψ i = 0 together with the assumed condition ψ 0 = 0 impose ∂ i ψ i = 0. As a consequence of this, we obtain the set of consistency conditions Observe that Eq. (15) can also be written in covariant forms as follows

Symplectic Reduction for Free Massless Rarita-Schwinger Theory
In this section, we study the FJ Hamiltonian reduction for the free massless Rarita-Schwinger theory which will lead us to the fundamental brackets of the theory. For this purpose, let us recast the Lagrangian in Eq.(1) in a more symmetric form: To study the theory in the first-order symplectic formalism, one needs to convert Eq.(18) into the desired symplectic form. That is, one needs to split the Lagrangian into its space and time components. After a straightforward decomposition, one gets where the symplectic coefficients are and the corresponding symplectic potential reads As expected, all the non-dynamical components have been relegated into the Hamiltonian part of the system. In analyzing the theory, one could also choose the conjugate momenta ofψ k as a dynamical variable. But in our analysis, we will not work with it. Instead, we consider ψ µ and ψ µ as the independent variables. Note that ψ 0 andψ 0 are not dynamical components, so they are Lagrange multipliers. Following [11,12], the elimination of constraints give the equations To solve the constraint equations, one can decompose the independent fields into its local transverse and γ-traceless parts as where "T " and "ˆ" stand for the transverse and traceless parts, respectively. Here the γ-traceless parts areψ such that γ iψ i = 0 and γ iψ i = 0. Then, by using the identity as well as the constraints in Eq. (22), one can show that the transverse and traceless decomposition of the fields in Eq.(23) can actually be written as follows As a side comment, one should note that as is done in [12], without addressing the transverse and γ-traceless parts (23), one could also directly start with the (26). Here, we further provide what the explicit form of the Longitudinal part is. (See Appendix A for the derivation of Eq.(26)). Accordingly, the constraint equations in Eq.(22) turn into completely transverse ones Finally, by inserting Eq.(26) and Eq.(27) in the Eq. (19), up to a boundary term, one gets a completely transverse Lagrangian Here the transverse symplectic coefficients and potential are Thus, by defining the symplectic variables as (ξ 1 , ξ 2 ) = (ψ T k ,ψ T k ), one gets the corresponding symplectic matrix which is clearly non-singular. Notice that the minus sign in the sub-block is due to the antisymmetric ǫ tensor. Therefore, by taking care of the epsilons contraction in the current signature, one can easily show that the inverse symplectic matrix is Once the inverse symplectic matrix is found, one can evaluate the fundamental brackets. That is, by using the definition of the FJ equal-time brackets for the Grassmann variables one gets the fundamental brackets for free massless Rarita-Schwinger theory as follows Note that, with the help of the identity in Eq.(25), the non-vanishing bracket can also be rewritten as which is identical with the one found in [18].

III. FREE MASSIVE RARITA-SCHWINGER THEORY
The Lagrangian that describes the 3 + 1-dimensional free massive Rarita-Schwinger theory is where Recall that the fermionic fields are anti-commuting Grassmannian variables. Accordingly, the field equations of the independent variables respectively read In dealing with the fundamental properties of the theory, as we did in the massless theory, we will work only with the first field equation in Eq.(34). Notice that by using the identity in Eq.(5), one can recast the field equation as follows Observe that the contraction of Eq.(35) with γ λ yields and the contraction of Eq.
Combining both contracted field equations Eq.(36) and Eq.(37), one obtains With these gauge-fixing conditions, the equation in Eq.(35) turns into the Dirac field equation for massive spin-3 2 vector-spinor field Note that, unlike the massless theory, one obtains the Dirac field equation in Eq.(39) without addressing the space and time decompositions of the field equations. On the other hand, due to the mass term, the Rarita-Schwinger gauge-invariance is inevitably lost.

Symplectic Reduction for Free Massive Rarita-Schwinger Lagrangian
Let us now study the symplectic Hamiltonian reduction of the free massive Rarita-Schwinger theory. For this purpose, let us recall that the Lagrangian in Eq.(33), up to a boundary term, can be written as In order to proceed the FJ symplectic reduction of Eq.(40), one needs to separate the dynamical components from the non-dynamical ones so that the non-dynamical components can be relegated to Hamiltonian part of the Lagrangian. Therefore, by splitting the Lagrangian into its space and time components, one will obtain where the coefficient of the dynamical parts are and the explicit form of the symplectic potential is Like the free massless theory, ψ 0 andψ 0 are zero modes of the system whose eliminations give rise the constraints As was done in the previous section, by decomposing the fields into the local transverse and γtraceless parts as in the Eq.(23), the constraints in Eq.(44) turn into completely transverse ones In this case, the longitudinal part reads ζ = ( / ∂ + im)γ · ψ T . Thus, by plugging the Eq.(23) and the transverse constraints Eq.(45) into the Eq.(41), up to a boundary term, the Lagrangian turns into where the transverse symplectic coefficients and potential respectively are Observe that the middle two terms in Eq.(46) are not in the symplectic forms. Therefore, by assuming the Darboux transformation with an additional assumption of the undesired terms in Eq.(46) drop and thus we are left with a completely transverse Lagrangian Note that the extra condition Eq.(49) is enforced by the Darboux transformation and the constraint equations; otherwise, the coupled terms in the symplectic part could not be decoupled. In fact, it seems there is a lack in the physical interpretation of Eq.(49). Therefore, it will be particularly interesting if one can show that it has a relation with the real constraints or not. Here, as is mentioned in Eq.(117), the remaining variables (i.e., the longitudinal components) are denoted as the Lagrange multipliers such that As noted in [11,12], since the last two terms in the Eq.(50) cannot be dropped via elimination of constraints anymore, Eq.(52) corresponds to the true constraints of the system. Note also that the true constraints cannot be rewritten as linear combinations of the ones that are obtained during the eliminations of the constraints; otherwise, they would also drop when the eliminations of constraint was performed. These are the constraints that cannot be eliminated anymore. Therefore, setting φ k (ψ T k ,ψ T k ) andφ i (ψ T k ,ψ T k ) to zero provides an unconstrained fully traceless Lagrangian Thus, with the definition of the dynamical variables (ξ 1 , ξ 2 ) = (ψ T k ,ψ T k ), the non-vanishing equaltime FJ bracket for the free massive Rarita-Schwinger theory becomes which is same as the one found in [19].

IV. GAUGED MASSLESS RARITA-SCHWINGER THEORY
In this section, we study the massless Rarita-Schwinger field minimally coupled to an external electromagnetic field which is described by the Lagrangian Here the gauge covariant derivative is D ν = ∂ ν + gA ν , where g is the relevant coupling constant and A µ is an Abelian gauge field. The field equations read As in the free massless and massive theories, while deducing the some basic properties of the theory, we will only deal with the first of Eq.(56). Notice that with the help of the identity in Eq.(5), the Eq.(56) turns into Moreover, contracting the Eq.(57) with γ λ yields Finally, substituting the Eq.(58) in Eq. (57) gives On the other side, contracting the Eq.(56) with D λ becomes which is a secondary constraint in the theory and does not provide any further simplification in the field equation in Eq.(59).

Symplectic Reduction for Gauged Massless Rarita-Schwinger Theory
Let us now apply the first-order symplectic formalism to the massless Rarita-Schwinger fields minimally coupled to an external electromagnetic field. For this purpose, let us note that the Lagrangian of the theory in Eq.(55) can be recast in a more symmetric form as follows Similarly, by splitting the Lagrangian in Eq.(61) into its space and time components, one gets where the symplectic coefficients are and the related symplectic potential is Note that although the gauge fields are non-dynamical variables, due to being external potentials, one cannot vary and then impose these variations to be vanished. Otherwise, as in the Quantum Electromagnetic Dynamics with external potential, the gauge field current would be enforced to be zero which is not a desired situation. Hence, as in the free theories, here ψ 0 andψ 0 are the only zero modes of the theory: Therefore, variations with respect to ψ 0 andψ 0 respectively give the following constraint equations As was done in the free theories, by decomposing the fields into the local transverse and γ-traceless parts as in Eq.(23) 3 and using the constraints in Eq.(65) as well as by assuming the Darboux transformation (48), with an additional assumption of the Lagrangian (62) turns into a completely transverse one where the transverse symplectic coefficients and potential read Note that the symplectic potential also contains gauge field parts. Furthermore, as is given in (117), the remaining variables (i.e., the longitudinal components) are denoted as the Lagrange multipliers such that which cannot be dropped via elimination of constraints anymore so, according to [11,12], they are the true constraint of the system. Thus, by setting φ k (ψ T k ,ψ T k ) andφ i (ψ T k ,ψ T k ) to zero, one arrives at a completely transverse Lagrangian Finally, with the definition of the symplectic dynamical variables (ξ 1 , ξ 2 ) = (ψ T k ,ψ T k ), one obtains the non-vanishing equal-time FJ basic bracket for the gauged massless Rarita-Schwinger theory as follows which is consistent with the Pauli-spin-part of the fundamental bracket obtained in [10] in which Adler studies the Dirac quantization of the non-Abelian gauged Rarita-Schwinger theory via the left-chiral component of the fermionic field. One should notice that such a difference is expected because in [10], the corresponding gauge fields are non-Abelian variables; however here the gauge fields are Abelian vector fields.

V. GAUGED MASSIVE RARITA-SCHWINGER
In this section, we study the massive Rarita-Schwinger theory minimally coupled to an external electromagnetic field which is described by the Lagrangian where the gauge-covariant derivative is D ν = ∂ ν + gA ν . Accordingly, the field equations for the independent anti-commuting fermionic fields are which with the help of the identity in Eq.(5) turns into Moreover, contraction of the equation in Eq.(75) with γ λ gives And contraction of field equation in Eq.(74) with D λ becomes which with the additional redefinition turns into Combining Eq.(76) and Eq.(79), one gets the secondary constraint that determines the equation of motion of ψ 0 component as follows Observe that using Eq.(80) in Eq.(79) gives the relation Finally, by plugging Eq.(80) and Eq.(81) into the field equation in Eq.(75), one obtains which is the equation that is used by Velo and Zwanziger in deducing the acausal wave propagation of the solution by finding the future-directed normals to the surfaces at each point [3].

Symplectic Reduction for Gauged Massive Rarita-Schwinger Theory
Finally, let us apply FJ symplectic Hamiltonian reduction to the massive Rarita-Schwinger field minimally coupled to an external electromagnetic field. In order to do so, let us rewrite the Lagrangian in Eq.(73) in a more symmetric form: Subsequently, by splitting Lagrangian in Eq.(83) into its space and time components, one gets where the symplectic coefficients are and the relevant Hamiltonian Note that as is emphasized in the massless gauged part, since the gauge fields are external potentials, one is not allowed to set their variation to zero. Hence, here ψ 0 ,ψ 0 are the only Lagrange multipliers that induce constraints on the system. Therefore, eliminations of constraint yield Like the free massive theory, by decomposing the dynamical components into the local transverse and traceless parts as in Eq.(23) 4 as well as using constraints in Eq.(87) and the Darboux transformation (48), with an additional assumption of the Lagrangian, up to a boundary term, turns into Here the transverse symplectic coefficients and potential are (91) 4 In this case, from the constraint equation, one finds Notice that, different from the free cases, the symplectic potential involves mass and gauge potentials. Here in the Eq.(90), as in the previous sections, the Lagrange multipliers are the Longitudinal parts of the vector-spinor field and the corresponding constraints read which are same as Eq.(70). Similarly, by setting Eq.(92) to zero [11,12], one arrives at a completely transverse Lagrangian whose symplectic part is same as the ones found so far. Thus, with the definition of the dynamical variables (ξ 1 , ξ 2 ) = (ψ T k ,ψ T k ), the non-vanishing equal-time bracket for the gauged massive Rarita-Schwinger theory becomes which is identical to the one found in [20].

VI. CONCLUSIONS
In this work, we studied 3 + 1-dimensional free and Abelian gauged Grassmannian Rarita-Schwinger theories for their massless and massive extensions in the context of Faddeev-Jackiw first-order symplectic formalism. We have obtained the fundamental brackets of theories which are consistent with the some results that we found in the literature but obtained in a more simpler way. The brackets are independent of whether the theories contain mass or gauge field or not, and thus the structure of constraints and symplectic potentials determine characteristic behaviors of the theories. It will be particularly interesting to find proper transformations that will relate the constraints obtained via the Faddeev-Jackiw symplectic method with the ones that are obtained via Dirac method. But since the constraints obtained in both methods are rather complicated, in this paper, we restrict ourselves only to the Faddeev-Jackiw analysis of Rarita-Schwinger theories and leave this as a future work. With the comparison with the literature, we concluded that the Faddeev-Jackiw symplectic approach provides a more economical way in deriving the fundamental brackets for the Rarita-Schwinger theories. In addition to these, we notice that, in contrast to the massive theory, the Dirac field equations for free massless Rarita-Schwinger theory cannot be covariantly deduced.

VIII. APPENDIX A: TRANSVERSE AND TRACELESS DECOMPOSITION OF FIELDS
In this section, let us give the derivations of (26) and (27): To solve the constraint equations, one can decompose the independent fields into its local transverse and γ-traceless parts as where "T " and "ˆ" stand for the transverse and traceless parts, respectively. Here the γ-traceless parts areψ Therefore, we have To find how the constraint equations in Eq.(22) decomposes under Eq.(95), let us focus on the following constraint equation Note that with the identity ǫ ijk γ 5 γ k = −γ 0 σ ij and Eq.(95), the Eq.(98) turns into Furthermore, by using the relation and the transverse and traceless properties of the fields in Eq.(97), one gets Notice that after contraction with iγ 0 and relabeling of the dummy indices, it becomes This structure is also valid for the other theories. The only difference arises in the definition of ζ which we give its explicit form in each section. Finally, by substituting this result into the constraint in Eq.(98), it turns into Because of the symmetric and anti-symmetric contraction in "j, k" indices, the second term drops, and we are left with the transverse constraint equation

IX. APPENDIX B: FADDEEV-JACKIW HAMILTONIAN REDUCTION FOR CONSTRAINED AND UNCONSTRAINED SYSTEMS
In this section, we review the Faddeev-Jackiw symplectic first-order formalism which was introduced particularly to quantize the constrained systems [11,12]. The method works on the first-order Lagrangian and does not require any classification of constraints. To better understand how the method works, let us consider With the definition of 2n-component phase-space coordinates ξ α = p α , α = 1, · · · , n and ξ β = q β , β = n + 1, · · · , 2n, Eq.(106) can be rewritten as a Lagrangian one-form Here the symplectic 2n × 2n matrix is [11,12]. But, in general, the symplectic twoform does not have to be constant. Therefore, let us now consider the following generic Lagrangian where A α is an arbitrary one-form. The variation of Eq.(109) with respect to ξ yields In the case of when the symplectic matrix is nonsingular, Eq.(110) becomeṡ Thus, by using Eq.(111) and the Poisson brackets for the bosonic variables, one obtains the FJ fundamental brackets as follows Note that, in the case of the Grassmannian variables, using the anti-commutation property of the variables as well as the Poisson brackets for the Grassmannian variables [21], one haṡ and the corresponding fundamental brackets become On the other side, when there are constraints in the system which are induced by the existence of the zero-modes, then the symplectic matrix cannot be inverted. In that case, according to the Darboux's theorem which states that for any given one-form A = A α dξ α where α = 1, · · · , N , one can always do the following changes in the variables ξ α → (p β , q γ z ρ ), β, γ = 1, · · · , n, ρ = 1, · · · , N − 2n, so that A turns into A = A α dq α . As is seen above, when there is no constraint, Eq.(115) diagonalizes f αβ . However when there are constraints, only a 2n × 2n sub-block of f αβ diagonalizes and the remaining N − 2n degrees of freedom (corresponding the zero-modes z ρ ) will not be in the symplectic form [11,12]; yet they occur in the rest of the Lagrangian: The equations ∂Φ ∂z α = 0 can be used to eliminate the zero-modes of z's only if ∂ 2 Φ ∂z ρ ∂z β is nonsingular. In the generic case, after diagonalization and elimination of z's as many as possible, one ultimately arrives at where the remaining z's are denoted by λ ρ (namely, Lagrange multipliers) and the φ ρ are the only true constraints in the system A. Symplectic Reduction for Dirac Theory of spin- 1 2 fields In this section, to see how the method works, we provide FJ Hamiltonian reduction for the Dirac theory for the spin-1 2 theory as an example. For this purpose, let us note that the Lagrangian can be written: As mentioned above, we assume that the independent dynamical variables are anti-commuting Grassmann variables. In order to pass to the symplectic analysis of the system, one needs to separate the dynamical components from the non-physical ones by splitting the Lagrangian (119) into its time and space components. In doing so, one arrives at whose variation, up to a boundary term, yields from which one gets the Dirac field equations as follows As is seen from Eq.(121) and Eq.(122), the symplectic matrix for the Dirac theory and its inverse are One should observe that, in contrast to the bosonic case, the symplectic matrix for the Grassmannian variables is symmetric and the fundamental brackets are defined as follows from which one gets the basic bracket for the Dirac theory This is also valid for the massless theory. Note that since the theory does not have any gauge redundancy, one does not need to assume any gauge-fixing.