Hamiltonian, Path Integral and BRST Formulations of Large N Scalar $QCD_{2}$ on the Light-Front and Spontaneous Symmetry Breaking

Recently Grinstein, Jora, and Polosa have studied a theory of large-$N$ scalar quantum chromodynamics in one-space one-time dimension. This theory admits a Bethe-Salpeter equation describing the discrete spectrum of quark-antiquark bound states. They consider gauge fields in the adjoint representation of $SU(N)$ and scalar fields in the fundamental representation. The theory is asymptotically free and linearly confining. The theory could possibly provide a good field theoretic framework for the description of a large class of diquark-antidiquark (tetra-quark) states. Recently we have studied the light-front quantization of this theory without a Higgs potential. In the present work, we study the light-front Hamiltonian, path integral and BRST formulations of the theory in the presence of a Higgs potential. The light-front theory is seen to be gauge-invariant, possessing a set of first-class constraints. The explicit occurrence of spontaneous symmetry breaking in the theory is shown in unitary gauge as well as in the light-front 't Hooft gauge.

Various possibilities for understanding hadron structure beyond the usual mesons and baryons [3,4] have been considered in the literature [1]- [26]. Some exotic states find a natural interpretation in terms of the four-quark or tetra-quark (qqqq) states [3], [9]- [26]. By now it is widely perceived that not only heavy states such as the X, Y, Z states have an exotic structure as tetra-quark states or diquark(Q)-antidiquark(Q) states [3], [9]- [19], but even some light scalar mesons could also be identified as diquark-antidiquark (QQ) or tetra-quark systems [20]- [25].
The work of Grinstein et al. [25] is seen to further support this hypothesis. In the work of Grinstein et al. [25], the gauge fields have been considered [25] in the adjoint representation of SU(N) and the scalar fields in the fundamental representation. The theory is asymptotically free and linearly confining [25]. Different aspects of this theory have been studied by several authors in various contexts [20]- [26].
Also, because there is no spin-statistics connection in one space and one-time dimension, the spinor QCD 2 is structurally similar to the scalar QCD 2 [5]- [8]. It is therefore enough to consider the scalar QCD 2 for a study of several aspects of QCD 2 [5]- [8]. The large-N behavior of scalar QCD 2 has been studied in details by 't Hooft and others [4]- [8].
In view of the above, the motivations for our present studies could be easily highlighted. In the first place, the work of 't Hooft, Isidori, Maiani, Polosa and Riquer [24] and others [20]- [25] has clearly shown as to how one could achieve a satisfactory explanation of light scalar meson decays by assuming a dominant diquark-antidiquark structure for the lightest scalar mesons [20]- [25] (where the diquark is being taken to be a spin zero antitriplet color state). In this work, a coherent picture of scalar mesons as a mixture of tetra-quark states (dominating the lightest mesons) and heavy quark-antiquark states (dominating the heavier mesons) emerges [24].
The studies of Grinstein, Jora and Polosa [25] on the large-N scalar QCD 2 [20]- [25], further support the hypothesis of 't Hooft, Isidori, Maiani, Polosa and Riquer [24] and others [20]- [25], about the assumption of a dominant diquark-antidiquark structure for the scalar mesons. The work of Grinstein, Jora and Polosa is based on the assumption that scalar QCD 2 with a large number of colors could be used to compute the mass spectrum as well as to estimate the mass of the first radial excitation of the lowest diquark-antidiquark scalar meson. They have applied a numerical procedure to solve the Bethe-Salpeter equations and compute the bound state discrete spectrum of this confining theory [25]. They have even obtained the possible masses of the spinor and scalar quarks by imposing that the ratio of the ground state eigenvalues of the spinor and scalar Bethe-Salpeter equations respectively, is equal to the ratio of the physical masses m π /m σ (cf. Ref. [25], for further details). They have even extended their discussion to the case of spin-one diquarks.
The above studies of Grinstein, Jora and Polosa [25], based on the scalar QCD with a large number of colors in one-space one-time dimension, clearly point towards some definite possibilities of gaining some insight, at least at the qualitative level, about the physical tetra-quark states in three-space one-time dimension. In addition to this, it may also be possible to study this theory, in three-space one-time dimension, at a somewhat later point of time.
In view of the above, it seems reasonable to pursue these studies further. In fact, in a recent paper [26], we have studied the light-front (LF) quantization (LFQ) [28]- [36] of this theory (with a mass term for the complex scalar (diquark) field but without the Higgs potential) on the hyperplanes defined by the equal light-cone time τ = x + = (x 0 + x 1 )/ √ 2 = constant [30]- [36], using the Hamiltonian [28] and path integral [29]- [31] formulations.
In the present work, we study the LF Hamiltonian, path integral, and Becchi-Rouet-Stora and Tyutin (BRST) [37]- [39] formulations of this theory in the presence of a Higgs potential [32]- [36] under appropriate light-cone gauge-fixing conditions. The LF theory is seen to be gauge-invariant (GI), possessing a set of first-class constraints. We absorb the mass term for the complex scalar (diquark) field φ in the definition of our Higgs potential [25], [26], and then we study the action of the theory.
One of the important motivations for introducing the Higgs potential is to study the aspects related to the spontaneous symmetry breaking (SSB) [32]- [34] in the theory. Another important motivation for introducing the Higgs potential in the theory is related to our long-term goal related to the study of this theory using the discrete light-cone (LC) quantization (DLCQ) along with the coherent state formalism [40]- [48], where we wish not only to study the aspects of the spontaneous symmetry breaking (SSB) but we also wish to make contact with the experimentally observational aspects of this theory using the LF Hamiltonian approach to study the two-and three-body relativistic bound state problems [25], [40]- [48]. This work therefore constitutes a part of our bigger project which involves a study of some aspects related to the spontaneous symmetry breaking as well as to a study of its DLCQ using the coherent state formalism [40]- [48], in the LF Hamiltonian approach to study the two-and three-body relativistic bound state problems [40]- [48].
In this sense one could think that the theory under our present consideration could perhaps provide a good basic field theoretic framework for a study of a large class of diquark-antidiquark or the tetra-quark states [2,3], [9]- [26] which have been investigated in various experiments. These are some of the motivations that necessitate our present studies. Now, because the theory is GI, we also study its BRST quantization [37]- [39] under appropriate BRST light-cone gauge-fixing. Usually in the Hamiltonian and the path integral quantization of a theory the gauge-invariance of the theory gets broken because the procedure of gauge-fixing converts the set of first-class constraints of the theory into a set of second-class constraints. A possible way to achieve the quantization of a GI theory, such that the gauge-invariance of the theory is maintained even under gauge-fixing is to use a generalized procedure called the BRST quantization [37]- [39], where the extended gauge symmetry, called the BRST symmetry, is maintained even under gauge-fixing.
In the next section, we briefly recap some basics of the instant-form (IF) quantization (IFQ) of this theory in the presence of a Higgs potential. Its LFQ in the presence of a Higgs potential is then considered in Sec. 3, using the Hamiltonian and path integral formulations. The light-front BRST formulation of the theory is studied in Sec. 4, under the appropriate BRST light-cone gauge-fixing. Finally, the summary and discussion are given in Sec. 5.

Some Basics of the Theory
In this section we recap some of the basics of this theory of large-N scalar QCD 2 in the presence of a Higgs potential, studied earlier by Grinstein, Jora and Polosa without a Higgs potential (but with a mass term for the complex scalar (diquark) field φ ) [25] (the mass term for the complex scalar (diquark) field φ is absorved in the definition of our Higgs potential [25,26]). The theory of large-N scalar QCD 2 that we propose to study is defined by the action [25]: Here α s is the QCD coupling constant. The covariant derivative in our considerations is defined as: where A µ (≡ A a µ T a ) are the gluon gauge fields and T a are the generators of Lie algebra corresponding to the group SU(N c ) obeying the commutation relations: with N c = 2 for SU (2) and N c = 3 for SU (3) . The structure constants f abc are antisymmetric in all indices. The gluon gauge field strength F a µν is defined as: Here (A µ × A ν ) a = f abc A b µ A c ν defines the cross product for any two "isotopic" vectors: A a µ and A a ν [8]. Further, the scalar fields φ and φ † transform as the N andN representations of the U(N) color group respectively [2]. Also, following the work of Grinstein, Jora and Polosa [25], we ignore all gluon self-coupling terms that arise from our chosen Lagrangian.
In the Lagrangian density of our theory (defined by Eq.(1)), the first term represents the kinetic energy of the gluon field, the second term represents the kinetic energy term for the scalar (diquark) field, the third term represents the interaction term for the scalar (diquark) field with the gluon field, and the last term represents the Higgs potential which is kept rather general, without making any specific choice for the parameters µ 2 and λ. However, they are chosen such that the potential remains a double well potential with the vacuum expectation value φ 0 =< 0|φ(x)|0 > = 0, so as to allow for the spontaneous symmetry breaking in the theory. Also, the mass term for the scalar (diquark) field has been absorbed in the definition of the Higgs potential. The values: µ 2 = m 2 and λ = 0 reproduce the theory of Grinstein, Jora and Polosa [25].

Instant-Form Quantization
We now consider the instant-form (IF) quantization (IFQ) of the theory. The action of the above theory in the IF of dynamics (with [25]: Here t = x 0 = x 0 and x = x 1 = −x 1 . Canonical momenta obtained from the above action are: Here π, π † , Π 0 (≡ Π 0a T a ) and E := Π 1 (≡ Π 1a T a ) are the momenta canonically conjugate respectively to φ, φ † , A 0 and A 1 . The above equations however, imply that the theory possesses only one primary constraint: The symbol ≈ here denotes a weak equality in the sense of Dirac [28], and it implies that the constraints hold as a strong equality only on the reduced hyper surface of the constraints and not in the rest of the phase space of the classical theory (and similarly one can consider it as a weak operator equality for the corresponding quantum theory). The canonical Hamiltonian density corresponding to L is: After including the primary constraint χ 1 in the canonical Hamiltonian density H c with the help of the Lagrange multiplier field u, the total Hamiltonian density H T could be written as : Hamilton's equations of motion of the theory that preserve the constraints of the theory in the course of time could be obtained from the total Hamiltonian: H T = H T dx 1 (and are omitted here for the sake of brevity). Demanding that the primary constraint χ 1 be preserved in the course of time, one obtains the secondary Gauss-law constraint of the theory as: The preservation of χ 2 for all times gives rise to one further constraint: The theory is thus seen to possess only three constraints χ i (with i = 1,2,3). The matrix R αβ of the Poisson brackets among the set of constraints χ i with (i = 1, 2, 3) is seen to be singular (the other details of the matrix R αβ are omitted here for the sake of brevity). This implies that the set of constraints χ i is first-class and that the theory under consideration is gauge-invariant (GI). Consequently the theory is seen to possess the local vector gauge symmetry defined by the local vector gauge transformations: where β ≡ β(x 0 , x 1 ) is an arbitrary real function of its arguments. This theory could now be quantized under some appropriate gauge-fixing conditions, e.g., under the time-axial or temporal gauge: A 0 ≈ 0. The details of this IFQ are however, outside the scope of the present work (what actually happens is that one of the matrix elements of the matrix R αβ involves a linear combination of a Dirac distribution function δ(x 1 − y 1 ) and its first derivative and finding its inverse is a rather non-trivial task). We now proceed with the LFQ of this theory in the next section.

Light-Front Hamiltonian and Path Integral Quantization
In this section we study the LF Hamiltonian and path integral formulations [25]- [31] of the above theory [25] under appropriate LC gauge-fixing. The action for the scalar theory In the work of Ref. [25], the authors have studied the above action, after implementing the gauge-fixing condition (GFC) A + ≈ 0 "strongly" in the above action. In contrast to this, we propose to study the theory defined by the above action, following the standard Dirac quantization procedure [28] and we do not fix any gauge at this stage. We instead consider this gauge-fixing condition (A + ≈ 0) as one of the gauge constraints [28]- [31] which becomes strongly equal to zero only on the reduced hyper surface of the constraints and remains non-zero in the rest of the phase space of the theory. It may be important to note here that one of the salient features of Dirac quantization procedure [28] is that in this quantization the gauge-fixing conditions should be treated on par with other gauge-constraints of the theory which are only weakly equal to zero in the sense of Dirac [28], and they become strongly equal to zero only on the reduced hyper surface of the constraints of the theory and not in the rest of the phase space of the classical theory (in the corresponding quantum theory these weak equalities become the weak operator equalities).
Another thing to be noted here is that we have introduced the Higgs potential in our present work and we have absorbed the mass term for the scalar (diquark) field in the definition of our Higgs potential [25,26], [32]- [34]. We now proceed to study the LF Hamiltonian and path integral formulations of the theory defined by the above action. The LF Euler-Lagrange equations of motion of the theory are: In the following, we consider the Hamiltonian formulation of the theory described by the above action. The canonical momenta obtained from the above action are: Here π, π † , Π + (≡ Π +a T a ) and Π − (≡ Π −a T a ) are the momenta canonically conjugate respectively to φ, φ † , A − and A + . The above equations however, imply that the theory possesses three primary constraints: (17) The canonical Hamiltonian density corresponding to L is: After including the primary constraints χ 1 , χ 2 and χ 3 in the canonical Hamiltonian density H c with the help of the Lagrange multiplier fields u, v and w, the total Hamiltonian density H T could be written as : The Hamilton's equations of motion of the theory that preserve the constraints of the theory in the course of time could be obtained from the total Hamiltonian (and are omitted here for the sake of brevity): H T = H T dx − . Demanding that the primary constraint χ 1 be preserved in the course of time, one obtains the secondary Gauss-law constraint of the theory as: The preservation of χ 2 , χ 3 and χ 4 , for all times does not give rise to any further constraints. The theory is thus seen to possess only four constraints χ i (with i = 1,2,3,4). The constraints χ 2 , χ 3 and χ 4 could however, be combined in to a single constraint: and with this modification, the new set of constraints of the theory could be written as: Further, the matrix of the Poisson brackets among the constraints Ω i , with (i = 1, 2) is seen to be a singular matrix implying that the set of constraints Ω i is first-class and that the theory under consideration is gauge-invariant. Expressions for the components of the vector gauge current density of the theory are obtained as: The divergence of the vector gauge current density of the theory could now be easily seen to vanish satisfying the continuity equation: ∂ µ j µ = 0 , implying that the theory possesses at the classical level, a local vector-gauge symmetry. The action of the theory is indeed seen to be invariant under the local vector gauge transformations: where β ≡ β(x + , x − ) is an arbitrary real function of its arguments and Π u , Π v and Π w are the momenta canonically conjugate to the Lagrange multiplier fields u, v and w respectively, which are treated here as dynamical fields. Using the Euler-Lagrange equations of motion of the theory and the expressions for the components of the vector gauge current density of the theory, one could now easily show that: It may be important to point out here that Grinstein, Jora and Polosa [25], have obtained an equation (under the gauge A + = 0) analogous to the above equation connecting ∂ − ∂ − A − and j + (cf. Eq. (5) of Ref. [25]) which has been shown [25], to admit a solution (in the absence of background fields) [25], which when substituted in to the Lagrangian density of the theory implies a linear potential between charges (for further details, we refer to the work of Ref. [25]). In order to quantize the theory using Dirac's procedure we now convert the set of first-class constraints of the theory η i into a set of second-class constraints, by imposing, arbitrarily, some additional constraints on the system called gauge-fixing conditions (GFC's) or the gauge-constraints [28]- [31]. For the present theory, we could choose, for example, the following set of GFC's: ζ 1 = A + ≈ 0 , ζ 2 = A − ≈ 0. Here the gauge A + ≈ 0 represents the LC time-axial or temporal gauge and the gauge A − ≈ 0 represents the LC coulomb gauge and both of these gauges are physically important gauges. Corresponding to this gauge choice, the theory has the following set of constraints under which the quantization of the theory could for example be studied: The matrix R αβ of the Poisson brackets among the set of constraints ξ i with (i = 1, 2, 3, 4) is seen to be nonsingular with the determinant given by The other details of the matrix R αβ are omitted here for the sake of brevity. Finally, following the Dirac quantization procedure, the nonvanishing equal light-cone-time commutators of the theory, under the GFC's: A + ≈ 0 and A − ≈ 0 are obtained as: The first-order Lagrangian density L I0 of the theory is: In the path integral formulation [29]- [31], the transition to quantum theory is made by writing the vacuum to vacuum transition amplitude for the theory called the generating functional Z[J k ]. For the present theory [25] under the GFC's: ζ 1 = A + ≈ 0 and ζ 2 = A − ≈ 0 and in the presence of the external sources J k it reads: Here, the phase space variables of the theory are: Φ k ≡ (φ, φ † , A − , A + , u, v, w) with the corresponding respective canonical conjugate momenta: Π k ≡ (π, π † , Π + , Π − , Π u , Π v , Π w ). The functional measure [dµ] of the generating functional Z[J k ] under the above gaugefixing is obtained as : The LF Hamiltonian and path integral quantization of the theory under the set of GFC's: A + ≈ 0 and A − ≈ 0 is now complete.

Spontaneous Symmetry Breaking
In this section, we consider the spontaneous symmetry breaking (SSB) in the theory in (i) the so-called unitary gauge and (ii) in the 't Hooft gauge and show explicitly the existence of SSB [32]- [34] in the theory in both the cases. Our Higgs potential possesses a local maximum at where the phase angle θ defines a direction in the complex φ− plane. Here the vacuum state (or the ground state) of the system is clearly non unique, and the SSB will occur for any particular choice of the value of θ. In our considerations we however, choose θ = 0 which in turn implies: where we choose v > 0. We now parametrize the field φ(x) in terms of its deviations from its vacuum expectation value (VEV): < 0| φ(x) |0 >= φ 0 = (v/ √ 2) > 0 in terms of two real fields σ(x)and η(x), which measure the deviations of the field φ(x) from the equilibrium ground state configuration φ(x) = φ 0 . For this we expand our complex scalar field φ(x), in terms of two real fields σ(x) and η(x) as: such that the real fields σ(x) and η(x) have vanishing vacuum expectation values. In fact, the term ϕ 0 here could be interpreted as the zero mode of the theory [32]- [34] and the fluctuation field ϕ(x) could be interpreted as the normal mode of the theory [32]- [34]. The Lagrangian density of our LF theory in terms of the real fields σ(x) and η(x) (after dropping the terms which are irrelevant for our discussions namely, a constant term, a term linear in the field σ(x) and all the quartic interaction terms in the fields) (with m σ = (2µ 2 + λv 2 )/2 and m v = |vρ| ) becomes: The first term in the above Lagrangian density represents the kinetic energy of the electromagnetic field; the second term represents the kinetic energy of the real scalar field σ(x); the third term represents the kinetic energy of the real scalar field η(x); the fourth term represents the mass term for the real scalar field σ(x); the fifth term, which involves the product of the fields A + (x) and A − (x) with the derivatives of the field η(x), represents a quadratic interaction term involving the fields A + (x), A − (x) and η(x), and it implies that the fields A + (x), A − (x) and η(x) are not independent normal coordinates and are therefore not free fields; consequently the sixth and seventh terms cannot be interpreted as the mass terms for the real scalar field η(x) and the electromagnetic field respectively. It also implies that the above Lagrangian density contains an unphysical field; we will eliminate it under some suitable gauge. The last four terms in the above Lagrangian density represent simply the cubic interaction terms of the theory which will be needed for our later discussions.

The Unitary Gauge and SSB
We now consider this theory in the so-called unitary gauge. In fact, for any complex field φ(x), a gauge transformation can be found which transforms φ(x) into a real field such as: The gauge in which the transformed field has this form is called as the unitary gauge. With this substitution the transformed Lagrangian density in the so-called unitary gauge (after dropping as before, the terms which are irrelevant for our discussions, namely, a constant term, a term linear in the field σ(x) and the quartic interaction terms) becomes: The interaction part of the above Lagrangian density however, does not contain any quadratic coupling terms involving the coupling of the fields σ(x), A + (x) and A − (x).
Hence treating the interaction part of the Lagrangian density L int U (x), in perturbation theory, one could interpret L 0 U as the free-field Lagrangian density of a real Klein-Gordon field σ(x) and a real massive vector field A µ (x). Upon quantizing the theory, the field σ(x) gives rise to neutral scalar bosons of mass m σ = (2µ 2 + λv 2 )/2 and the field A µ (x) gives rise to neutral vector bosons of mass m v = |vρ|. This is an explicit demonstration of the SSB in the theory through the Higgs mechanism where the massive spin 0 boson associated with the field σ(x) is the Higgs boson (or Higgs Scalar) of the theory. Here the vector field A µ (x) has become massive in the process of SSB through the Higgs mechanism.

The Light-Front 't Hooft Gauge and SSB
We consider the LF 't Hooft gauge defined by: and construct the LF 't Hooft-gauge-fixed Lagrangian density of the theoryL: by adding the LF 't Hooft-gauge-fixing term L′ tH : to the Lagrangian density of the theory expressed in terms of the real scalar fields σ and η given by Eq. (33) and ignoring the terms irrelevant for our discussion as explained in the forgoing. The LF 't Hooft gauge-fixed action of the theoryS could now be written after a partial integration (and with m η = (6µ 2 + λv 2 )/6 ) as: The fields σ(x), η(x) and A µ (x), could now be treated in perturbation theory as three independent fields which could be quantized in the usual manner. The LF 't Hooft gauge here reintroduces the field η(x) which gets eliminated in the so-called unitary gauge. However, there are no real particles corresponding to the quantized η(x) field and they appear in a manner akin to the longitudinal and scalar photons of QED theory.
which mixes operators with Bose and Fermi statistics. We then introduce new anticommuting variable c andc (Grassman numbers on the classical level and operators in the quantized theory) and a commuting variable b such that [37]- [39]: with the propertyδ 2 = 0. We now define a BRST-invariant function of the dynamical phase space variables of the theory to be a function f such thatδf = 0. Now the BRST gauge-fixed quantum Lagrangian density L BRST for the theory could be obtained by adding to the first-order Lagrangian density L I0 , a trivial BRST-invariant function, e.g. as follows: The last term in the above equation is the extra BRST-invariant gauge-fixing term. After one integration by parts, the above equation could now be written as: Proceeding classically, the Euler-Lagrange equation for b reads: the requirementδb = 0 then implies which in turn implies The above equation is also an Euler-Lagrange equation obtained by the variation of L BRST with respect toc. In introducing momenta one has to be careful in defining those for the fermionic variables. We thus define the bosonic momenta in the usual manner so that but for the fermionic momenta with directional derivatives we set implying that the variable canonically conjugate to c is ( ∂ +c ) and the variable conjugate toc is (∂ + c). For writing the Hamiltonian density from the Lagrangian density in the usual manner we remember that the former has to be Hermitian so that: The consistency of the last two equations could now be easily checked by looking at the Hamilton's equations for the fermionic variables. Also for the operators c,c, ∂ + c and ∂ +c , one needs to satisfy the anticommutation relations of ∂ + c withc or of ∂ +c with c, but not of c, withc. In general, c andc are independent canonical variables and one assumes that [37]- [39]: where { , } means an anticommutator. We thus see that the anticommulators in the above equation are non-trivial and need to be fixed. In order to fix these, we demand that c satisfy the Heisenberg equation: and using the property c 2 = c 2 = 0 one obtains The last three equations then imply : Here the minus sign in the above equation is nontrivial and implies the existence of states with negative norm in the space of state vectors of the theory. The BRST charge operator Q is the generator of the BRST transformations. It is nilpotent and satisfies Q 2 = 0. It mixes operators which satisfy Bose and Fermi statistics. According to its conventional definition, its commutators with Bose operators and its anti-commutators with Fermi operators for the present theory satisfy: All other commutators and anti-commutators involving Q vanish. In view of this, the BRST charge operator of the present theory can be written as: This equation implies that the set of states satisfying the conditions: belong to the dynamically stable subspace of states |ψ > satisfying Q|ψ >= 0, i.e., it belongs to the set of BRST-invariant states. In order to understand the condition needed for recovering the physical states of the theory we rewrite the operators c andc in terms of fermionic annihiliation and creation operators. For this purpose we consider Euler lagrange equation for the variable c derived earlier. The solution of this equation gives (for the light-cone time x + ≡ τ ) the Heisenberg operators c(τ ) and correspondinglyc(τ ) in terms of the fermionic annihilation and creation operators as: Which at the light-cone time τ = 0 imply By imposing the conditions (obtained earlier): we then obtain Now let |0 > denote the fermionic vacuum for which Defining |0 > to have norm one, the last three equations imply The theory is thus seen to possess negative norm states in the fermionic sector. The existence of these negative norm states as free states of the fermionic part of H BRST is however, irrelevant to the existence of physicsl states in the orthogonal subspace of the Hilbert space. In terms of annihilation and creation operators H BRST is: and the BRST charge operator is: Now because Q|ψ >= 0, the set of states annihiliated by Q contains not only the set for which the constraints of the theory hold but also additional states for which The Hamiltonian is also invariant under the anti-BRST transformation given by: with generator or anti-BRST chargē which in terms of annihilation and creation operators reads: We also have with and we further impose the dual condition that both Q andQ annihilate physical states, implying that: The states for which the constraints of the theory hold, satisfy both of these conditions and are in fact, the only states satisfying both of these conditions, since although with there are no states of this operator with G † |ψ >= 0 and F † |ψ >= 0, and hence no free eigenstates of the fermionic part of H BRST that are annihilated by each of G, G † , F , and F † . Thus the only states satisfying Q|ψ >= 0 andQ|ψ >= 0 are those that satisfy the constraints of the theory. Now because Q|ψ >= 0, the set of states annihilated by Q contains not only the set of states for which the constraints of the theory hold but also additional states for which the constraints of the theory do not hold. This situation is, however, easily avoided by additionally imposing on the theory, the dual condition: Q|ψ >= 0 andQ|ψ >= 0. By imposing both of these conditions on the theory simultaneously, one finds that the states for which the constraints of the theory hold are the only states satisfying both of these conditions. This is traced to the conditions on the fermionic variables c andc which constrain the solutions such that one cannot have simultaneously c, ∂ + c andc, ∂ +c , applied to |ψ > giving zero. Thus the only states satisfying Q|ψ >= 0 andQ|ψ >= 0 are those that satisfy the constraints of the theory and they belong to the set of BRST-invariant as well as to the set of anti-BRST-invariant states.
Alternatively, one can understand the above point in terms of fermionic annihiliation and creation operators as follows. The condition Q|ψ >= 0 implies the that the set of states annihilated by Q contains not only the states for which the constraints of the theory hold but also additional states for which the constraints do not hold. However,Q|ψ >= 0 guarantees that the set of states annihilated byQ contains only the states for which the constraints hold, simply because G † |ψ > = 0 and F † |ψ > = 0. This completes the BRST formulation of the theory.

Summary and Discussion
Theoretical and experimental studies of multi-quark states are challenging and a number of phenomenological models [1]- [26] have been proposed in order to provide interpretation and gain understanding.
In particular, as mentioned in the foregoing, Grinstein, Jora and Polosa [25] have studied a model of large-N scalar QCD 2 [25]. This theory of Grinstein et al. [25] admits a Bethe-Salpeter equation describing the discrete spectrum of qq bound states [20]- [25]. In the their work, the gauge fields have been considered [25] in the adjoint representation of SU(N) and the scalar fields in the fundamental representation. The theory is asymptotically free and linearly confining [25]. Different aspects of this theory have been studied by several authors in various contexts [20]- [25].
In Ref. [26], we have studied the LFQ of the theory of large-N scalar QCD 2 studied by Grinstein, Jora and Polosa [25], without Higgs potential [20]- [25] on the LF using the Hamiltonian [28] and path integral [29]- [31] formulations. In in the present work, we have studied this theory in the presence of a Higgs potential and we have studied its LFQ using the Hamiltonian, path integral and BRST formulations [37]- [39]. We have also shown explicitly the occurrence of the SSB in the theory in the unitary gauge as well as in the LF 't Hooft gauge [32]- [34].
In the Hamiltonian and path integral quantization of the theory the gauge-invariance of the theory gets broken because of the gauge-fixing. In view of this, we go to a more generalized quantization procedure called the BRST quantization [37]- [39], where the extended gauge symmetry of the theory is maintained even under gauge-fixing.
In the present work, we have studied the LF-BRST quantization of the theory under some specific LF-BRST gauge-fixing (where a particular but non-unique gauge has been chosen). In this procedure, we embed the original GI theory into a BRST system, the quantum Hamiltonian H BRST (which includes the gauge-fixing contribution) commutes with the BRST charge as well as with the anti-BRST charge. The new extended gauge symmetry which replaces the gauge invariance is maintained (even under the BRST gaugefixing) and projecting any state onto the sector of BRST and anti-BRST invariant states yields a theory which is isomorphic to the original GI theory.