BRST quantization

BRST quantization is a differential geometric approach to performing consistent, anomaly-free perturbative calculations in a non-abelian gauge theory. The analytical form of the BRST "transformation" and its relevance to renormalization and anomaly cancellation were described by Carlo Maria Becchi, Alain Rouet, and Raymond Stora in a series of papers culminating in the 1976 "Renormalization of gauge theories". The equivalent transformation and many of its properties were independently discovered by Igor Viktorovich Tyutin. Its significance for rigorous canonical quantization of a Yang–Mills theory and its correct application to the Fock space of instantaneous field configurations were elucidated by Taichiro Kugo and Izumi Ojima. Later work by many authors, notably Thomas Schücker and Edward Witten, has clarified the geometric significance of the BRST operator and related fields and emphasized its importance to topological quantum field theory and string theory.

In the BRST approach, one selects a perturbation-friendly gauge fixing procedure for the action principle of a gauge theory using the differential geometry of the gauge bundle on which the field theory lives. One then quantizes the theory to obtain a Hamiltonian system in the interaction picture in such a way that the "unphysical" fields introduced by the gauge fixing procedure resolve gauge anomalies without appearing in the asymptotic states of the theory. The result is a set of Feynman rules for use in a Dyson series perturbative expansion of the S-matrix which guarantee that it is unitary and renormalizable at each loop order—in short, a coherent approximation technique for making physical predictions about the results of scattering experiments.

This is related to a supersymplectic manifold where pure operators are graded by integral ghost numbers and we have a BRST cohomology.

From a practical perspective, a quantum field theory consists of an action principle and a set of procedures for performing perturbative calculations. There are other kinds of "sanity checks" that can be performed on a quantum field theory to determine whether it fits qualitative phenomena such as quark confinement and asymptotic freedom. However, most of the predictive successes of quantum field theory, from quantum electrodynamics to the present day, have been quantified by matching S-matrix calculations against the results of scattering experiments.

In the early days of QFT, one would have had to say that the quantization and renormalization prescriptions were as much part of the model as the Lagrangian density, especially when they relied on the powerful but mathematically ill-defined path integral formalism. It quickly became clear that QED was almost "magical" in its relative tractability, and that most of the ways that one might imagine extending it would not produce rational calculations. However, one class of field theories remained promising: gauge theories, in which the objects in the theory represent equivalence classes of physically indistinguishable field configurations, any two of which are related by a gauge transformation. This generalizes the QED idea of a local change of phase to a more complicated Lie group.

QED itself is a gauge theory, as is general relativity, although the latter has proven resistant to quantization so far, for reasons related to renormalization. Another class of gauge theories with a non-Abelian gauge group, beginning with Yang–Mills theory, became amenable to quantization in the late 1960s and early 1970s, largely due to the work of Ludwig D. Faddeev, Victor Popov, Bryce DeWitt, and Gerardus 't Hooft. However, they remained very difficult to work with until the introduction of the BRST method. The BRST method provided the calculation techniques and renormalizability proofs needed to extract accurate results from both "unbroken" Yang–Mills theories and those in which the Higgs mechanism leads to spontaneous symmetry breaking. Representatives of these two types of Yang–Mills systems—quantum chromodynamics and electroweak theory—appear in the Standard Model of particle physics.

It has proven rather more difficult to prove the existence of non-Abelian quantum field theory in a rigorous sense than to obtain accurate predictions using semi-heuristic calculation schemes. This is because analyzing a quantum field theory requires two mathematically interlocked perspectives: a Lagrangian system based on the action functional, composed of fields with distinct values at each point in spacetime and local operators which act on them, and a Hamiltonian system in the Dirac picture, composed of states which characterize the entire system at a given time and field operators which act on them. What makes this so difficult in a gauge theory is that the objects of the theory are not really local fields on spacetime; they are right-invariant local fields on the principal gauge bundle, and different local sections through a portion of the gauge bundle, related by passive transformations, produce different Dirac pictures.

What is more, a description of the system as a whole in terms of a set of fields contains many redundant degrees of freedom; the distinct configurations of the theory are equivalence classes of field configurations, so that two descriptions which are related to one another by a gauge transformation are also really the same physical configuration. The "solutions" of a quantized gauge theory exist not in a straightforward space of fields with values at every point in spacetime but in a quotient space (or cohomology) whose elements are equivalence classes of field configurations. Hiding in the BRST formalism is a system for parameterizing the variations associated with all possible active gauge transformations and correctly accounting for their physical irrelevance during the conversion of a Lagrangian system to a Hamiltonian system.

The principle of gauge invariance is essential to constructing a workable quantum field theory. But it is generally not feasible to perform a perturbative calculation in a gauge theory without first "fixing the gauge"—adding terms to the Lagrangian density of the action principle which "break the gauge symmetry" to suppress these "unphysical" degrees of freedom. The idea of gauge fixing goes back to the Lorenz gauge approach to electromagnetism, which suppresses most of the excess degrees of freedom in the four-potential while retaining manifest Lorentz invariance. The Lorenz gauge is a great simplification relative to Maxwell's field-strength approach to classical electrodynamics, and illustrates why it is useful to deal with excess degrees of freedom in the representation of the objects in a theory at the Lagrangian stage, before passing over to Hamiltonian mechanics via the Legendre transformation.

The Hamiltonian density is related to the Lie derivative of the Lagrangian density with respect to a unit timelike horizontal vector field on the gauge bundle. In a quantum mechanical context it is conventionally rescaled by a factor ${\displaystyle i\hbar }$ . Integrating it by parts over a spacelike cross section recovers the form of the integrand familiar from canonical quantization. Because the definition of the Hamiltonian involves a unit time vector field on the base space, a horizontal lift to the bundle space, and a spacelike surface "normal" (in the Minkowski metric) to the unit time vector field at each point on the base manifold, it is dependent both on the connection and the choice of Lorentz frame, and is far from being globally defined. But it is an essential ingredient in the perturbative framework of quantum field theory, into which the quantized Hamiltonian enters via the Dyson series.

For perturbative purposes, we gather the configuration of all the fields of our theory on an entire three-dimensional horizontal spacelike cross section of P into one object (a Fock state), and then describe the "evolution" of this state over time using the interaction picture. The Fock space is spanned by the multi-particle eigenstates of the "unperturbed" or "non-interaction" portion ${\displaystyle {\mathcal {H}}_{0}}$  of the Hamiltonian ${\displaystyle {\mathcal {H}}}$ . Hence the instantaneous description of any Fock state is a complex-amplitude-weighted sum of eigenstates of ${\displaystyle {\mathcal {H}}_{0}}$ . In the interaction picture, we relate Fock states at different times by prescribing that each eigenstate of the unperturbed Hamiltonian experiences a constant rate of phase rotation proportional to its energy (the corresponding eigenvalue of the unperturbed Hamiltonian).

Hence, in the zero-order approximation, the set of weights characterizing a Fock state does not change over time, but the corresponding field configuration does. In higher approximations, the weights also change; collider experiments in high-energy physics amount to measurements of the rate of change in these weights (or rather integrals of them over distributions representing uncertainty in the initial and final conditions of a scattering event). The Dyson series captures the effect of the discrepancy between ${\displaystyle {\mathcal {H}}_{0}}$  and the true Hamiltonian ${\displaystyle {\mathcal {H}}}$ , in the form of a power series in the coupling constant g; it is the principal tool for making quantitative predictions from a quantum field theory.

To use the Dyson series to calculate anything, one needs more than a gauge-invariant Lagrangian density; one also needs the quantization and gauge fixing prescriptions that enter into the Feynman rules of the theory. The Dyson series produces infinite integrals of various kinds when applied to the Hamiltonian of a particular QFT. This is partly because all usable quantum field theories to date must be considered effective field theories, describing only interactions on a certain range of energy scales that we can experimentally probe and therefore vulnerable to ultraviolet divergences. These are tolerable as long as they can be handled via standard techniques of renormalization; they are not so tolerable when they result in an infinite series of infinite renormalizations or, worse, in an obviously unphysical prediction such as an uncancelled gauge anomaly. There is a deep relationship between renormalizability and gauge invariance, which is easily lost in the course of attempts to obtain tractable Feynman rules by fixing the gauge.

The traditional gauge fixing prescriptions of continuum electrodynamics select a unique representative from each gauge-transformation-related equivalence class using a constraint equation such as the Lorenz gauge ${\displaystyle \partial ^{\mu }A_{\mu }=0}$ . This sort of prescription can be applied to an Abelian gauge theory such as QED, although it results in some difficulty in explaining why the Ward identities of the classical theory carry over to the quantum theory—in other words, why Feynman diagrams containing internal longitudinally polarized virtual photons do not contribute to S-matrix calculations. This approach also does not generalize well to non-Abelian gauge groups such as the SU(2)xU(1) of Yang–Mills electroweak theory and the SU(3) of quantum chromodynamics. It suffers from Gribov ambiguities and from the difficulty of defining a gauge fixing constraint that is in some sense "orthogonal" to physically significant changes in the field configuration.

More sophisticated approaches do not attempt to apply a delta function constraint to the gauge transformation degrees of freedom. Instead of "fixing" the gauge to a particular "constraint surface" in configuration space, one can break the gauge freedom with an additional, non-gauge-invariant term added to the Lagrangian density. In order to reproduce the successes of gauge fixing, this term is chosen to be minimal for the choice of gauge that corresponds to the desired constraint and to depend quadratically on the deviation of the gauge from the constraint surface. By the stationary phase approximation on which the Feynman path integral is based, the dominant contribution to perturbative calculations will come from field configurations in the neighborhood of the constraint surface.

The perturbative expansion associated with this Lagrangian, using the method of functional quantization, is generally referred to as the R_ξ gauge. It reduces in the case of an Abelian U(1) gauge to the same set of Feynman rules that one obtains in the method of canonical quantization. But there is an important difference: the broken gauge freedom appears in the functional integral as an additional factor in the overall normalization. This factor can only be pulled out of the perturbative expansion (and ignored) when the contribution to the Lagrangian of a perturbation along the gauge degrees of freedom is independent of the particular "physical" field configuration. This is the condition that fails to hold for non-Abelian gauge groups. If one ignores the problem and attempts to use the Feynman rules obtained from "naive" functional quantization, one finds that one's calculations contain unremovable anomalies.

The problem of perturbative calculations in QCD was solved by introducing additional fields known as Faddeev–Popov ghosts, whose contribution to the gauge-fixed Lagrangian offsets the anomaly introduced by the coupling of "physical" and "unphysical" perturbations of the non-Abelian gauge field. From the functional quantization perspective, the "unphysical" perturbations of the field configuration (the gauge transformations) form a subspace of the space of all (infinitesimal) perturbations; in the non-Abelian case, the embedding of this subspace in the larger space depends on the configuration around which the perturbation takes place. The ghost term in the Lagrangian represents the functional determinant of the Jacobian of this embedding, and the properties of the ghost field are dictated by the exponent desired on the determinant in order to correct the functional measure on the remaining "physical" perturbation axes.

Intuition for the BRST formalism is provided by describing it geometrically, in the setting of fiber bundles. This geometric setting contrasts with and illuminates the older traditional picture, that of algebra-valued fields on Minkowski space, provided in (earlier) quantum field theory texts.

In this setting, a gauge field can be understood in one of two different ways. In one, the gauge field is a local section of the fiber bundle. In the other, the gauge field is little more than the connection between adjacent fibers, defined on the entire length of the fiber. Corresponding to these two understandings, there are two ways to look at a gauge transformation. In the first case, a gauge transformation is just a change of local section. In general relativity, this is referred to as a passive transformation. In the second view, a gauge transformation is a change of coordinates along the entire fiber (arising from multiplication by a group element g) which induces a vertical diffeomorphism of the principal bundle.

This second viewpoint provides the geometric foundation for the BRST method. Unlike a passive transformation, it is well-defined globally on a principal bundle, with any structure group over an arbitrary manifold. That is, the BRST formalism can be developed to describe the quantization of any principle bundle on any manifold. For concreteness and relevance to conventional QFT, much of this article sticks to the case of a principal gauge bundle with compact fiber over 4-dimensional Minkowski space.

A principal gauge bundle P over a 4-manifold M is locally isomorphic to U × F, where U ⊂ R⁴ and the fiber F is isomorphic to a Lie group G, the gauge group of the field theory (this is an isomorphism of manifold structures, not of group structures; there is no special surface in P corresponding to 1 in G, so it is more proper to say that the fiber F is a G-torsor). The most basic property as a fiber bundle is the "projection to the base space" π : P → M, which defines the vertical directions on P (those lying within the fiber π⁻¹(p) over each point p in M). As a gauge bundle it has a left action of G on P which respects the fiber structure, and as a principal bundle it also has a right action of G on P which also respects the fiber structure and commutes with the left action.

The left action of the structure group G on P corresponds to a change of coordinate system on an individual fiber. The (global) right action R_g : P → P for a fixed g in G corresponds to an actual automorphism of each fiber and hence to a map of P to itself. In order for P to qualify as a principal G-bundle, the global right action of each g in G must be an automorphism with respect to the manifold structure of P with a smooth dependence on g, that is, a diffeomorphism P × G → P.

The existence of the global right action of the structure group picks out a special class of right invariant geometric objects on P—those which do not change when they are pulled back along R_g for all values of g in G. The most important right invariant objects on a principal bundle are the right invariant vector fields, which form an ideal ${\displaystyle {\mathfrak {E}}}$  of the Lie algebra of infinitesimal diffeomorphisms on P. Those vector fields on P which are both right invariant and vertical form an ideal ${\displaystyle V{\mathfrak {E}}}$  of ${\displaystyle {\mathfrak {E}}}$ , which has a relationship to the entire bundle P analogous to that of the Lie algebra ${\displaystyle {\mathfrak {g}}}$  of the gauge group G to the individual G-torsor fiber F.

The "field theory" of interest is defined in terms of a set of "fields" (smooth maps into various vector spaces) defined on a principal gauge bundle P. Different fields carry different representations of the gauge group G, and perhaps of other symmetry groups of the manifold such as the Poincaré group. One may define the space ${\displaystyle Pl}$  of local polynomials in these fields and their derivatives. The fundamental Lagrangian density of one's theory is presumed to lie in the subspace ${\displaystyle Pl_{0}}$  of polynomials which are real-valued and invariant under any unbroken non-gauge symmetry groups. It is also presumed to be invariant not only under the left action (passive coordinate transformations) and the global right action of the gauge group but also under local gauge transformations—pullback along the infinitesimal diffeomorphism associated with an arbitrary choice of right-invariant vertical vector field ${\displaystyle \epsilon \in V{\mathfrak {E}}}$ .

Identifying local gauge transformations with a particular subspace of vector fields on the manifold P provides a better framework for dealing with infinite-dimensional infinitesimals: differential geometry and the exterior calculus. The change in a scalar field under pullback along an infinitesimal automorphism is captured in the Lie derivative, and the notion of retaining only the term linear in the vector field is implemented by separating it into the interior derivative and the exterior derivative. In this context, "forms" and the exterior calculus refer exclusively to degrees of freedom which are dual to vector fields on the gauge bundle, not to degrees of freedom expressed in (Greek) tensor indices on the base manifold or (Roman) matrix indices on the gauge algebra.

The Lie derivative on a manifold is a globally well-defined operation in a way that the partial derivative is not. The proper generalization of Clairaut's theorem to the non-trivial manifold structure of P is given by the Lie bracket of vector fields and the nilpotence of the exterior derivative. This provides an essential tool for computation: the generalized Stokes theorem, which allows integration by parts and then elimination of the surface term, as long as the integrand drops off rapidly enough in directions where there is an open boundary. (This is not a trivial assumption, but can be dealt with by renormalization techniques such as dimensional regularization as long as the surface term can be made gauge invariant.)

Central to the BRST formalism is the BRST operator ${\displaystyle s_{B}}$ , defined as the tangent to the Ward operator ${\displaystyle W(\delta \lambda )}$ . The Ward operator on each field may be identified (up to a sign convention) with the Lie derivative along the vertical vector field associated with the local gauge transformation ${\displaystyle \delta \lambda }$  appearing as a parameter of the Ward operator. The BRST operator ${\displaystyle s_{B}}$  on fields resembles the exterior derivative on the gauge bundle, or rather to its restriction to a reduced space of alternating forms which are defined only on vertical vector fields. The Ward and BRST operators are related (up to a phase convention introduced by Kugo and Ojima, whose notation we will follow in the treatment of state vectors below) by ${\displaystyle W(\delta \lambda )X=\delta \lambda \;s_{B}X}$ . Here, ${\displaystyle X\in {Pl}_{0}}$  is a zero-form (scalar). The space ${\displaystyle {Pl}_{0}}$  is the space of real-valued polynomials in the fields and their derivatives that are invariant under any (unbroken) non-gauge symmetry groups.

Like the exterior derivative, the BRST operator is nilpotent of degree 2, i. e., ${\displaystyle (s_{B})^{2}=0}$ . The variation of any "BRST-exact form" ${\displaystyle s_{B}X}$  with respect to a local gauge transformation ${\displaystyle \delta \lambda }$  is given by the interior derivative ${\displaystyle \iota _{\delta \lambda }.}$  It is

${\displaystyle {\begin{aligned}\left[\iota _{\delta \lambda },s_{B}\right]s_{B}X&=\iota _{\delta \lambda }(s_{B}s_{B}X)+s_{B}\left(\iota _{\delta \lambda }(s_{B}X)\right)\\&=s_{B}\left(\iota _{\delta \lambda }(s_{B}X)\right)\end{aligned}}}$ 

Note that this is also exact.

The Hamiltonian perturbative formalism is carried out not on the fiber bundle, but on a local section. In this formalism, adding a BRST-exact term to a gauge invariant Lagrangian density preserves the relation ${\displaystyle s_{B}X=0.}$  This implies that there is a related operator ${\displaystyle Q_{B}}$  on the state space for which ${\displaystyle [Q_{B},{\mathcal {H}}]=0.}$  That is, the BRST operator on Fock states is a conserved charge of the Hamiltonian system. This implies that the time evolution operator in a Dyson series calculation will not evolve a field configuration obeying ${\displaystyle Q_{B}|\Psi _{i}\rangle =0}$  into a later configuration with ${\displaystyle Q_{B}|\Psi _{f}\rangle \neq 0}$  (or vice versa).

The nilpotence of the BRST operator can be understood as saying that its image (the space of BRST exact forms) lies entirely within its kernel (the space of BRST closed forms). The "true" Lagrangian, presumed to be invariant under local gauge transformations, is in the kernel of the BRST operator but not in its image. This implies that the universe of initial and final conditions can be limited to asymptotic "states" or field configurations at timelike infinity, where the interaction Lagrangian is "turned off". These states lie in the kernel of ${\displaystyle Q_{B},}$  but as the construction is invariant, the scattering matrix remains unitary. BRST-closed and exact states are defined similarly to BRST-closed and exact fields; closed states are annihilated by ${\displaystyle Q_{B},}$  while exact states are those obtainable by applying ${\displaystyle Q_{B}}$  to some arbitrary field configuration.

When defining the asymptotic states, the states that lie inside the image of ${\displaystyle Q_{B}}$  can also be suppressed, but the reasoning is a bit subtler. Having postulated that the "true" Lagrangian of the theory is gauge invariant, the true "states" of the Hamiltonian system are equivalence classes under local gauge transformation; in other words, two initial or final states in the Hamiltonian picture that differ only by a BRST-exact state are physically equivalent. However, the use of a BRST-exact gauge breaking prescription does not guarantee that the interaction Hamiltonian will preserve any particular subspace of closed field configurations that are orthogonal to the space of exact configurations. This is a crucial point, often mishandled in QFT textbooks. There is no a priori inner product on field configurations built into the action principle; such an inner product is constructed as part of the Hamiltonian perturbative apparatus.

The quantization prescription in the interaction picture is to build a vector space of BRST-closed configurations at a particular time, such that this can be converted into a Fock space of intermediate states suitable for Hamiltonian perturbation. As is conventional for second quantization, the Fock space is provided with ladder operators for the energy-momentum eigenconfigurations (particles) of each field, complete with appropriate (anti-)commutation rules, as well as a positive semi-definite inner product. The inner product is required to be singular exclusively along directions that correspond to BRST-exact eigenstates of the unperturbed Hamiltonian. This ensures that any pair of BRST-closed Fock states can be freely chosen out of the two equivalence classes of asymptotic field configurations corresponding to particular initial and final eigenstates of the (unbroken) free-field Hamiltonian.

The desired quantization prescriptions provide a quotient Fock space isomorphic to the BRST cohomology, in which each BRST-closed equivalence class of intermediate states (differing only by an exact state) is represented by exactly one state that contains no quanta of the BRST-exact fields. This is the appropriate Fock space for the asymptotic states of the theory. The singularity of the inner product along BRST-exact degrees of freedom ensures that the physical scattering matrix contains only physical fields. This is in contrast to the (naive, gauge-fixed) Lagrangian dynamics, in which unphysical particles are propagated to the asymptotic states. By working in the cohomology, each asymptotic state is guaranteed to have one (and only one) corresponding physical state (free of ghosts).

The operator ${\displaystyle Q_{B}}$  is Hermitian and non-zero, yet its square is zero. This implies that the Fock space of all states prior to the cohomological reduction has an indefinite norm, and so is not a Hilbert space. This requires that a Krein space for the BRST-closed intermediate Fock states, with the time reversal operator playing the role of the "fundamental symmetry" relating the Lorentz-invariant and positive semi-definite inner products. The asymptotic state space is then the Hilbert space obtained by quotienting BRST-exact states out of the Krein space.

To summarize: no field introduced as part of a BRST gauge fixing procedure will appear in asymptotic states of the gauge-fixed theory. However, this does not imply that these "unphysical" fields are absent in the intermediate states of a perturbative calculation! This is because perturbative calculations are done in the interaction picture. They implicitly involve initial and final states of the non-interaction Hamiltonian ${\displaystyle {\mathcal {H}}_{0}}$ , gradually transformed into states of the full Hamiltonian in accordance with the adiabatic theorem by "turning on" the interaction Hamiltonian (the gauge coupling). The expansion of the Dyson series in terms of Feynman diagrams will include vertices that couple "physical" particles (those that can appear in asymptotic states of the free Hamiltonian) to "unphysical" particles (states of fields that live outside the kernel of ${\displaystyle s_{B}}$  or inside the image of ${\displaystyle s_{B}}$ ) and vertices that couple "unphysical" particles to one another.

T. Kugo and I. Ojima are commonly credited with the discovery of the principal QCD color confinement criterion. Their role in obtaining a correct version of the BRST formalism in the Lagrangian framework seems to be less widely appreciated. It is enlightening to inspect their variant of the BRST transformation, which emphasizes the hermitian properties of the newly introduced fields, before proceeding from an entirely geometrical angle.

The ${\displaystyle {\mathfrak {g}}}$ -valued gauge fixing conditions are taken to be ${\displaystyle G=\xi \partial ^{\mu }A_{\mu },}$  where ${\displaystyle \xi }$  is a positive number determining the gauge. There are other possible gauge fixings, but are outside of the present scope. The fields appearing in the Lagrangian are:

• The QCD color field, that is, the ${\displaystyle {\mathfrak {g}}}$ -valued connection form ${\displaystyle A_{\mu }.}$ 
• The Faddeev–Popov ghost ${\displaystyle c^{i}}$ , which is a ${\displaystyle {\mathfrak {g}}}$ -valued scalar field with fermionic statistics.
• The antighost ${\displaystyle b_{i}={\bar {c}}_{i}}$ , also a ${\displaystyle {\mathfrak {g}}}$ -valued scalar field with fermionic statistics.
• The auxilliary field ${\displaystyle B_{i}}$  which is a ${\displaystyle {\mathfrak {g}}}$ -valued scalar field with bosonic statistics.

The field ${\displaystyle c}$  is used to deal with the gauge transformations, wheareas ${\displaystyle b}$  and ${\displaystyle B}$  deal with the gauge fixings. There actually are some subtleties associated with the gauge fixing due to Gribov ambiguities but they will not be covered here.

The BRST Lagrangian density is

${\displaystyle {\mathcal {L}}={\mathcal {L}}_{\textrm {matter}}(\psi ,\,A_{\mu }^{a})-{1 \over 4g^{2}}\operatorname {Tr} [F^{\mu \nu }F_{\mu \nu }]+{1 \over 2g^{2}}\operatorname {Tr} [BB]-{1 \over g^{2}}\operatorname {Tr} [BG]-{\xi \over g^{2}}\operatorname {Tr} [\partial ^{\mu }bD_{\mu }c]}$ 

Here, ${\displaystyle D_{\mu }}$  is the covariant derivative with respect to the gauge field (connection) ${\displaystyle A_{\mu }.}$  The Faddeev–Popov ghost field ${\displaystyle c}$  has a geometrical interpretation as a version of the Maurer–Cartan form on ${\displaystyle V{\mathfrak {E}}}$ , which relates each right-invariant vertical vector field ${\displaystyle \delta \lambda \in V{\mathfrak {E}}}$  to its representation (up to a phase) as a ${\displaystyle {\mathfrak {g}}}$ -valued field. This field must enter into the formulas for infinitesimal gauge transformations on objects (such as fermions ${\displaystyle \psi }$ , gauge bosons ${\displaystyle A_{\mu }}$ , and the ghost ${\displaystyle c}$  itself) which carry a non-trivial representation of the gauge group.

While the Lagrangian density isn't BRST invariant, its integral over all of spacetime, the action is. The transformation of the fields under an infinitessimal gauge transformation ${\displaystyle \delta \lambda }$  is given by

${\displaystyle {\begin{aligned}\delta \psi _{i}&=\delta \lambda D_{i}c\\\delta A_{\mu }&=\delta \lambda D_{\mu }c\\\delta c&=\delta \lambda {\tfrac {i}{2}}[c,c]\\\delta b=\delta {\bar {c}}&=\delta \lambda B\\\delta B&=0\end{aligned}}}$ 

Note that ${\displaystyle [\cdot ,\cdot ]}$  is the Lie bracket, NOT the commutator. These may be written in an equivalent form, using the charge operator ${\displaystyle Q_{B}}$  instead of ${\displaystyle \delta \lambda }$ . The BRST charge operator ${\displaystyle Q_{B}}$  is defined as

${\displaystyle Q_{B}=c^{i}\left(L_{i}-{\frac {1}{2}}{{f_{i}}^{j}}_{k}b_{j}c^{k}\right)}$ 

where ${\displaystyle L_{i}}$  are the infinitesimal generators of the Lie group, and ${\displaystyle f_{ij}{}^{k}}$  are its structure constants. Using this, the transformation is given as

${\displaystyle {\begin{aligned}Q_{B}A_{\mu }&=D_{\mu }c\\Q_{B}c&={i \over 2}[c,c]\\Q_{B}b&=B\\Q_{B}B&=0\end{aligned}}}$ 

The details of the matter sector ${\displaystyle \psi }$  are unspecified, as is left the form of the Ward operator on it; these are unimportant so long as the representation of the gauge algebra on the matter fields is consistent with their coupling to ${\displaystyle \delta A_{\mu }}$ . The properties of the other fields are fundamentally analytical rather than geometric. The bias is towards connections with ${\displaystyle \partial ^{\mu }A_{\mu }=0}$  is gauge-dependent and has no particular geometrical significance. The anti-ghost ${\displaystyle b={\bar {c}}}$  is nothing but a Lagrange multiplier for the gauge fixing term, and the properties of the scalar field ${\displaystyle B}$  are entirely dictated by the relationship ${\displaystyle \delta {\bar {c}}=i\delta \lambda B}$ . These fields are all Hermitian in Kugo–Ojima conventions, but the parameter ${\displaystyle \delta \lambda }$  is an anti-Hermitian "anti-commuting c-number". This results in some unnecessary awkwardness with regard to phases and passing infinitesimal parameters through operators; this can be resolved with a change of conventions.

We already know, from the relation of the BRST operator to the exterior derivative and the Faddeev–Popov ghost to the Maurer–Cartan form, that the ghost ${\displaystyle c}$  corresponds (up to a phase) to a ${\displaystyle {\mathfrak {g}}}$ -valued 1-form on ${\displaystyle V{\mathfrak {E}}}$ . In order for integration of a term like ${\displaystyle -i(\partial ^{\mu }{\bar {c}})D_{\mu }c}$  to be meaningful, the anti-ghost ${\displaystyle {\bar {c}}}$  must carry representations of these two Lie algebras—the vertical ideal ${\displaystyle V{\mathfrak {E}}}$  and the gauge algebra ${\displaystyle {\mathfrak {g}}}$ —dual to those carried by the ghost. In geometric terms, ${\displaystyle {\bar {c}}}$  must be fiberwise dual to ${\displaystyle {\mathfrak {g}}}$  and one rank short of being a top form on ${\displaystyle V{\mathfrak {E}}}$ . Likewise, the auxiliary field ${\displaystyle B}$  must carry the same representation of ${\displaystyle {\mathfrak {g}}}$  (up to a phase) as ${\displaystyle {\bar {c}}}$ , as well as the representation of ${\displaystyle V{\mathfrak {E}}}$  dual to its trivial representation on ${\displaystyle A_{\mu }.}$  That is, ${\displaystyle B}$  is a fiberwise ${\displaystyle {\mathfrak {g}}}$ -dual top form on ${\displaystyle V{\mathfrak {E}}}$ .

The one-particle states of the theory are discussed in the adiabatically decoupled limit g → 0. There are two kinds of quanta in the Fock space of the gauge-fixed Hamiltonian that lie entirely outside the kernel of the BRST operator: those of the Faddeev–Popov anti-ghost ${\displaystyle {\bar {c}}}$  and the forward polarized gauge boson. This is because no combination of fields containing ${\displaystyle {\bar {c}}}$  is annihilated by ${\displaystyle s_{B}}$  and the Lagrangian has a gauge breaking term that is equal, up to a divergence, to

${\displaystyle s_{B}\left({\bar {c}}\left(i\partial ^{\mu }A_{\mu }-{\tfrac {1}{2}}\xi s_{B}{\bar {c}}\right)\right).}$ 

Likewise, there are two kinds of quanta that will lie entirely in the image of the BRST operator: those of the Faddeev–Popov ghost ${\displaystyle c}$  and the scalar field ${\displaystyle B}$ , which is "eaten" by completing the square in the functional integral to become the backward polarized gauge boson. These are the four types of "unphysical" quanta which do not appear in the asymptotic states of a perturbative calculation.

The anti-ghost is taken to be a Lorentz scalar for the sake of Poincaré invariance in ${\displaystyle -i(\partial ^{\mu }{\bar {c}})D_{\mu }c}$ . However, its (anti-)commutation law relative to ${\displaystyle c}$  i.e. its quantization prescription, which ignores the spin–statistics theorem by giving Fermi–Dirac statistics to a spin-0 particle—will be given by the requirement that the inner product on our Fock space of asymptotic states be singular along directions corresponding to the raising and lowering operators of some combination of non-BRST-closed and BRST-exact fields. This last statement is the key to "BRST quantization", as opposed to mere "BRST symmetry" or "BRST transformation".

(Needs to be completed in the language of BRST cohomology, with reference to the Kugo–Ojima treatment of asymptotic Fock space.)

While the BRST symmetry, with its corresponding charge Q, elegantly captures the essence of gauge invariance, it presents a challenge for path integral quantization. The naive path integral, summing over all gauge configurations, vastly overcounts physically distinct states due to the redundancy introduced by gauge transformations. This overcounting manifests as a divergence in the path integral arising from integrating over the gauge orbits. To address this, we introduce a gauge-fixing procedure within the BRST framework.

The core idea is to restrict the path integral to a representative set of gauge configurations, eliminating the redundant gauge degrees of freedom. This is achieved by introducing a gauge-fixing function, denoted f(A), where A represents the gauge field. The specific choice of f(A) determines the gauge. Different choices lead to different representations of the same physical theory, though the final physical results must be independent of this choice.

The gauge-fixing procedure within BRST quantization is implemented by adding a term to the Lagrangian density that depends on both the gauge-fixing function and the ghost fields. This term is constructed to be BRST-exact, meaning it can be written as the BRST variation of some quantity. This ensures that the modified action still possesses BRST symmetry.

A general form for the gauge-fixing Lagrangian density is:

${\displaystyle L_{gf}=-iQ(f(A)*{\bar {c}})}$ 

where ${\displaystyle {\bar {c}}}$  is the antighost field. The factor of -i is a convention. Since Q² = 0, the BRST variation of L_gf is zero, preserving the BRST invariance of the total action.

Let's illustrate this with two common examples:

1. Gupta-Bleuler (Lorenz) Gauge in Electromagnetism:

In this gauge, the gauge-fixing function is ${\displaystyle f(A)=\partial _{\mu }A^{\mu }}$ . The gauge-fixing Lagrangian density becomes:

${\displaystyle L_{gf}=-iQ((\partial _{\mu }A^{\mu })*{\bar {c}})=-i((\partial _{\mu }\partial ^{\mu }c)*{\bar {c}}-(\partial _{\mu }A^{\mu })*B)}$ 

where B is an auxiliary Nakanishi-Lautrup field introduced to rewrite the gauge condition. After integrating out B in the path integral, we obtain the familiar form:

${\displaystyle L_{gf}=-{\frac {1}{2\xi }}(\partial _{\mu }A^{\mu })^{2}+(\partial _{\mu }{\bar {c}})(\partial ^{\mu }c)}$ 

where ξ is a gauge parameter. The Lorenz gauge corresponds to the Feynman gauge (ξ = 1). Note that the ghost fields remain coupled to the gauge field through the BRST variation.

2. ξ-Gauges in Yang-Mills Theories:

For non-Abelian gauge theories, a generalized class of ξ-gauges can be defined with the gauge-fixing function ${\displaystyle f(A)=\partial _{\mu }A^{\mu a}+\xi B^{a}}$ , where a is the gauge group index. The gauge-fixing Lagrangian density then becomes:

${\displaystyle L_{gf}=-iQ((\partial _{\mu }A^{\mu a}+\xi B^{a})*{\bar {c}}^{a})=B^{a}(\partial _{\mu }A^{\mu a})+(\xi /2)B^{a}B^{a}+{\bar {c}}^{a}(\partial _{\mu }D^{\mu }c)^{a}}$ 

where D^μ is the covariant derivative. The auxiliary field B^a can be integrated out, resulting in:

${\displaystyle L_{gf}=-{\frac {1}{2\xi }}(\partial _{\mu }A^{\mu a})^{2}+{\bar {c}}^{a}(\partial _{\mu }D^{\mu }c)^{a}}$ 

Again, ξ is a gauge parameter, and different choices of ξ correspond to different gauges within this family.

The introduction of the gauge-fixing term L_gf modifies the action and consequently the path integral. Crucially, the BRST symmetry is preserved, ensuring that physical observables remain independent of the gauge choice. Furthermore, the gauge-fixing procedure breaks the original gauge symmetry of the classical action, making the path integral well-defined. The ghost fields, originally introduced to compensate for the unphysical degrees of freedom, now play a crucial role in maintaining the unitarity of the theory in the quantized version.

This section only applies to classical gauge theories. i.e. those that can be described with first class constraints. The more general formalism is described using the Batalin–Vilkovisky formalism.

The BRST construction ^[1] applies to a situation of a Hamiltonian action of a gauge group ${\displaystyle G}$  on a phase space ${\displaystyle M}$ . Let ${\displaystyle {\mathfrak {g}}}$  be the Lie algebra of ${\displaystyle G}$  and ${\displaystyle 0\in {\mathfrak {g}}^{*}}$  a regular value of the moment map ${\displaystyle \Phi :M\to {\mathfrak {g}}^{*}}$ . Let ${\displaystyle M_{0}=\Phi ^{-1}(0)}$ . Assume the ${\displaystyle G}$ -action on ${\displaystyle M_{0}}$  is free and proper, and consider the space ${\displaystyle {\tilde {M}}}$  of ${\displaystyle G}$ -orbits on ${\displaystyle M_{0}}$ .

The Hamiltonian mechanics of a gauge theory is described by ${\displaystyle r}$  first class constraints ${\displaystyle \Phi _{i}}$  acting upon a symplectic space ${\displaystyle M}$ . ${\displaystyle M_{0}}$  is the submanifold satisfying the first class constraints. The action of the gauge symmetry partitions ${\displaystyle M_{0}}$  into gauge orbits. The symplectic reduction is the quotient of ${\displaystyle M_{0}}$  by the gauge orbits.

According to algebraic geometry, the set of smooth functions over a space is a ring. The Koszul-Tate complex (the first class constraints aren't regular in general) describes the algebra associated with the symplectic reduction in terms of the algebra ${\displaystyle C^{\infty }(M)}$ .

First, using equations defining ${\displaystyle M_{0}}$  inside ${\displaystyle M}$ , construct a Koszul complex

${\displaystyle ...\to K^{1}(\Phi )\to C^{\infty }(M)\to 0}$ 

so that ${\displaystyle H^{0}(K(\Phi ))=C^{\infty }(M_{0})}$  and ${\displaystyle H^{p}(K(\Phi ))=0}$  for ${\displaystyle p\neq 0}$ .

Then, for the fibration ${\displaystyle M_{0}\to {\tilde {M}}}$  one considers the complex of vertical exterior forms ${\displaystyle (\Omega _{vert}^{\cdot }(M_{0}),d_{vert})}$ . Locally, ${\displaystyle \Omega _{vert}^{\cdot }(M_{0})}$  is isomorphic to ${\displaystyle \Lambda ^{\cdot }V^{*}\otimes C^{\infty }({\tilde {M}})}$ , where ${\displaystyle \Lambda ^{\cdot }V^{*}}$  is the exterior algebra of the dual of a vector space ${\displaystyle V}$ . Using the Koszul resolution defined earlier, one obtains a bigraded complex

${\displaystyle K^{i,j}=\Lambda ^{i}V^{*}\otimes \Lambda ^{j}V\otimes C^{\infty }(M).}$ 

Finally (and this is the most nontrivial step), a differential ${\displaystyle s_{B}}$  is defined on ${\displaystyle K=\oplus _{i,j}K^{i,j}}$  which lifts ${\displaystyle d_{vert}}$  to ${\displaystyle K}$  and such that ${\displaystyle (s_{B})^{2}=0}$  and

${\displaystyle H_{s_{B}}^{0}=C^{\infty }({\tilde {M}})}$ 

with respect to the grading by the ghost number : ${\displaystyle K^{n}=\oplus _{i-j=n}K^{i,j}}$ .

Thus, the BRST operator or BRST differential ${\displaystyle s_{B}}$  accomplishes on the level of functions what symplectic reduction does on the level of manifolds.

There are two antiderivations, ${\displaystyle \delta }$  and ${\displaystyle d}$  which anticommute with each other. The BRST antiderivation ${\displaystyle s_{B}}$  is given by ${\displaystyle \delta +d+\mathrm {more} }$ . The operator ${\displaystyle s_{B}}$  is nilpotent; ${\displaystyle s^{2}=(\delta +d)^{2}=\delta ^{2}+d^{2}+(\delta d+d\delta )=0}$ 

Consider the supercommutative algebra generated by ${\displaystyle C^{\infty }(M)}$  and Grassman odd generators ${\displaystyle {\mathcal {P}}_{i}}$ , i.e. the tensor product of a Grassman algebra and ${\displaystyle C^{\infty }(M)}$ . There is a unique antiderivation ${\displaystyle \delta }$  satisfying ${\displaystyle \delta {\mathcal {P}}_{i}=-\Phi _{i}}$  and ${\displaystyle \delta f=0}$  for all ${\displaystyle f\in C^{\infty }(M)}$ . The zeroth homology is given by ${\displaystyle C^{\infty }(M_{0})}$ .

A longitudinal vector field on ${\displaystyle M_{0}}$  is a vector field over ${\displaystyle M_{0}}$  which is tangent everywhere to the gauge orbits. The Lie bracket of two longitudinal vector fields is itself another longitudinal vector field. Longitudinal ${\displaystyle p}$ -forms are dual to the exterior algebra of ${\displaystyle p}$ -vectors. ${\displaystyle d}$  is essentially the longitudinal exterior derivative defined by

${\displaystyle {\begin{aligned}d\omega (V_{0},\ldots ,V_{k})=&\sum _{i}(-1)^{i}d_{{}_{V_{i}}}(\omega (V_{0},\ldots ,{\widehat {V}}_{i},\ldots ,V_{k}))\\&+\sum _{i<j}(-1)^{i+j}\omega ([V_{i},V_{j}],V_{0},\ldots ,{\widehat {V}}_{i},\ldots ,{\widehat {V}}_{j},\ldots ,V_{k})\end{aligned}}}$ 

The zeroth cohomology of the longitudinal exterior derivative is the algebra of gauge invariant functions.

The BRST construction applies when one has a Hamiltonian action of a compact, connected Lie group ${\displaystyle G}$  on a phase space ${\displaystyle M}$ .^[2][3] Let ${\displaystyle {\mathfrak {g}}}$  be the Lie algebra of ${\displaystyle G}$  (via the Lie group–Lie algebra correspondence) and ${\displaystyle 0\in {\mathfrak {g}}^{*}}$  (the dual of ${\displaystyle {\mathfrak {g}})}$  a regular value of the momentum map ${\displaystyle \Phi :M\to {\mathfrak {g}}^{*}}$ . Let ${\displaystyle M_{0}=\Phi ^{-1}(0)}$ . Assume the ${\displaystyle G}$ -action on ${\displaystyle M_{0}}$  is free and proper, and consider the space ${\displaystyle {\widetilde {M}}=M_{0}/G}$  of ${\displaystyle G}$ -orbits on ${\displaystyle M_{0}}$ , which is also known as a symplectic reduction quotient ${\displaystyle {\widetilde {M}}=M/\!\!/G}$ .

First, using the regular sequence of functions defining ${\displaystyle M_{0}}$  inside ${\displaystyle M}$ , construct a Koszul complex

${\displaystyle \Lambda ^{\bullet }{\mathfrak {g}}\otimes C^{\infty }(M).}$ 

The differential, ${\displaystyle \delta }$ , on this complex is an odd ${\displaystyle C^{\infty }(M)}$ -linear derivation (differential algebra) of the graded ${\displaystyle C^{\infty }(M)}$ -algebra ${\displaystyle \Lambda ^{\bullet }{\mathfrak {g}}\otimes C^{\infty }(M)}$ . This odd derivation is defined by extending the Lie algebra homomorphism ${\displaystyle {\mathfrak {g}}\to C^{\infty }(M)}$  of the Hamiltonian action. The resulting Koszul complex is the Koszul complex of the ${\displaystyle S({\mathfrak {g}})}$ -module ${\displaystyle C^{\infty }(M)}$ , where ${\displaystyle S({\mathfrak {g}})}$  is the symmetric algebra of ${\displaystyle {\mathfrak {g}}}$ , and the module structure comes from a ring homomorphism ${\displaystyle S({\mathfrak {g}})\to C^{\infty }(M)}$  induced by the Hamiltonian action ${\displaystyle {\mathfrak {g}}\to C^{\infty }(M)}$ .

This Koszul complex is a resolution of the ${\displaystyle S({\mathfrak {g}})}$ -module ${\displaystyle C^{\infty }(M_{0})}$ , that is,

${\displaystyle H^{j}(\Lambda ^{\bullet }{\mathfrak {g}}\otimes C^{\infty }(M),\delta )={\begin{cases}C^{\infty }(M_{0})&j=0\\0&j\neq 0\end{cases}}}$ 

Then, consider the Chevalley–Eilenberg complex for the Koszul complex ${\displaystyle \Lambda ^{\bullet }{\mathfrak {g}}\otimes C^{\infty }(M)}$  considered as a differential graded module over the Lie algebra ${\displaystyle {\mathfrak {g}}}$ :

${\displaystyle K^{\bullet ,\bullet }=C^{\bullet }\left({\mathfrak {g}},\Lambda ^{\bullet }{\mathfrak {g}}\otimes C^{\infty }(M)\right)=\Lambda ^{\bullet }{\mathfrak {g}}^{*}\otimes \Lambda ^{\bullet }{\mathfrak {g}}\otimes C^{\infty }(M).}$ 

The "horizontal" differential ${\displaystyle d:K^{i,\bullet }\to K^{i+1,\bullet }}$  is defined on the coefficients

${\displaystyle \Lambda ^{\bullet }{\mathfrak {g}}\otimes C^{\infty }(M)}$ 

by the action of ${\displaystyle {\mathfrak {g}}}$  and on ${\displaystyle \Lambda ^{\bullet }{\mathfrak {g}}^{*}}$  as the exterior derivative of right-invariant differential forms on the Lie group ${\displaystyle G}$ , whose Lie algebra is ${\displaystyle {\mathfrak {g}}}$ .

Let Tot(K) be a complex such that

${\displaystyle \operatorname {Tot} (K)^{n}=\bigoplus \nolimits _{i-j=n}K^{i,j}}$ 

with a differential D = d + δ. The cohomology groups of (Tot(K), D) are computed using a spectral sequence associated to the double complex ${\displaystyle (K^{\bullet ,\bullet },d,\delta )}$ .

The first term of the spectral sequence computes the cohomology of the "vertical" differential ${\displaystyle \delta }$ :

${\displaystyle E_{1}^{i,j}=H^{j}(K^{i,\bullet },\delta )=\Lambda ^{i}{\mathfrak {g}}^{*}\otimes C^{\infty }(M_{0})}$ , if j = 0 and zero otherwise.

The first term of the spectral sequence may be interpreted as the complex of vertical differential forms

${\displaystyle (\Omega ^{\bullet }{\operatorname {vert} }(M_{0}),d_{\operatorname {vert} })}$ 

for the fiber bundle ${\displaystyle M_{0}\to {\widetilde {M}}}$ .

The second term of the spectral sequence computes the cohomology of the "horizontal" differential ${\displaystyle d}$  on ${\displaystyle E_{1}^{\bullet ,\bullet }}$ :

${\displaystyle E_{2}^{i,j}\cong H^{i}(E_{1}^{\bullet ,j},d)=C^{\infty }(M_{0})^{g}=C^{\infty }({\widetilde {M}})}$ , if ${\displaystyle i=j=0}$  and zero otherwise.

The spectral sequence collapses at the second term, so that ${\displaystyle E_{\infty }^{i,j}=E_{2}^{i,j}}$ , which is concentrated in degree zero.

Therefore,

${\displaystyle H^{p}(\operatorname {Tot} (K),D)=C^{\infty }(M_{0})^{g}=C^{\infty }({\widetilde {M}})}$ , if p = 0 and 0 otherwise.

• Batalin–Vilkovisky formalism
• Quantum chromodynamics

• Chapter 16 of Peskin & Schroeder (ISBN 0-201-50397-2 or ISBN 0-201-50934-2) applies the "BRST symmetry" to reason about anomaly cancellation in the Faddeev–Popov Lagrangian. This is a good start for QFT non-experts, although the connections to geometry are omitted and the treatment of asymptotic Fock space is only a sketch.
• Chapter 12 of M. Göckeler and T. Schücker (ISBN 0-521-37821-4 or ISBN 0-521-32960-4) discusses the relationship between the BRST formalism and the geometry of gauge bundles. It is substantially similar to Schücker's 1987 paper.^[1]

• Figueroa-O'Farrill, J. M.; Kimura, T. (1991). "Geometric BRST Quantization I. Prequantization". Commun. Math. Phys. 136 (2). Springer-Verlag: 209–229. Bibcode:1991CMaPh.136..209F. doi:10.1007/BF02100022. ISSN 0010-3616. MR 1096113. S2CID 120119621.
• Kostant, B.; Sternberg, S. (1987). "Symplectic reduction, BRS cohomology, and infinite-dimensional Clifford algebras". Ann. Phys. 176 (1). Elsevier: 49–113. Bibcode:1987AnPhy.176...49K. doi:10.1016/0003-4916(87)90178-3.

Original BRST papers:

• Brandt, Friedemann; Barnich, Glenn; Henneaux, Marc (2000), "Local BRST cohomology in gauge theories", Physics Reports, 338 (5): 439–569, arXiv:hep-th/0002245, Bibcode:2000PhR...338..439B, doi:10.1016/S0370-1573(00)00049-1, ISSN 0370-1573, MR 1792979, S2CID 119420167
• Becchi, C.; Rouet, A.; Stora, R. (1974). "The abelian Higgs Kibble model, unitarity of the S-operator". Physics Letters B. 52 (3). Elsevier BV: 344–346. Bibcode:1974PhLB...52..344B. doi:10.1016/0370-2693(74)90058-6. ISSN 0370-2693.
• Becchi, C.; Rouet, A.; Stora, R. (1975). "Renormalization of the abelian Higgs-Kibble model". Communications in Mathematical Physics. 42 (2). Springer Science and Business Media LLC: 127–162. Bibcode:1975CMaPh..42..127B. doi:10.1007/bf01614158. ISSN 0010-3616. S2CID 120552882.
• Becchi, C; Rouet, A; Stora, R (1976). "Renormalization of gauge theories". Annals of Physics. 98 (2). Elsevier BV: 287–321. Bibcode:1976AnPhy..98..287B. doi:10.1016/0003-4916(76)90156-1. ISSN 0003-4916.
• I.V. Tyutin, "Gauge Invariance in Field Theory and Statistical Physics in Operator Formalism", Lebedev Physics Institute preprint 39 (1975), arXiv:0812.0580.
• Kugo, Taichiro; Ojima, Izumi (1979). "Local Covariant Operator Formalism of Non-Abelian Gauge Theories and Quark Confinement Problem". Progress of Theoretical Physics Supplement. 66. Oxford University Press (OUP): 1–130. Bibcode:1979PThPS..66....1K. doi:10.1143/ptps.66.1. ISSN 0375-9687.
• A more accessible version of Kugo–Ojima is available online in a series of papers, starting with: Kugo, T.; Ojima, I. (1978-12-01). "Manifestly Covariant Canonical Formulation of the Yang-Mills Field Theories. I: -- General Formalism --". Progress of Theoretical Physics. 60 (6). Oxford University Press (OUP): 1869–1889. Bibcode:1978PThPh..60.1869K. doi:10.1143/ptp.60.1869. ISSN 0033-068X. This is probably the single best reference for BRST quantization in quantum mechanical (as opposed to geometrical) language.
• Much insight about the relationship between topological invariants and the BRST operator may be found in: E. Witten, "Topological quantum field theory", Commun. Math. Phys. 117, 3 (1988), pp. 353–386

• BRST systems are briefly analyzed from an operator theory perspective in: S. S. Horuzhy and A. V. Voronin, "Remarks on Mathematical Structure of BRST Theories", Comm. Math. Phys. 123, 4 (1989) pp. 677–685
• A measure-theoretic perspective on the BRST method may be found in Carlo Becchi's 1996 lecture notes.

• Brst cohomology on arxiv.org