---
abstract: 'To localize the scalar, fermion, and abelian gauge fields on our 3-brane, a simple mechanism with a hypothetical “magnetic field” in the bulk is proposed. This mechanism is to treat all fields in the equal footing without ad hoc consideration. In addition, the machanism can be easily realized in a flat dimension six Minkowski space and it works even in the weak coupling limit.'
author:
- 'Shih-Hao Ho[^1]'
- 'C. Q. Geng[^2]'
title: Localizing fields on brane in magnetized backgound 
---

It is known that the brane-world scenario allows us to use infinite extra dimensions without compactification. Since the low energy effective theory should be four-dimensional, the mechanism to localize matter fields naturally on a brane still is a main topic in the brane-world scenario [@Shifman:2009df]. In Ref. [@Rubakov:1983bb], Rubakov and Shaposhinikov proposed to localize fields on our 3-brane by topological defects[^3]. The mechanism to localize scalar and fermion is essentially given in Refs. [@Dashen:1974cj; @Jackiw:1975fn]. The feature of the mechanism is the existence of zero modes, identified to be our standard model fields. In the literature, there have been many proposals to trap fields in our four-dimesional world [@Shifman:2009df]. In this paper, we try to illustrate a simple and old way to localize particles on a plane by applying a strong magnetic field.

Let us consider the Lagrangian in a theory with a massive complex scalar field $\phi$, and a (classical) background gauge field $a_M$ associated with a $U(1)_X$ in a flat dimension six spacetime,[^4] \[scalar lagrangian\] \_S=(\_M+i g a\_M)\^\* (\^M-ig a\^M)- m\_s\^2\^\*. The equation of motion for these fields can be easily derived to be \[eom s\] &&=0.

For the case in which “magnetic field” is present along the extra dimensions, i.e. $f_{45}=-f_{54}\equiv B$, $f_{\mu\nu}=0$ for $\mu$,$\nu = 0,1,2$ and 3 in the gauge $a_M=(0,0,0,0,-B x_5/2,B x_4/2)$ with $f_{MN}\equiv \partial_M a_N - \partial_N a_M$,[^5] Eq.(\[eom s\]) becomes [^6] \[eom s1\] =0. Clearly, it reduces to a non-relativistic quantum mechanics problem. Using the raising and lowering operator method, we obtain the energy spectrum and ground state wave function: && p\_0\^2=p\_1\^2+p\_2\^2+p\_3\^2+(2n+1)|gB|+m\_s\^2  , n=0,1,2,\
&& \[eom s01\] \_[n=0]{}= N\_s e\^[ip\_x\^]{} r\^l e\^[i l]{} e\^[-|gB|r\^2/4]{}  , l=0,1,2,, representing the Landau levels for spin-zero particles with the lowest Landau level (LLL) for n=0 of the ground state energy $E_0^2=p_1^2+p_2^2+p_3^2+|gB|+m_s^2$, where $N_s$ is the normalization constant, $r^2=x_4^2+x_5^2$ and $tan \theta =x_5/x_4$. This shows that the ground state with respect to the bulk dynamics is localized near $r^2 \equiv r_0^2 = 2 l/|gB|$ as long as the bulk magnetic field is large enough. However, we observe that, in the presence of the term $|gB|$, these “zero modes” [^7] carry information from extra dimensions even if the bulk mass is absent.

We emphasis that the “center” of wavefunction in Eq.(\[eom s01\]) is $r=r_0$, which $depends$ on the angular quantum number $l$. Hence the field is concentrated at $r=r_0(l)$ with characteristic length $1/|gB|$ for a fixed $l$[^8], we cannot argue which value of $l$ is a preferred one since all values of $l$ correspond to the LLL. The absolute position in extra dimensions, however, is irrelevant in our consideration, the relevant quantity is always the relative positions of different fields along the extra dimensions. We believe this will become much clearer when the interactions between fields are taken into account. For present purpose, we emphasis the $dynamics$ of fields is effectively confined on 3-brane even if we consider the free fields in an magnetized background. This can be seen as we consider the Green’s function in the second part of this paper.

We now turn to the fermion case, the goal is still to find a zero mode localized on our brane [@Jackiw:1984ji]. The Lagrangian for a massive Dirac fermion $\Psi$ on a background gauge field $a_M$ in a flat dimension six spacetime is \[dirac lagrangian\] \_D=|(i D\_[M]{}\^[M]{}-m\_f), where $D_M\equiv \partial_M-ig a_M$ with the coupling constant g and the representations of 8 by 8 Dirac matrices are chosen to be: \[gamma\] \^=\^ \_1 = (

[cc]{}0 & \^\
\^ & 0

), &&   \^4=i\^5 \_1=(

[cc]{}0 & i \^[5]{}\
 i \^[5]{} & 0

),\
\^5=<span style="font-variant:small-caps;">1</span>i \_2=(

[cc]{}0 & 1\
 -1 & 0

), &&    \^7=\^0\^1\^2\^3\^4\^5= (

[cc]{}1 &0\
0& -1

), with the four dimensional gamma matrices in the chiral representation. The gamma matrices in Eq. (\[gamma\]) satisfy the Clifford algebra {\^M ,\^N} &= & 2 \^[MN]{},  \^[MN]{}=diag(+1,-1,-1,-1,-1,-1),\
{\^[M]{}, \^7 } &=&0,  (\^7)\^2=1, with M,N =$0,1,2,3,4,5$. By defining the 6D chiral states for the eight-componet spinor $\Psi$: = (

[c]{}\_+\
\_-

), with  \^7\_=\_, we can decompose equations of motion for $\Psi_{\pm}$ into && (i\^D\_-\^5 D\_4+i **1** D\_5)\_-=m\_f \_+ ,\
&& (i\^D\_-\^5 D\_4-i **1** D\_5)\_+=m\_f \_- . Consequently, these equations can be rewritten as && (i\^D\_-\^5 D\_4-i **1** D\_5)(i\^D\_-\^5 D\_4+i **1** D\_5)\_-=m\_f\^2 \_-,\
&& (i\^D\_-\^5 D\_4+i **1** D\_5)(i\^D\_-\^5 D\_4-i **1** D\_5)\_+=m\_f\^2 \_+ , which lead to &&\_-=0 \[eom f1\],\
&&\_+=0 .\[eom f2\] Since $\gamma^5$ commutes with the whole operator in the brackets in Eqs. (\[eom f1\]) and (\[eom f2\]), it can be replaced by the eigenvalue $\gamma$ which is the chirality of a 4D fermion: \_=(

[c]{}\_[\_L]{}\
\_[\_R]{}

),  \^5 \_[\_[L,R]{}]{}=\_[\_[L,R]{}]{}, with $\gamma=-1 (+1)$ for L(R). We see that both Eqs. (\[eom f1\]) and (\[eom f2\]) take the same form as Eq. (\[eom s1\]). Therefore, we can write down the energy spectrum and ground state wavefunction immediately: && p\_0\^2=p\_1\^2+p\_2\^2+p\_3\^2+(2n+1+sgn(gB))|gB|+m\_f\^2  , n=0,1,2,   for  \_- \[eom f3\]\
&& p\_0\^2=p\_1\^2+p\_2\^2+p\_3\^2+(2n+1-sgn(gB))|gB|+m\_f\^2  , n=0,1,2,   for   \_+ \[eom f4\]\
&&\_[\_[n=0]{}]{}= N\_f e\^[ip\_x\^]{} r\^l e\^[i l]{} e\^[-|gB|r\^2/4]{}  \_  , l=0,1,2,3.... where $N_f$ is the normalization constant, $\omega_{\gamma}$ denotes the eigenvector of $\gamma^5$, $r^2=x_4^2+x_5^2$ and $tan \theta =x_5/x_4$ as mentioned early. Now, we observe that for the lowest energy state with $n=0$ and $\gamma=\pm sgn(gB)$ for $\Psi_\pm$, the energy is $E_0^2=p_1^2+p_2^2+p_3^2$ if the bulk mass is absent. Clearly, there exist two zero modes with respect to the bulk dynamics, corresponding to $\Psi_{-L}$ and $\Psi_{+R}$ if $sgn(gB)>0$.

Next, we consider the simplest case of complex massive vector gauge fields $A_M$ associated with $U(1)_A$, which also carry a charge of $U(1)_X$, i.e. $A_M \rightarrow A_M+\partial_M \Lambda(x)$ and $A_M \rightarrow e^{i \alpha(x)} A_M$ under $U(1)_A$ and $U(1)_X$ gauge transformations, respectively. Based on Ref. [@Jackiw:1997vw], we can construct a gauge invariant Lagrangian under these two $U(1)$ groups with the help of a scalar Stueckelberg field $\rho$ whose transformation laws are $\rho \rightarrow \rho+ m_v \Lambda^{\prime}(x)$ and $\rho \rightarrow e^{i \alpha(x)} \rho$ under $U(1)_A$ and $U(1)_X$, respectively [@Ruegg:2003ps]. Then, the Langrangian reads: \[v\] \_V=\^[\*]{}\_[MN]{}\^[MN]{}-ig f\^[NM]{}V\^\*\_N V\_M+m\_v\^2 V\^\*\_M V\^M-|D\_M A\^M+ m\_v|\^2, where $V_M \equiv A_M-\frac{1}{m}D_M \rho$, $\mathbb{F}_{MN} \equiv D_M A_N-D_N A_M$ and $\xi$ is the gauge parameter. The gauge invariance of $U(1)_X$ is obvious and the invariance under $U(1)_A$ is guaranteed by the nontrivial condition of $D_M \Lambda^{\prime} =\partial_M \Lambda$.[^9] We can derive equations of motion by varying $A_M$ and $\rho$, given by (D\^2 \^[MN]{} && -D\^M D\^N)A\_M +D\_M( f\^[MN]{})- i f\^[NM]{}A\_M+f\^[NM]{}D\_M\
&&+m\_v\^2(A\^N-D\^N ) +(D\^N D\^M A\_M+m\_v D\^N )=0 \[v1\] and &&D\_N ( D\_M \^[MN]{}- f\^[NM]{}(A\_M -D\_M )+m\_v(A\^N-D\^N ))-m\_v\^2=0 \[v2\] . We will concentrate on the case with $f^{MN}=$constant and choose the Feynman gauge $\xi=1$. Consequently, Eq.(\[v1\]) is reduced to[^10] A(x) \_[M]{}=0\[v3\]. Here, we still have degree of freedoms to set $A_0=0$ and write $A_M \equiv A(x) \epsilon_M$ with the polarization vector $\epsilon_M$ \_=(

[c]{}\_[1]{}\
\_[2]{}\
\_[3]{}\
0\
0

),   =1,2,3;   \_= (

[c]{}0\
0\
0\
1\
i

). In this basis, Eq. (\[v3\]) becomes diagonal, given by A(x) \_[M]{}=0 \[v4\], where $\chi=0$ for $M=1,2,3$ and $\chi=\mp 1$ for $M=\pm$. Hence, the energy spectrum and coordinate space component of $A_M$ is the solution of Eq.(\[eom s1\]) with the replacement $m^2 \rightarrow m^2+2gB\chi$, found to be && p\_0\^2=p\_1\^2+p\_2\^2+p\_3\^2+(2n+1+2 sgn(gB))|gB|+m\_v\^2  , n=0,1,2,\
&& A(x)\_[n=0]{}= N\_v e\^[ip\_x\^]{} r\^l e\^[i l]{} e\^[-|gB|r\^2/4]{}  , l=0,1,2,, where $N_v$ is the normalization constant, $r^2=x_4^2+x_5^2$ and $tan \theta =x_5/x_4$. The ground state energy is $E_0^2=p_1^2+p_2^2+p_3^2+(1+2 sgn(gB)\chi)|gB|+m_v^2$. Therefore, in the viewpoint of 4D, we have a localized massive vector field $A_{\mu}$ with an effective mass $\sqrt{|gB|}$ in 4D and two localized scalar field $A_{\pm}$ which carry an extra degree of freedom with respect to the extra dimensional “handness” with effective squared masses of $m_v^2+(1+2sgn(gB))|gB|$ and $m_v^2+(1-2sgn(gB))|gB|$ for “right-handed” and “left-handed” states, respectively.

There are three possible cases of $A_{\pm}$, which are listed in Table I.

\[t\]

        Mass parameter                     Effective masses of zero modes                          Type of modes
  --------------------------- -------------------------------------------------------- -------------------------------------
         $m_v^2>|gB|$          $m_{eff}^2|_{\chi=\pm 1}=m_v^2+(1\mp 2sgn(gB))|gB|>0$        $A_{\pm}$ are both massive
   $m_v^2=|gB|$, $sgn(gB)>0$   $m_{eff}^2|_{\chi=+1}=4|gB|$, $m_{eff}^2|_{\chi=-1}=0$      $A_+$ massless, $A_-$ massive
   $m_v^2=|gB|$, $sgn(gB)<0$   $m_{eff}^2|_{\chi=-1}=4|gB|$, $m_{eff}^2|_{\chi=+1}=0$      $A_-$ massless, $A_+$ massive
   $m_v^2<|gB|$, $sgn(gB)>0$     $m_{eff}^2|_{\chi=+1}>0$, $m_{eff}^2|_{\chi=-1}<0$     $A_+$ tachyonic mode, $A_-$ massive
   $m_v^2<|gB|$, $sgn(gB)<0$     $m_{eff}^2|_{\chi=-1}>0$, $m_{eff}^2|_{\chi=+1}<0$     $A_-$ tachyonic mode, $A_+$ massive

  : Types of zero modes for $A_{\pm}$

\[T1\]

We note that in the third case, there is a heavy scalar for $\chi=-sgn(gB)$ and a unstable mode for $\chi=sgn(gB)$. This is the so-called zero mode problem and a general feature when we consider the spectrum of a vector field in the presence of magnetic field [@Nielsen:1978rm; @Hashimoto:1997gm].[^11] We also note that as a special case, to have a massless vector field $A_{\mu}$ in 4D, we may re-consider the Lagrangian in Eq. (\[v\]) with an “imaginary” mass term as \_V=\^[\*]{}\_[MN]{}\^[MN]{}-ig f\^[NM]{}V\^\*\_N V\_M-|D\_M A\^M+m\_v|\^2+(i m\_v\^)\^2 V\^\*\_M V\^M, where $(m^{\prime}_v)^2=|gB|$.

Despite of the localized wavefunctions, we emphasize that the effective dynamics is confined on the 3-brane when the strong background field (and massless limit) is applied. This effective dimensional reduction can be seen in either the coordinate/momentum space reduction [@Gusynin:1999pq] or the phase space reduction [@Dunne:1989hv]. The phase space reduction is examined in detail in Ref. [@Dunne:1989hv]. Here, we would like to see the propagators for scalar and fermion fields and show that their propagations are indeed localized on the brane with the characteristic length $l_B^{-2} \equiv |gB|$.

Using Schwinger’s proper time method [@Schwinger:1951nm], we can calculate the propagators of scalar and fermion fields in the presence of a constant background. Since the presence of the constant background field breaks the translational invariance, we recover the translation invariance by factoring out the Schwinger phase: &&G\_s(x,y)=e\^[i g \^x\_y a\_M(\_M) d\^M ]{} \_s(x-y) ,\
&&G\_f(x,y)=e\^[i g \^x\_y a\_M(\_M) d\^M ]{} \_f(x-y), where $G_s$ and $G_f$ stand for the propagators for scalars and fermions, respectively, and the path of the integral is a straight line.

The translational parts of the propagators $\tilde{G}_s$ and $\tilde{G}_f$ in the momentum space are given by \[sp\] (p)&=&d\^6x (x) e\^[ipx]{} =\^\_0 ds e\^[-is(m\^2+p\_\^2+(p\_4\^2+p\_5\^2))]{} ,\
\[fp\] (p)&=&d\^6x (x) d\^[ipx]{} =\^\_0 ds e\^[-is(m\^2+p\_\^2+(p\_4\^2+p\_5\^2))]{}\
&& \[(p\_ \^+m)(1+sgn(gB)\^4\^5 tan(|gB|s)) -(p\_4\^4+p\_5\^5)(1+tan\^2(|gB|s)\],which can be decomposed to be [@Gusynin:1999pq; @Chodos:1990vv]: &&(p)=i e\^[-]{} \_[n=0]{}\^(-1)\^n , \[sp1\]\
&&(p)=i e\^[-]{} \_[n=0]{}\^, \[fp1\] respectively, with D\_n(gB,p)=&&(p\_\^+m)\
&&+4(p\_4\^4+p\_5\^5) L\^1\_[n+1]{}(x), where $p_{\bot}^2\equiv p_4^2+p_5^2$, $x \equiv \frac{2(p_4^2+p_5^2)}{|gB|}$ and $L_n^{(\alpha)}(x)$ are the (associated) Laguerre polynomials. In the LLL for $n=0$, the propagators take the forms: &&(p)=i e\^[-]{} \[sp2\],\
&&(p)= i e\^[- ]{} (1-i\_4\_5) \[fp2\], which are the exact free propagators in (3+1) dimension in the region of $p_{\bot}^2<<|gB|$. We can also transform them back to the coordinate space: &&\_s(x-y)=e\^[-]{} g\_s(x\_-y\_),\
&&\_f(x-y)=e\^[-]{} g\_f(x\_-y\_)(1-i\_4\_5), where $g_s(x_{\mu})$ and $g_f(x_{\mu})$ are the ordinary scalar and fermion propagators in (3+1) spacetime, given by && g\_s(x\_-y\_)= e\^[-i p\_(x-y)\^]{},\
&& g\_f(x\_-y\_)= e\^[-i p\_(x-y)\^]{}. It is clear that the propagations along the extra dimension are indeed suppressed by the Gaussian function with the characteristic length $l_B$.

To consider in Feynman gauge $\xi=1$, we start from the differential equation G\^[NP]{}\_v(x,y)=i \^P\_M \^[(6)]{}(x-y) \[v5\]. By applying a rotation matrix $\mathbf{P}$ with ()\_M\^N=(

[cccccc]{}1 & 0 & 0 & 0 & 0 & 0\
0 & 1 & 0 & 0 & 0 & 0\
0 & 0 & 1 & 0 & 0 & 0\
0 & 0 & 0 & 1 & 0 & 0\
0 & 0 & 0 & 0 & &\
0 & 0 & 0 & 0 & &

), Eq.(\[v5\]) in the matrix form becomes i \^[(6)]{} (x-y)&&=\^[-1]{} \^[-1]{}\
&&= , where $\mathbf{\hat{G}}\equiv \mathbf{P}^{-1} \mathbf{G}\mathbf{P} $ and $\chi=0,-1,1$ for $M=0\sim3, 4,5$, respectively. It is clear that we can obtain the scalar propagator again with the matrix structure $g_{MN}$ and the replacement $m_s^2\rightarrow m_v^2+2gB\chi$. The translational invariant part of the propagator for the vector field (up to a Schiwinger phase) is given by: \_[MN]{}(p)=\_M\^P\_[PQ]{}(p)()\^Q\_N=\_M\^P ()\^Q\_N , where $G_s(p)$ is given in Eq. (\[sp\]).

Some remarks are in order. First, we emphasize that we can take the limit $g\rightarrow 0$ but keep the combination $|gB|>>p^2_{\bot}$ to guarantee the use of the LLL. Secondly, as shown explicitly in Eqs. (\[sp2\]) and (\[fp2\]), when we consider a scattering amplitude, the momentum integral is effectively four-dimensional since the transverse momentum integral is factored out and becomes a Gaussian-type integral which is only a finite constant. This is still true even if we go beyond the LLL as we can see from Eqs. (\[sp1\]) and (\[fp1\]).

In summary, we have considered a simple way to localize scalar, fermion and abelian gauge fields on our 3-brane. The advantage of this consideration is that it is a simple and well-studied setup. We note that the remnant from the “magnetic field” gives effective mass terms for scalars even if we set the the bulk mass to zero, leading to an effective massive scalar theory in (3+1) spacetime. On the other hand, if we start from a massless bulk fermion theory, the effective theory in (3+1) spacetime will still be massless.

To localize a massive gauge field, an additional Stuckelberg field is needed to preserve the gauge invariance. The choice of this bulk mass term is crucial for us to avoid the “zero-mode” problem. However, if we would like to reduce to an effective massless abelian gauge theory in (3+1) spacetime, an imaginary mass could be considered, leading to a tachyonic mode in the theory. This unavoidable tachyonic mode shall have interesting phenomenology and deserve further studies.

Acknowledgments {#acknowledgments .unnumbered}
===============

We are grateful to Professor Roman Jackiw for many useful discussions and reading the manuscript. We would like to thank Yuhsin Tsai for helpful discussions. The work by S.H.H. is supported by the National Science Council of R.O.C. under Grant number: NSC98-2917-I-564-122 and that by C.Q.G. is supported in part by the National Science Council of R.O.C. under Grant number: NSC-98-2112-M-007-008-MY3.

[99]{}

For a review, see M. Shifman, Int. J. Mod. Phys.  A [**25**]{}, 199 (2010). V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett.  B [**125**]{}, 136 (1983). G. A. Palma, Phys. Rev.  D [**73**]{}, 045023 (2006) \[arXiv:hep-th/0505170\].

R. F. Dashen, B. Hasslacher and A. Neveu, Phys. Rev.  D [**10**]{}, 4130 (1974). R. Jackiw and C. Rebbi, Phys. Rev.  D [**13**]{}, 3398 (1976).

T.  Haugset, J.  Aa.  Ruud and F.  Ravndal, Phys. Scr. [**47**]{}, 715 (1993).

R. Jackiw, Phys. Rev.  D [**29**]{}, 2375 (1984) \[Erratum-ibid.  D [**33**]{}, 2500 (1986)\] R. Jackiw, Phys. Rev.  D [**57**]{}, 2635 (1998). For a review on the Stuckelberg formalism, see H. Ruegg and M. Ruiz-Altaba, Int. J. Mod. Phys.  A [**19**]{}, 3265 (2004). N. K. Nielsen and P. Olesen, Nucl. Phys.  B [**144**]{}, 376 (1978). A. Hashimoto and W. Taylor, Nucl. Phys.  B [**503**]{}, 193 (1997).

This dimensional reduction effect has been widely studied in (3+1) and (2+1) dimensions. See, for example, V. P. Gusynin, V. A. Miransky and I. A. Shovkovy, Nucl. Phys.  B [**563**]{}, 361 (1999). and reference therein.

G. V. Dunne, R. Jackiw and C. A. Trugenberger, Phys. Rev.  D [**41**]{}, 661 (1990); G. V. Dunne and R. Jackiw, Nucl. Phys. Proc. Suppl.  [**33C**]{}, 114 (1993); R. Jackiw, Nucl. Phys. Proc. Suppl.  [**108**]{}, 30 (2002) \[Phys. Part. Nucl.  [**33**]{}, S6 (2002)\] \[Lect. Notes Phys.  [**616**]{}, 294 (2003)\]. J. S. Schwinger, Phys. Rev.  [**82**]{}, 664 (1951). A. Chodos, K. Everding and D. A. Owen, Phys. Rev.  D [**42**]{}, 2881 (1990).

[^1]: E-mail address: shho@mit.edu

[^2]: E-mail address: geng@phys.nthu.edu.tw

[^3]: In Ref.  [@Palma:2005xv], the author attempted to reproduce SM phenomenology exploiting the type of fermion local- ization by the mechanism introduced in Ref.  [@Rubakov:1983bb]

[^4]: Our convention is $\eta^{MN}=diag(+1,-1,-1,-1,-1,-1)$, while Greek and Latin indices represent the (3+1) spacetime and all the dimensions indices, respectively.

[^5]: The wavefunctions $depend$ on the gauge choice and they differ only by a phase factor as expected [@Haugset].

[^6]: Here we note the location of this 3-brane along the extra dimensions is $arbitrary$ and it also depends on the gauge choice of classical field $a_M$. For instance, the choice of $a_M(x)= (0,0,0,0,-B(x_5+a)/2, B(x_4+b)/2)$ with some constants $a$, $b$ will give us the same background “magnetic field" but it will also change the location of the brane from $x_4=x_5=0$ to $x_4=-b$, $x_5=-a$. But the important property is that the wavefunctions and the propagators fall off exponentially outside the brane.

[^7]: The “zero mode" referred by us means “zero mode of Landau level", i.e. the Lowest Landau Level (LLL).

[^8]: This can be seen more clearly in Landau gauge.

[^9]: Note that $(D_M D^M+\xi m^2)\Lambda^{\prime}=0$.

[^10]: In the original Stuckelberg formalism, the vector field $A_M$ and scalar field $\rho$ satisfy the same field equation [@Ruegg:2003ps]: &&(\_M \^M+m\^2)A\_M =0\
    &&(\_M \^M+m\^2) =0 In our consideration, the field equations are not the same owing to an external background field. We can derive the equation of motion for scalar field $\rho$ from Eq. (\[v2\]) by using the relation $[D_M , D_N]= - i g f_{MN}$: (D\_[M]{}D\^[M]{}+m\^2)=0.

[^11]: In Ref. [@Hashimoto:1997gm], a gauge theory with U(2) on a $T^2$ has been considered. There exists an instability in this guage theory, which corresponds to the instability in a brane-anti-brane configuration.
