Anomaly and Sign problem in $\mathcal{N}=(2,2)$ SYM on Polyhedra : Numerical Analysis
Abstract
We investigate the twodimensional supersymmetric YangMills (SYM) theory on the discretized curved space (polyhedra). We first revisit that the number of supersymmetries of the continuum SYM theory on any curved manifold can be enhanced at least to two by introducing an appropriate gauge background associated with the symmetry. We then show that the generalized Sugino model on the discretized curved space, which was proposed in our previous work, can be identified to the discretization of this SUSY enhanced theory, where one of the supersymmetries remains and the other is broken but restored in the continuum limit. We find that the anomaly exists also in the discretized theory as a result of an unbalance of the number of the fermions proportional to the Euler characteristics of the polyhedra. We then study this model by using the numerical MonteCarlo simulation. We propose a novel phasequench method called “anomalyphasequenched approximation” with respect to the anomaly. We numerically show that the WardTakahashi (WT) identity associated with the remaining supersymmetry is realized by adopting this approximation. We figure out the relation between the sign (phase) problem and pseudozeromodes of the Dirac operator. We also show that the divergent behavior of the scalar onepoint function gets milder as the genus of the background increases. These are the first numerical observations for the supersymmetric lattice model on the curved space with generic topologies.
Contents:
1 Introduction
One of the most important and nontrivial aspects of gauge field theories is the dynamics. Even if the theory is strongly restricted by some symmetry like supersymmetry, it is in general not sufficient to determine the whole dynamics. This leads to a strong motivation to pursue a way to examine the dynamical aspect of the supersymmetric gauge theories nonperturbatively. Among a great deal of approaches in this direction, the application of the lattice technique to supersymmetric gauge theories has been a longstanding theme in theoretical high energy physics [1, 2, 3, 4]. Although the lattice regularization breaks the Poincare invariance to its discrete subgroup and the supersymmetry cannot be straightforwardly realized on the lattice, several ways of bypassing the problem have been developed to date. In particular, for lowdimensional gauge theories with extended supersymmetries, lattice models could avoid finetuning for the continuum limit due to partially preserved supercharges. In [5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21], the lattice supersymmetric gauge models preserving one or two supercharges are proposed based on the discretized topologically twisted gauge theories (for relations among several lattice formulations, see [22, 23, 24, 25, 26]). Since the supersymmetries are partially preserved in the models, the numerical simulation can be carried out on the basis similar to lattice QCD [27, 28, 29, 30, 31, 32, 33, 34]. The problem of the vacuum degeneracy of lattice gauge fields is also shown to be avoided without using an admissibility condition [35]. For discretization of higher dimensional supersymmetric gauge theories without finetuning, see [36, 37]. Relevance of the lattice supersymmetric models to fermion discretizations is discussed in [38].
The lattice formulations of supersymmetric theories are in general discretized on a periodic hypercubic lattice, that is, the topology is restricted only on the torus. However, the topology plays an significant role in supersymmetric gauge theories especially in the context of topological field theory [39, 40, 41, 42, 43, 44, 45]. The significance of such theories has recently been rerecognized in relation to the localization technique in supersymmetric gauge theories [46].
In [47], a discretization of the topologically twisted twodimensional supersymmetric YangMills (SYM) theory on an arbitrary Riemann surface is proposed (generalized Sugino model). There, the Riemann surface is discretized to a polygon, on which the supersymmetric gauge theory preserving a single supercharge is defined. It is shown that the theory can be defined on any decomposition of the twodimensional surface and one can take its continuum limit without any finetuning. In [48], the analytical study based on the localization technique for the model is performed and it was shown that the partition function of the model mainly depends on the Euler characteristics of the background Riemann surface.
In this paper, we investigate the twodimensional SYM on a curved spacetime, where two supersymmetries survive under an appropriate gauge backgrounds associated with the gauged symmetry [49, 50], both from the theoretical and numerical viewpoints. What we will show is summarized in the followings:

We show that the generalized Sugino model is nothing but a discretization of the SUSY enhanced theory, where one of the supersymmetries is preserved and the other is broken but restored in the continuum limit. Both these continuum and discretized supersymmetric theories are not topological but physical in the sense that we can consider any kinds of operators as observables, that is, we do not need to restrict the observables to cohomology unlike the traditional topological field theories.

We show that the anomaly exists also in the generalized Sugino model as a result of the unbalance of the fermion numbers which is related to the Euler characteristics of background.

We study the generalized Sugino model based on the numerical MonteCarlo simulation. We show that the flat directions of the scalar fields can be controlled by adding a mass term. By investigating the behavior of the expectation values of the onepoint function of scalar fields, we find that the divergent behavior of the expectation value in a small mass region gets milder as the genus of background geometry increases.

We divide the phase of the Pfaffian of the Dirac operator into the anomalyinduced phase and the residual phase by introducing a specific operator called the ”compensator”. We approximate the theory by ignoring the residual phase in the Monte Carlo simulation, which we call the “anomalyphasequenched approximation”. We show that the WardTakahashi (WT) relations expected from the remaining supersymmetry of the model are satisfied based on this approximation.

We show that the sign problem in the model originates in the pseudozeromodes of the Dirac operators, which causes the anomaly. In other words, the Pfaffian phase except the anomalyinduced phase (the residual phase) does not contribute the path integral regardless of the topology of the background
^{5} . This fact guarantees the validity of the anomalyphase quenched approximation.
This paper is organized as follows: In Sec. 2, we discuss SUSY enhancement in the continuum SYM theory on the curved space with introducing an appropriate gauge backgrounds. We also discuss its relation to the generalized Sugino model on the discretized curved space. In Sec. 3, we discuss the anomaly in the discretized model and propose the anomalyphasequenched approximation. In Sec. 4, we show the results of Monte Carlo simulations for the scalar onepoint functions and the WT identity in the model. We also show the origin of the anomaly and the origin of the sign problem based on the numerical methods.
2 Continuum field theory and Generalized Sugino model
2.1 Rigid supersymmetry on a curved space
We first review a general aspect of the rigid supersymmetry on a curved manifold. The rigid supersymmetry on the curved manifold is developed to apply the localization to various field theories, whose general constructions are discussed in [49, 50] and more detailed discussion for twodimensional case is given in [51].
The rigid supersymmetry on the flat spacetime is the usual supersymmetry, that is, it is realized as a fermionic variation of the fields via
(2.1) 
where the index runs the number of the supercharges which we are considering, are the supercharges and are globally constant fermionic parameters on the whole flat spacetime. The supersymmetry of the theory is guaranteed by the invariance of the action and the measure under the variation by (2.1). In particular the Lagrangian of the theory is varied up to the total derivatives. However, if we consider the same kind of fermionic transformation on a curved spacetime, it is in general hard to make an action invariant under the transformation (2.1) since the fixed spinor does not commute with the covariant derivative in general even if it has a position dependence as long as it is a fixed function. Thus we sometimes make the supersymmetry local and build a supergravity by introducing extra fields and symmetries.
One exception is when the curved spacetime admits covariantly constant Killing spinors,
(2.2) 
In this case, we can make an invariant action under the transformation (2.1) with the parameter satisfying (2.2), since they commute with the covariant derivatives and act as a constant parameter in the variation of the Lagrangian. The existence of the covariantly constant Killing spinor restricts the structure of the manifolds. The Killing spinors exist only on Ricci flat Kähler manifolds for example. This is the reason why the rigid supersymmetry on the curved manifold has not been considered for a long time.
However, there is a route to avoid this obstacle. The point is that many supersymmetric theories have global symmetries called the Rsymmetries and the parameter is charged under them. So if we gauge one of the global Rsymmetries, the Killing equation (2.2) is modified to
(2.3) 
where is a vector field corresponding to the gauged Rsymmetry. The solutions of (2.3) include broader possibilities than (2.2) depending on a choice of . In particular, if the field strength from the gauge field cancels the effect of the curvature, we can obtain the covariantly constant Killing spinor on the curved manifold in the deformed sense. In this way, the rigid supersymmetry on the curved space is constructed by gauging the Rsymmetry. This construction of the rigid supersymmetry is also understood from a point of view of fixed background fields in the supergravity theory as discussed in [49, 50].
2.2 rigid supersymmetry on the Riemann surface
Now let us apply the above strategy to the 2D SYM theory. This theory on the twodimensional flat spacetime is obtained by the dimensional reduction from the 4D SYM theory, which has four supercharges. We denote the parameters of the supersymmetric variation associated with four supercharges as and (. This theory possesses two abelian Rsymmetries: one comes from the original Rsymmetry of the fourdimensional theory which we call and the other comes from a rotational symmetry of the dimensionally reduced space which is called .
We consider the SYM theory on a Riemann surface. It is significant that the topology of the twodimensional curved manifolds without boundary is completely classified by the number of handles (genus) , and once the topology of the Riemann surface is given, we can always choose the metric of the Riemann surface in a specific coordinate patch as conformally flat:
(2.4) 
We thus denote the Riemann surface as and use the conformal metric (2.4) in the following. All the objects of the differential geometry like covariant derivative, spin connection, curvature, and so on are derived from the function in this coordinate patch.
Here we gauge the global symmetry. Recalling that the supersymmetry parameters and have charges and under , respectively, the modified Killing equations (2.3) for the components of become
(2.5) 
in the complex coordinates and . These equations do not have a nonvanishing solution for general , but if we choose
(2.6) 
by using the conformal factor , the spinor
(2.7) 
with a constant Grassmann value becomes a solution of (2.5) since the spin connection and the gauge potential are canceled with each other for . Similarly,
(2.8) 
with a constant Grassmann parameter is a solution of the modified Killing equation. Thus, associated with the Killing spinors (2.7) and (2.8), we have two conserved supercharges on the curved Riemann surface by choosing the background gauge field as (2.6).
For the conserved Killing spinors (2.7) and (2.8), the supersymmetric transformation of the vector multiplets in twodimensions is written as
(2.9) 
In particular, we define a specific linear combination of the supercharges,
(2.10) 
Defining new fermionic fields by
(2.11) 
we see that the supersymmetric transformations under can be simply written as
(2.12) 
It is important that, from the viewpoint of the modified Killing equation, the fermions and can be regarded as 1form and 0forms on , respectively. Using the supercharge , the action can be written in the exact form:
(2.13) 
with . It is easy to see satisfies where is the infinitesimal gauge transformation by the parameter . Therefore the invariance of the action under the transformation is manifest in the expression (2.13).
Note that this choice of the supercharge and the field components is also known as the topological Atwist of theory on . Thus the supersymmetric theory on the curved space with the background Rgauge field can be naturally identified to the topologically twisted theory. However, we emphasize that it does not mean that we must restrict the observables only to the cohomology. We can also regard the theory as a “physical” supersymmetric gauge theory.
Since the theory preserves both of the supercharges and , the action (2.13) is invariant under not only but also another linear combination . The supersymmetry algebra of is obtained roughly by exchanging the role of and from that of reflecting the symmetry, which rotates and , of the theory.
On the other hand, the symmetry acts on the fields as
(2.14) 
and the action (2.13) is manifestly invariant under this symmetry. However symmetry is broken quantummechanically by the anomaly. In fact there is a mismatch in the number of the zero modes of the fermions on unless the Euler characteristic is not equal to zero. This anomaly plays an important role even in the discretized theory discussed in the next subsection.
2.3 Generalized Sugino model as a discretization of the rigid SUSY theory
Let us next consider a discretization of the continuum theory (2.13). Actually, it is already achieved in [47] where the authors gave a theory on an arbitrary discretized Riemann surface (the generalized Sugino model) whose continuum limit is nothing but the theory given by (2.13). Although the original derivation in [47] is based on the topological twisting, we can obtain the same formulation by taking into account the background vector field as discussed in the previous subsections. In fact, we can derive the same formulation by simply replacing all derivatives in the supersymmetric algebra (2.12) and the action (2.13) with difference operators since they do not contain the spin connection in the spin representation as the consequence of the background vector field. Therefore, as we have emphasized, we can regard the discretized Sugino model as a discretized model of the “physical” theory, that is, we do not need to restrict the observables to the cohomology.
Let us briefly review the generalized Sugino model to fix the notation.
This model is constructed on a discretized Riemann surface
which consists of a set of sites, links with directions and faces.
We denote
the number of the sites, the links and the faces as
, and , respectively.
The bosonic variables are written as
and the fermionic variables are written as
where the indices , and
are the labels of the sites, links and faces,
respectively, which stand for constituents of the polygon with which
the variable is associated.
We assume that there is an gauge symmetry
and all the fields take the form of matrix.
In particular, are unitary matrices, are hermitian matrices
and the others are general complex matrices.
In this paper, we regard and as the complex
conjugate with each other
(2.15) 
where is a unitary matrix on the site , and stand for the site labels of the origin and the tip of the link , respectively, and the index in stands for a representative site of the face . We sometimes use the same character to express the representative site of the face . The way of the gauge transformation (2.15) shows that only lives on the link and the other variables live on the sites. We should note that this simply means that the link and face variables except for live on the representative sites of the links and faces.
The supersymmetric transformations of the variables are
(2.16) 
which are the discrete versions of (2.12). Note that satisfies in a parallel way to the continuum transformation. The action of the generalized Sugino model is written in the exact form,
(2.17) 
where
(2.18) 
where is set as
(2.19) 
with
(2.20) 
in order to single out only the physical vacuum of the gauge fields [35]. Here is the face variable corresponding to the face defined as
(2.21) 
where is the number of the links which surround the face , are the labels of those links which construct the face in this order, and () express the directions of the links in the face, which is determined recursively as follows (see Fig. 1): We set to be the representative site of the face . If then and , otherwise and . For , if then and otherwise and .
It is easy to see that the (treelevel) continuum limit of the discretized action (2.17) becomes (2.13) by choosing the parameter , , and appropriately. We can also see that any relevant operator which breaks the symmetry does not appear quantum mechanically by the power counting argument. Therefore we do not need any finetuning to take the continuum limit. For more detail, see [47, 48].
In the actual Monte Carlo simulation, we also have to control the flat directions of the scalar fields and . For this purpose, we add a mass term,
(2.22) 
to the discretized action (2.17). This term explicitly breaks the symmetry, but it is a softbreaking term. We can explicitly evaluate the WT identity associated with the symmetry with dependence as shown in Section 3.2. The explicit form of the action is summarized in the appendix A.
We note that, although the transformation (2.14) is still a symmetry of the discretized action (2.17), the symmetry and the symmetries are broken by the discretization. This is because the discretized version of in (2.19) does not admit the exchange of the role of and . The symmetry and the symmetry will be recovered only in the continuum limit.
To sum up, the action of the generalized Sugino model is given by the summation of the exact part given by (2.17) and the mass term given by (2.22). We divide the exact part into the bosonic part and the fermionic part as
(2.23) 
We further separate the terms including the auxiliary field from as
(2.24) 
with
(2.25) 
Therefore the bosonic part of the whole generalized Sugino model is given by
(2.26) 
Since we use the action after integrating out the auxiliary field in the Monte Carlo simulation, it is convenient to define
(2.27) 
On the other hand, the fermionic part of the action can be written in the form of a bilinear of the fermionic variables:
(2.28) 
where is a vector with the size of which consists of all the elements of the fermionic variables and is an antisymmetric matrix (Dirac operator) with the same size of . By integrating out the fermionic variables , there appears the Pfaffian of ;
(2.29) 
which in general takes a complex value since we are considering a spacetime with the Euclidean signature. We thus write the phase of the Pfaffian as ;
(2.30) 
Since the strategy of the Monte Carlo simulation is to regard the Boltzmann factor of the theory as a probability density, it must be real positive. We thus often approximate the Pfaffian to its absolute value in producing a probability density in the Monte Carlo simulation. Therefore the corresponding partition function in the Monte Carlo simulation of the generalized Sugino model with the standard phasequenched approximation is given by
(2.31) 
where is the measure of the bosonic variables except for the auxiliary field :
(2.32) 
3 anomaly and a novel phasequenched method
In the Monte Carlo computation for the present model, that is, the generalized Sugino model on a discretized Riemann surface, we need to take care of the two deeplyrelated properties: Pfaffian phase and anomaly. In this section, we examine the anomaly in detail and introduce a novel phasequenched method with respect to the anomaly.
3.1 Anomalyphasequench method
Let us consider the measure of the generalized Sugino model
(3.1) 
where
(3.2)  
(3.3) 
Recalling the transformation of the variables (2.14), we see that, even if the action is invariant under the rotation, the measure of the functional integral can have a net charge:
(3.4) 
where denotes the charge of an operator and is the Euler characteristics of the discretized Riemann surface. This nonzero charge of the integral measure corresponds to the lattice counterpart of the anomaly. In the continuum limit, the Euler characteristics is of course intact and the charge of the measure results in the anomaly term in the continuum theory. Since the Pfaffian appears by integrating out the fermionic variables, the Pfaffian has exactly the same charge with that of the integration measure (3.4). Therefore the Pfaffian phase always includes a phase of the origin unless . This is not an artificial phase coming from the Wick rotation but from the topology of the background. This motivates us to decompose the Pfaffian phase into two parts; the anomaly phase and the residual phase,
(3.5) 
where stands for the phase originated in anomaly and is the residual phase. Note that there is always an ambiguity in defining , which we will fix later in (3.11).
An immediate consequence of this observation is that, when the Euler characteristics is not equal to zero, the partition function of the generalized Sugino model without any quenched approximation trivially vanishes:
(3.6) 
because of the unbalance of the number of fermions between the measure and the Boltzmann factor or the nonvanishing background charge. This means that the ordinary definition of the expectation value of an operator ,
(3.7) 
is illdefined.
Although the denominator of the expectation value is usually set so that the total probability becomes unity, it is more important that (the absolute value of) the Boltzmann factor is proportional to the probability density. In this sense the denominator is just a normalization. We thus use (2.31) instead of (3.7) as the denominator of the definition of the expectation value of our model in this paper:
(3.8) 
which is suitable to the Monte Carlo simulation
Incidentally, in the Monte Carlo simulation, we often “approximate” the expectation value (3.8) by replacing the Pfaffian of the numerator to its absolute value:
(3.9) 
which we call the naive phasequenched approximation in this paper.
However, this is not an appropriate approximation of the expectation value
(3.8).
In fact,
from the same reason why the partition function vanishes,
the expectation value (3.8) vanishes unless the charge of
the operator is equal to .
In particular,
if we compute the expectation value of a neutral operator,
it must vanish unless the Euler characteristics of the background is zero ().
However the expectation value of such a neutral operator
in the naive phasequenched approximation (3.9) is apparently
nonzero even if .
Therefore we cannot use (3.9) as
an approximation of the expectation value (3.8)
of the generalized Sugino model
In order to overcome this problem, we introduce a gauge invariant operator with a specific charge . We assume that the operator is invariant not only under any (bosonic) global symmetry transformation of the theory other than but also under the transformation:
(3.10) 
for the later purpose. A notable fact is that if we insert the operator into the path integral, it exactly cancels the charge of the integration measure or the Pfaffian. In this sense, we call the operator the “compensator”. We can define the part of the Pfaffian phase in (3.5) through the compensator as
(3.11) 
Once we define the anomaly phase , we can introduce a phasequench method of a novel type by ignoring only the residual phase of the Pfaffian phase in (3.5):
(3.12) 
which we call the anomalyphasequench method. We see that vanishes unless the operator has the charge , which is exactly the property of the expectation value of the generalized Sugino model.
Furthermore, we expect that there is no sign problem in the present
discretized theory as long as the anomalyinduced phase is cancelled by the compensator.
In other words, the residual phase of the Pfaffian in (3.5)
will not contribute to the results of the expectation values.
As shown in [33],
this is actually the case in the torus background
where it was shown that the Pfaffian becomes real positive
in the continuum limit.
The point is that,
for the continuum 2D SYM in the flat background,
the eigenvalues of the Dirac operator always appear
as complex conjugate pairs [33].
We can apply the same logic to the continuum theory on a curved
background discussed in Section 2,
and the product of the nonzero eigenvalues of the Dirac operator
is expected be real positive for the continuum theory in a curved background as well.
Recalling that the anomaly of the continuum theory comes from
the fermionic zeromodes,
this strongly suggests that there is no sign problem in the present theory
as long as we eliminate the anomalyinduced phase
We close this subsection by making three comments. First, let us consider the expectation value
(3.13) 
for an arbitrary operator with zero charge . Since the combination has the charge , it has a nontrivial value in general. We see that the anomalyphasequenched approximation of (3.13) can be expressed in the naivephasequenched approximation by
(3.14) 
which is useful when we evaluate (3.13) in the Monte Carlo simulation.
Second, the compensator is not unique but has many possible choices. We propose the following three different compensators among them; the tracetype compensator,
(3.15) 
the IzyksonZuber(IZ) type compensator,
(3.16) 
and the determinanttype compensator,
(3.17) 
Here we should note that there is an ambiguity in the branch of the exponents when they are fractional numbers. We choose the branch for and in this paper. It is notable that an appropriate type of compensators depends on topology of the background space, which we will discuss in details in the next section.
Third, there is no natural principle to fix the overall sign of the Pfaffian. In fact, we can easily see that the sign of the Pfaffian can flip by changing the order of the fermionic variables in defining the fermionic action. (Note that this is not a sign problem, but just ambiguity of a common sign of the Pfaffian.)
3.2 WardTakahashi identity
As a candidate of the measurements in the Monte Carlo simulation,
we derive a WT identity corresponding to the symmetry
Let us consider the integral,
(3.18) 
where in given in (2.23), is a compensator and is a constant Grassmann number with the charge . Since the combination is neutral under the transformation, the integral does not have the charge. The integral itself is not altered even if we change all the variables in the expression (3.18) to . The expansions of and by are
(3.19) 
while the measure and the compensator are invariant. We then obtain the relation among the expectation values:
(3.20) 
where we have integrated out the auxiliary field and divided the integral by . Since the expectation value of the fermion action is related to the degrees of freedom of the fermions of the system. In fact, we can show in general that
(3.21) 
where is an arbitrary function, is a vector of all the fermion degrees of freedom and is an antisymmetric matrix with the size of the fermion degrees of freedom. Therefore, we can further estimate as
(3.22) 
for the tracetype and determinanttype compensators and
(3.23) 
for the IZtype compensator. We then finally obtain the WT identity,
(3.24) 
for the tracetype and determinanttype compensators and
(3.25) 
for the IZtype compensator. We emphasize that this relation is exact even if it includes the parameter of the SUSY breaking term. In the Monte Carlo simulation, we check if these identities are satisfied in the anomalyphasequenched approximation (3.14).
In the next section, we evaluate the left hand side of the equation (3.24) for () and () and the equation (3.25) for in the anomalyphasequenched approximation. It gives a nontrivial check for the following claims:

The generalized Sugino model on the curved space is a valid extension of the Sugino model on the flat space and the model correctly reproduces the anomaly through a unbalance of the fermion number due to the topology of the background spacetime.

The quench of the phase among the whole Pfaffian phase gives a correct result in the lattice simulation on the generic background space. It means that is negligible and the sign problem due to the phase is absent in the generalized Sugino model on the generic background spacetime.
4 Monte Carlo Simulations
In this section, we show the results of the Monte Carlo simulation of the generalized Sugino model on several polyhedra with and . We first check if the flat directions of the scalar fields are properly controlled by the mass term (2.22). We next show that the anomalyphasequenched approximation properly works by evaluating the lefthand sides of (3.24) and (3.25) in this approximation. We investigate the behavior of the Pfaffian phase in details and show that we have no further sign problem in this model as long as the anomalyinduced phase is cancelled by the compensator. We also check that the origin of the anomaly originates from the light modes (pseudozeromodes) of the Dirac operator in the discretized model.
4.1 Algorithm and data analysis
We have used the rational hybrid Monte Carlo algorithm [52], where we ignore the phase of the Pfaffian and estimate the absolute value of the Pfaffian using the pseudofermion method by approximating the matrix by rational functions,
(4.1) 
We have set and the coefficients and are determined by the Remez algorithm. For each set of parameters, we have generated 15000 – 20000 configurations in general, while we have generated 80000 – 100000 configurations for some parameters to obtain a sufficient statistics. The simulation code is written in FORTRAN 90 and is not yet parallelized, which had run on personal computers with Intel Core i7. After measuring the autocorrelations, we have estimated the standard error by using the Jack Knife method.
4.2 Simulation parameters
In this paper, we consider only polyhedra with the same shape of the face as discretization of Riemann surfaces. For , we consider the five regular polyhedra, namely, the regular tetrahedron, the regular octahedron, the regular cube, the regular icosahedron and the regular dodecahedron. For , we consider the regular square lattice with the size of , and . For , we consider such a double torus which is realized by gluing two regular square lattice with the size by one surface (see Fig.2).
We fix the physical ’tHoot coupling as , thus the dimensionless coupling can be expressed as with the lattice spacing (we note and ). We identify the lattice spacing with length of a link and fix the surface area of background as . Hence, the lattice spacing is given by
(4.2) 
where is area of a face with unit sides,
(4.3) 
The continuum limit can be realized by . We summarize geometries used for our simulations in Table 1.
genus  Euler ch.  geometry  shape of face  lattice spacing  

2  tetra  4  6  4  T  0.7598  
octa  6  12  8  T  0.5373  
cube  8  12  6  S  0.4082  
icosa  12  30  20  T  0.3398  
dodeca  20  30  12  P  0.2201  
0  reg.lat.  9  18  9  S  0.3333  
reg.lat.  16  32  16  S  0.2500  
reg.lat.  25  50  25  S  0.2000  
2  Fig.2  14  32  16  S  0.2500 
The gauge group is restricted to and thus the parameter appeared in (2.19) to single out the physical vacuum was set to . For the scalar mass to stabilize the flat directions of , we have chosen the dimensionless mass parameter as appropriate values for each topology.
Recall that the role of the compensator is to cancel the nonvanishing charge of the pathintegral measure. Since the unbalance of the charge is caused by the fermions, we can effectively regard the insertion of the compensator as the insertion of the fermions. On the background, the number of the fermion modes of and in the measure exceeds that of . Since the phase of and is related with that of the bosonic field by the symmetry, it is a natural and good choice to consider the compensator composed only of in terms of numerical efficiency. We then adopt for (sphere). On the background, the situation is opposite. In this case, we need to choose a compensator containing itself, so we adopt in order to cancel the phase efficiently for (double torus). Since and are identical up to sign for the gauge group , we do not use in this paper.
<