Matrix model and Yukawa couplings on the noncommutative torus

The IKKT model is proposed as a non-perturbative formulation of superstring theory. We propose a Dirac operator on the noncommutative torus, which is consistent with the IKKT model, based on noncommutative geometry. Next, we consider zero-mode equations of the Dirac operator with magnetic fluxes. We find that zero-mode solutions have the chirality and the generation structures similar to the commutative case. Moreover, we compute Yukawa couplings of chiral matter fields.


Introduction
Superstring theory is a promising candidate for a unified theory of all forces in nature. Superstring theory is only defined perturbatively and has infinite degenerate vacua. Therefore, it is said that superstring theory has no predictions for our world, and we need a more fundamental theory.
Matrix models are proposed as a non-perturbative formulation of superstring theory. In this paper, we focus on Ishibashi-Kawai-Kitazawa-Tsuchiya (IKKT) model [1]. This model is derived from the matrix regularization of the Green-Schwarz action in the Schild gauge or large-N reduced model of ten-dimensional (10D) super Yang-Mills (SYM) theory.
In the IKKT model, matter fields and degree of freedom of spacetime are considered to be embedded in matrices. Several attempts have been made to show it. In ref. [2], the authors considered an intersecting flat D-branes system based on the IKKT model. To analyze the chirality and the generation structures of the system, they used the analogy of the harmonic oscillator in quantum mechanics. They showed the existence of a chiral JHEP04(2019)079 zero-mode which are coming from a string connecting two different D-branes. However, they considered a simple configuration, e.g., two D-branes are orthogonal. Therefore, they could not realize multiple chiral zero-modes. In addition, in more general configurations, we cannot easily find the chiral zero-modes by using the same method because the complicated contribution comes from the mixing term between two D-branes. In ref. [3], the authors considered intersecting fuzzy spheres as compact D-branes and realized two chiral zeromodes. To admit a fuzzy sphere as a classical solution, we must introduce a new term. However, the origin of such a term is not clear. In either case, it is difficult to obtain the number of chiral zero-modes that we expect. In ref. [4], the authors challenged the realization of three generations by using numerical analysis. By considering the squashed fuzzy sphere in addition to the fuzzy spheres, they succeeded in realizing three generations numerically. On the other hand, the authors of ref. [5], considered a fuzzy torus with a magnetic flux based on the finite-matrix formulation of gauge theories. 1 They computed the index of the overlap Dirac operator on the fuzzy torus by the Monte Carlo simulations and showed the numerical results that are consistent with the index theorem. In this formulation, the exact relationship with the IKKT model is not clear because we must introduce the special type of the Dirac operator to obtain the non-trivial index of the Dirac operator. Moreover, it is still difficult to compute physically important quantities like Yukawa couplings.
The purposes of this paper are to define a Dirac operator that is consistent with the IKKT model, to analyze the chirality and the generation structures and to compute Yukawa couplings of chiral matter fields.
To analyze concretely, we consider the noncommutative (NC) torus as a classical solution of the IKKT model. It is known that the NC torus as an irrational rotation ring can be realized in the IKKT model [6]. In addition, we consider the analogy of the toroidal compactifications of 10D SYM theory. In ref. [7], the authors considered the toroidal compactifications of 10D SYM theory with magnetic fluxes. A key concept is the twisted bundle. The twisted bundle can be interpreted as a compatibility condition between the periodic boundary conditions on the torus and the gauge transformations. This compatibility condition implies that the magnetic fluxes are quantized, then the zero-mode solution of the Dirac equation can have the chirality and the degeneracies. The authors of ref. [7] identified these results as the chirality and the generation structures in the four-dimensional effective theory and computed the Yukawa couplings by computing the overlap integral over the torus. From the viewpoint of phenomenology in superstring theory, this result is very important. Therefore, we consider the twisted bundle on the NC torus.
In the following, we propose a Dirac operator on the NC torus that is consistent with the IKKT model, and we solve the zero-mode equation of this Dirac operator with magnetic fluxes. Then, we compute the normalization factors of the zero-mode solutions and Yukawa couplings of chiral matter fields.

JHEP04(2019)079
The organization of this paper is as follows. In section 2, we briefly review the IKKT model and the realization of the NC torus in the IKKT model. In section 3, we review the basic results of the toroidal compactifications of 10D SYM theory. In section 4, we introduce differential operators on the NC torus based on noncommutative geometry. In addition, we propose a Dirac operator that is consistent with the IKKT model. In section 5, we solve the zero-mode equation by using the analogy of the Fourier transformation. In addition, we compute the normalization factors of the zero-mode solutions and the Yukawa couplings of chiral matter fields. To compute, we define the trace that is consistent with the gauge transformations and the torus translations. Section 6 contains conclusions and discussion.

IKKT model
The action of the IKKT model is defined as follows [1], where X M (M = 0 ∼ 9) is a 10D vector and N ×N Hermitian matrix, ψ is a 10D Majorana-Weyl spinor whose components are N × N matrices and g is a scale factor. Indices are contracted by the Minkowski metric. On the other hand, we will see that this action admits an infinite dimensional representation like linear operators in the next subsection. The action (2.1) has some symmetries: where i (i = 1, 2) is 10D Majorana-Weyl spinor as a Grassmann odd parameter, c M is a 10D constant vector, 1 N is the N × N identity matrix, and Λ is a N × N Hermitian matrix. Eqs. (2.2) and (2.3) can be identified N = 2 supersymmetry (SUSY) as follows. If we consider a linear combination of δ (1) and δ (2) as where we used eq. (2.5) and the equation of motion of ψ. If we identify X M as a 10D spacetime coordinate, eq. (2.6) is N = 2 on-shell SUSY algebra. In this sense, we often say that the degree of freedom of spacetime is embedded in matrices. From this, there is a possibility of analysis for the dynamics of spacetime, e.g., a mechanism of compactifications.

JHEP04(2019)079
When ψ = 0, the equation of motion of X M is The simplest solution is ∀ [X M , X N ] = 0. However, we can show that attractive forces act between the eigenvalues in the one-loop effective potential around this vacuum. Therefore, they do not spread, and the correspondence with the original theory does not hold. In this case, we need more conditions on gauge groups [8,9]. The second simplest, however, an interesting solution is that some commutators are proportional to the identity matrix, i.e., where θ M N is a real anti-symmetric matrix with respect to the Lorentz indices. We omit the identity matrix. When each matrix has a finite size, eq. (2.8) is a contradiction. We should interpret eq. (2.8) as it is satisfied at large matrix size. This is the correspondence between a function algebra and a matrix algebra in the matrix regularization. On the other hand, the infinite dimensional representation, which we are interested, strictly respects eq. (2.8).
We can also interpret eq. (2.8) as a D-branes (BPS) configuration as follows. We can find that the transformation by δ (1) is proportional to the identity matrix and cancel the transformation by δ (2) . Explicitly, if we set 2 Therefore, the half of N = 2 SUSY is preserved in this background, and the one-loop effective potential exactly vanishes due to SUSY. It is an advantage of this model that a many-body system of D-branes can be realized. When we consider two sets of eq. (2.8), is also a solution of eq. (2.7). In general, eq. (2.9) is not a D-branes configuration. However, some configurations become a D-branes configuration again, e.g., the case of X M is interpreted as two coincident D-branes, and other cases are considered in ref. [1].

NC torus in the IKKT model
In this subsection, we briefly review the NC torus in the IKKT model based on refs. [6,10]. A unitary transformation is defined by where U is a unitary matrix (operator). If the trace in the action (2.1) has the cyclic property, this unitary transformation becomes a symmetry of the IKKT model. The infinitesimal form of eq. (2.10) corresponds to eq. (2.5). The authors of ref. [6] showed that the NC torus can be realized by this unitary transformation.

JHEP04(2019)079
We restrict the action (2.1) to the subspace of (X M , ψ) belonging to the same gauge class before and after translations in the directions X 4 and X 5 . In other word, we consider conditions as follows, where R i (i = 4, 5) is the scalar matrix (operator) of a real coefficient corresponding to the periods of the torus. In the following, let us define M = 0 ∼ 9, µ = 0 ∼ 3 and i = 4 ∼ 9.
We can see that any finite dimensional representations do not satisfy the conditions (2.11). However, we can find solutions if we interpret X M and ψ as operators on an infinite dimensional Hilbert space. 2 From eq. (2.11), we can confirm that U 4 U 5 U −1 4 U −1 5 commutes with X M and ψ. Therefore, we assume (or we may be able to apply Schur's lemma for countable dimensional cases [11,12] because we can also realize U 4 and U 4 as operators on l 2 (R)) that is a scalar operator, i.e., where θ is a real parameter. Eq. (2.12) is the algebra of the NC torus in mathematics. Precisely, we can restrict θ ∈ [0, 1/2] because of some isomorphisms. Depending on whether the parameter θ is a rational number or an irrational number, the mathematical structure differs. In fact, the algebra of the NC torus admits a finite dimensional representation if θ ∈ Q, and we cannot realize such algebra by the finite dimensional representation if θ ∈ R\Q. In this paper, we assume θ ∈ R\Q because we are interested in the infinite dimensional representation of the NC torus. In this sense, X M , ψ, U 4 and U 5 are no longer finite dimensional matrices but operators. We use "hat" to indicate that it is an operator (we omit the identity operator).
In the following, we show a concrete representation ofX M ,Û 4 andÛ 5 . Let us start from a Hilbert space H = L 2 (R) ⊗ C m , m ∈ N (m corresponds to Z m in the usual construction [6,10,13]). We can realizeÛ 4 andÛ 5 as operators on H, i.e., where n is some integer. For simplicity, we assume that m and n are positive and co-prime each other (or m = 1 and n = 0). We defineX 4 andX 5 aŝ where 1 m is the m×m identity matrix,x andp are the position operator and the momentum operator on one-dimensional quantum mechanics (where =1), respectively. We can see

JHEP04(2019)079
that θ 45 = 2πm n−mθ R 4 R 5 and easily confirm that eqs. (2.11) and (2.12) are satisfied. In the following, since the part acting on C m is the identity matrix, we use the Hilbert space From an analogy of quantum mechanics, since the periodicity is realized by the unitary operatorsÛ 4 andÛ 5 , and we restrict the fundamental (or physical) region to 0 ≤ X 4 ≤ 2πR 4 and 0 ≤ X 5 ≤ 2πR 5 . 3

Toroidal compactifications of 10D SYM theory
In this section, we review the basic results of toroidal compactifications of 10D SYM theory based on refs. [7,14].

Twisted bundle on the torus
Let us start from U(1) gauge theory on the torus T 2 . In this paper, we take a lattice Λ Z 2 that is generated by v 1 = (2πR 4 , 0) and v 2 = (0, 2πR 5 ), and T 2 is defined by R 2 /Λ.

A constant magnetic flux is introduced by
then the magnetic flux F 45 = F . We use the same gauge choice (axial gauge in ref. [14]) in section 5. We must confirm that the gauge theory is well-defined on T 2 . Obviously, we should confirm only about x 4 -direction. The covariant derivatives transform like If this translation can be absorbed as the gauge symmetry, U(1) gauge theory on T 2 is well-defined. Actually, we can realize eq. (3.2) by In general, we need to consider another gauge transformation Ω 5 (x 4 , x 5 ) associated with the translation along the x 5 -direction. In the above case, Ω 5 (x 4 , x 5 ) = 1. The above concept is called the twisted bundle.
In addition, we must confirm a consistency condition 3 Let Hm (m ∈ Z) be a Hilbert space such that the spectrum of the Hermitian operatorX 4 lies in the interval 2πmR4 ≤ X 4 ≤ 2π(m + 1)R4. In this case, Hm is unitary equivalent to ∀ H n∈Z Therefore, we interpret the fundamental (or physical) region with respect to the X 4 -direction is 0 ≤ X 4 ≤ 2πR4. We interpret the X 5 -direction as well.

JHEP04(2019)079
where A = (2π) 2 R 4 R 5 is the area of the torus. Namely, the magnetic flux on the torus is quantized.

Zero-modes of Dirac operator on T 2
We can analytically solve the zero-mode equations of the Dirac operator with the background gauge field (3.1). The Dirac operator is defined by where σ 4 and σ 5 are the Pauli matrices σ 1 and σ 2 , respectively. For convenience, we rewrite the constant magnetic flux as F = 2π A N ν, where N is a positive integer and ν = sign(F ). The zero-mode equation of the Dirac operator is If ψ is labeled by the eigenvalue of the chirality matrix σ 3 , i.e., the zero-mode equations of each component can be written as We can easily factorize ψ s satisfying eq. (3.6) as is periodic with respect to the x 5 -direction. Therefore, we can use the Fourier expansion as From the first condition of eq. (3.7), the Fourier coefficients satisfy

JHEP04(2019)079
If we decompose n into n = N p + q (p ∈ Z, q = 0, . . . , N−1), thus we can rewrite eq. (3.8) as From the above, f s (x 4 + isx 5 ) can be written as If νs = +1, f s (x 4 +isx 5 ) converges. In other words, the sign of the magnetic flux determines the chirality of the zero-mode solution. In addition, the summation with respect to q means that there are N -independent zero-mode solutions. In ref. [7], these results are interpreted as the chirality and the generation structures in the four-dimensional effective theory. Therefore, each degenerated solution of eq. (3.6) is Each C s q is determined by the normalization condition, i.e., From eq. (3.9), Eventually, where the function ϑ is the Jacobi theta-function where Imτ > 0. We can easily generalize the above discussion to non-Abelian, bi-fundamental or higher dimensional torus.

JHEP04(2019)079 4 Dirac operator on the NC torus and the IKKT model
In this section, we introduce differential operators on the NC torus and construct them from the IKKT model.

Differential operators on the NC torus
Let A be a C * -algebra over C and d : A → A be a linear map. d is called a derivation on A if d satisfies the Leibniz rule, i.e., In mathematics, the NC torus is the C * -algebra generated byÛ 4 andÛ 2 satisfying eq. (2.12). Therefore, derivations are completely defined by how they act on the two generators. We define and satisfy the linearity and the Leibniz rule. δ i (i = 4, 5) is called the basic derivation. The factors 1 R 4 and 1 R 5 come from the definition of the torus, cf., subsection 3.1. We can confirm that the basic derivations commute each other.
The Dirac operator on the NC torus is defined by the basic derivations where σ 4 and σ 5 are the Pauli matrices σ 1 and σ 2 , respectively.

Dirac operator on the NC torus based on the IKKT model
In this subsection, we construct the differential operators, which are introduced in the previous subsection, based on the IKKT model. First, we focus on noncommutative SYM theory based on the IKKT model [15]. We consider the expansion of the action (2.1) around the specific background: space-timeX M bg satisfying eq. (2.8), gauge fieldÂ bg,M will be defined in section 5. Namely, We refer to appendix A for the details. A point is that the partial derivatives are defined as where (θ −1 ) M N is the inverse matrix of θ M N . We can confirm that the partial derivatives satisfy ∂ MX N bg = δ N M . Here, we assume that θ M N is non-degenerate. In the following, we identifyX M bg with eq. (2.14) (i.e., we mainly focus on the NC torus (M = 0 ∼ 9 → i = 4, 5)), and we omit the subscription "bg".

JHEP04(2019)079
The action of ∂ i on the algebra of the NC torus is Therefore, the basic derivations on the NC torus based on the IKKT model are defined as (4.4) We can confirm that the basic derivations (4.4) satisfy eq. (4.1). In addition, eq. (2.8) and the Jacobi identity assures the commutativity of the basic derivations (4.4).
From the above, we propose a Dirac operator on the NC torus that is consistent with the IKKT model as follows, 4 In the fermionic part of the action (2.1), is identified as a Dirac operator on the extra dimensions, e.g., [17]. As mentioned above, we focus on the first two extra dimensions that are the NC torus (i.e., 9 i=4 → i=4,5 ). We can verify that eqs. (4.5) and (4.6) are equivalent. Namely, Therefore, the operator (4.5) is a suitable Dirac operator on the NC torus based on the IKKT model. The partial derivatives are defined by eq. (4.3), and the gauge field is introduced by eq. (4.2). The covariant derivatives are naturally defined by However, we cannot introduce any fermionsψ in a fundamental representation, i.e.,ψ → Ωψ, whereΩ is a unitary operator as a gauge transformation. This comes from eq. (2.10).

JHEP04(2019)079
On the other hand, we can introduce fermions in a fundamental representation if we define the covariant derivatives as However, the original action (2.1) does not have the field strength which is defined by the covariant derivatives (4.8). Therefore, if we want to consider fermions in a fundamental representation, we should realize as a part of the action (2.1). Let us consider, for example, U(2) gauge theory with an adjoint matter based on the action (2.1), i.e., We assume that the gauge group breaks such that U(2) → U(1) × U(1) → U(1): top left component ofÂ i . If we focus on onlyψ 12 and ignore other components ofÂ i , then we can derive the action with the fermion in U(1) fundamental representation.
5 Zero-mode analysis on the magnetized NC torus

Twisted bundle on the NC torus
To consider the twisted bundle on the NC torus, we introduce a background gauge field A 4 (X 4 ,X 5 ) = 0,Â 5 (X 4 ,X 5 ) = FX 4 , (5.1) and the field strength isF 45 = F . The gauge transformation is obtained bŷ whereΩ is a unitary operator. The background gauge field is varied by the torus translations (2.11), i.e., A 4 (X 4 + 2πR 4 ,X 5 ) =Â 4 (X 4 ,X 5 + 2πR 5 ) = 0, We can realize eq. (5.2) as the gauge transformations, i.e., where ∝ represents an action on C m part. This part does not affect eq. (5.2) since the C m part ofX 4 andX 5 is the identity matrix. Therefore, in the following, we assume this part is the identity matrix and omit.

Zero-modes of the Dirac operator / D phys
In this subsection, we consider the zero-mode equations of the Dirac operator / D phys with the background gauge field (5.1). In the following, we show the zero-mode solutions in (i) fundamental representation (ii) bifundamental representation.
Case (i). In this case, the zero-mode equation is whereψ s (s = ±1) is in the fundamental representation. Simply, with the periodic boundary conditionŝ ψ s (X 4 + 2πR 4 ,X 5 ) = exp 2πiR 4 1 + θ 45 F · FX 5 ψ s (X 4 ,X 5 ) ψ s (X 4 ,X 5 + 2πR 5 ) =ψ s (X 4 ,X 5 ) Although eq. (5.6) is written by the operators, this equation has the same form with eq. (3.6). Therefore, we expect that the zero-mode equation (5.6) can be constructed from that of eq. (3.6). Here, we consider the Fourier transformations. The Fourier transformation of the whole zero-mode solution of eq. (3.6) is obtained by where C is the normalization constant and the magnetic flux F satisfies the quantization condition (3.4). We expect that the whole solution of eq. (5.6) is obtained bŷ However, eq. (5.8) does not satisfy eq. (5.6) because the quantization conditions for the magnetic flux are different.
Here, we focus on that p and q do not appear alone in eq. (5.7). This fact allows us to vary the decomposition of the label from n = |N |p + q to n = |N |p + q. In addition,

JHEP04(2019)079
we should replace the magnetic flux F with F satisfying (5.5). Then, we can obtain the zero-mode solutions, 5 i.e., where q = 0 ∼ |N | − 1. We can confirm that the above expansions satisfy eq. (5.6) if F s > 0. The normalization constant C will be computed in subsection 5.4. However, this form is difficult to use when we compute the normalization constant and Yukawa couplings. Fortunately, we can rewrite by using the Jacobi theta-function If we require sN > 0, eq. (5.10) can be interpreted as the operator form of eq. (3.11) up to the normalization constant. Therefore, in the following, we restrict ourselves to the case sN > 0.
In this case, we consider a fermion in U(2) adjoint representation, The zero-mode equation for each chirality is written as First, we consider the top right componentψ s 12 . Precisely, the zero-mode equation of ψ s 12 is

Eigenvalues of the Laplacian
In the IKKT model, ) is identified as the Laplacian (e.g., [2,3] or eq. (A.6)). The eigenvalue problem of the Laplacian relates to that of the square of the Dirac operator. Let us rewrite the Dirac operator as . (5.14) The action of (5.14) on (1, 2) component of the fermion is written as From the above, D (D) has the zero-mode solutions if F 12 > 0 (F 12 < 0). Therefore, we can see that eq. (5.12) is not only the zero-modes of the Dirac operator but also the lightest mode of the Laplacian. In addition, we can construct the eigenmodes of the Laplacian corresponding (1, 2) component. We focus on the commutation relation of D and D , If we select F 12 > 0, D and −D can play roles of the creation operator and the annihilation operator, respectively. By considering the analogy of the harmonic oscillator, i.e., then the eigenmodes of the Laplacian are obtained by ψ +,n 12 = λ nψ +,n 12 , ψ +,n 12 := D nψ+ 12 , λ n = (θ 45 ) 2 F 12 (2n + 1).
If we select F 12 < 0, it is sufficient to reverse the roles of D and D .
In the commutative case, 6 the spectrum of the Laplacian is obtained by λ n = F 12 (2n + 1), where F i (i = 1, 2) is the magnetic flux satisfying (3.4) (cf., [7]) and F 12 = F 1 − F 2 . It seems that the spectrum vanishes in the limit θ 45 → 0. However, the bosonic part of the effective action is the fourth order for the NC parameter if we ignore the order included in the definition of the partial derivatives. Therefore, (θ 45 ) 2 of the spectrum is necessary, but the limit θ 45 → 0 is non-trivial.

Normalizations and Yukawa couplings
Normalizations. Let us start from the normalization constant of eq. (5.12) since F 2 = 0 corresponds to the case (i).

JHEP04(2019)079
where is a function satisfies I 12 (0) = π/|F 12 | (the Gaussian integral). Therefore, the normalization constant C 12 is We should note the definition of the trace on the infinite dimensional space. For example, In addition, we should confirm that the trace is well-defined with respect to the compactification conditions (2.11). In this paper, target operatorsf are periodic with respect to the X 5 -direction. This means On the other hand, the target operatorsf are quasiperiodic with respect to the X 4direction. Since the unitary operator for the quasiperiodicityΩ 4 is written byX 5 (Ω 4

JHEP04(2019)079
is also periodic with respect to the X 5 -direction), then where we used the equivalence (5.15) in the fourth line. For general gauge transformationsÛ , we can show the equivalence between before and after the gauge transformations if we assume the existence of the completeness relation ofÛ . Therefore, the cyclic property of the trace is held, at least, for the gauge transformations. Therefore, the gauge symmetry of the action (2.1) is still held. Similarly, we can verify that all traces defined by the appropriate completeness relation are equivalent.
Yukawa couplings. In the following, for simplicity, we consider the magnetic fluxes break the gauge group U(N ) → 3 a=1 U(N a ), where 3 a=1 N a = N to compute the Yukawa couplings. The background gauge field is obtained bŷ where F i (i = 1, 2, 3) satisfies (5.5) and 1 N i is the N i × N i identity matrix. We should consider a sign assignment of the magnetic fluxes. In the following, we select F 23 , F 21 , F 13 > 0 (This implies s 23 , s 21 , s 13 = +1 and N 23 , N 21 , N 13 > 0). This sign assignment is justified by the relation F 12 + F 23 + F 31 = 0.
The component of the zero-mode fermions corresponding the above sign assignment is obtained byψ where i, j and k denote the degeneracies and const. =1 (corresponds to the wavefunctions of the gauginos in the commutative case). The lightest mode bosons have the same matrix structure with the (anti-) chiral zero-mode fermions. Let us denote by Φ ab,I (a, b = 1, 2, 3) the (a, b) block component of the zero-mode fermions or the lightest mode bosons since their function form is the same. In the action (2.1), the Yukawa couplings are described as the product of three matrices, i.e., where we used s 21 , s 13 = +1 and the product formula of the Jacobi theta-function [18] ϑ On the other hand, the orthogonality ofΦ ab,I is assured in the above result. Therefore, the Yukawa couplings (5.16) is obtained by where we assume ∃ m ∈ Z |N 21 |+|N 31 | such that K = I + J + |N 21 |m.
The Yukawa couplings (5.17) differ from the commutative case (cf., [7]) by the overall factor only if we fix the generation numbers N 23 , N 21 and N 13 . On the other hand, we can confirm that the Yukawa couplings (5.17) go back to the commutative case in the limit θ 45 → 0 since the normalization constant goes back to the commutative case. 7 However, the limit θ 45 → 0 is non-trivial since the NC parameter θ 45 remains in the Yukawa couplings, cf.,Γ i := θ ij Γ j in eq. (A.9). 7 The normalization constants are slightly different between the commutative case and our case. This is because the integral of 1 is normalized to obtain the area of the torus in the commutative case. In our case, the integral of1 is normalized by considering the analogy of quantum mechanics (2.15). This normalization is natural from the viewpoint of the NC torus without magnetic fluxes. Therefore, if we consider the integral of 1 is normalized to obtain 1 in the commutative case, the normalization constants are the same between the commutative case and our case.

JHEP04(2019)079 6 Conclusions and discussions
In this paper, we performed the analysis of the chirality and the generation structures on the magnetized NC torus based on the IKKT model. In subsection 4.2, we proposed the suitable Dirac operator on the NC torus by considering noncommutative geometry. In section 5, we analyzed the zero-mode solutions of the Dirac operator we proposed. We showed that zero-mode solutions have the chirality and the generations structures. In addition, we computed the Yukawa couplings of chiral matter fields. Compared with the commutative case, the difference of the Yukawa couplings is the overall factor only. Advantages of our method are (i) we can consider geometric conditions such as periodic boundary conditions (ii) we can write down the analytic form of zero-mode solutions which can easily be compared with the commutative case. This is important to observe NC effects from the IKKT model.
When we consider the microscopic world, we compute the physical quantities through functions on spacetime like wavefunctions. On the other hand, in noncommutative geometry, we consider a function algebra on a certain space which has a NC product. For example, the star-product is a NC product in the context of deformation quantization. In the relationship between analytical mechanics and quantum mechanics, this corresponds to the replacement of the coordinate of the phase space by operators. Therefore, in the sense of noncommutative geometry, we can admit that the chirality and the generation structures of our zero-mode solutions have the physical meanings even though zero-mode solutions are written by the operators.
From a Phenomenological point of view, our results may not be new because the difference of the Yukawa couplings is the overall factor only. However, the important point of this paper is that the IKKT model can describe the string-motivated model including the NC effect. The IKKT model is considered as a non-perturbative formulation of superstring theory. Therefore, we expect that the IKKT model should describe the results of stringmotivated models known so far. From this point of view, our results are important.
We are interested in generalizations: orbifolds, complex structure + Wilson line and the origin of the magnetic flux.
The toroidal orbifolds are typical models in string phenomenology. In magnetized toroidal orbifolds, the generation structure differs from the toroidal compactifications without orbifold projections. In refs. [19,20], the NC toroidal orbifolds are considered.
On the other hand, in general, the Yukawa couplings in the toroidal compactifications are functions of the complex structure moduli and the Wilson lines. Therefore, the values of the complex structure moduli and the Wilson lines are important to compare with the observed values (cf., [21]). Our results correspond to the case whose complex structure τ is τ = iR 4 /R 5 . We can consider the complexification when we introduce the basic derivations. We expect that the Yukawa couplings on the magnetized NC torus are obtained by the general complex structure τ instead of τ = iR 4 /R 5 and the overall factor including τ .
In our results, the magnetic flux played an important role. However, we introduced the magnetic flux by hand. We expect that the magnetic flux is also generated from the dynamics of the IKKT model. Recently, in refs. [22,23], the authors showed that the JHEP04(2019)079 magnetic flux may come from the tachyon condensations induced from the dynamics of D-branes and non-BPS D-branes. We expect that our results can be described by the full dynamics of the IKKT model.
Gauge selection. In this paper, we selected the axial gauge (5.1). Our method, especially the ansatz of the zero-mode solutions which are similar to the Fourier transformation, depends on the gauge selection. Therefore, we should confirm the gauge invariance of our results.
In the commutative case, the background gauge field with the fixed magnetic flux F is obtained by where t ∈ [0, 1]. We can realize the gauge transformation On the other hand, the background gauge field on the NC torus with the fixed magnetic flux F is obtained byÂ where t ∈ [0, 1]. We expect that the gauge transformation from ∀ t 1 to ∀ t 2 by the unitary operator, at least,Û = exp iα(X 4X 5 +X 5X 4 ) , where α ∈ R. However, we can find a condition such that If we fix the starting point t = t 1 , then t 2 must be t 1 or 1 − t 1 since eq. (6.1) is a quadratic equation with respect to the t 2 . We need to confirm whether the gauge transformations from ∀ t 1 to ∀ t 2 exist. This situation is the same with a gauge theory on a NC space with the star-product formulation. This is an open question. If we cannot find, different gauge backgrounds may correspond to physically different theories.
The limit: θ 45 → 0. The zero-mode solutions (5.10) and (5.12) become those of the commutative case if we consider the limit θ 45 → 0 and assume the operators (X 4 ,X 5 ) correspond to the coordinate on the torus (x 4 , x 5 ). In addition, the Yukawa couplings (5.17) are the same with the commutative case in this limit. However, from the discussion in ref. [15] (or appendix A), the NC parameter remains in the effective action, e.g., the effective metric. Therefore, in terms of the effective action, the limit θ 45 → 0 is non-trivial. This is also an open question for us.

A Effective action of the IKKT model
We derive the effective action of the IKKT model by considering eq. (4.2). We refer to [7] for the basic techniques. In this appendix, we consider the whole spacetime again.
The action of the IKKT model is We introduced the gauge field as a fluctuation (4.2). By substituting,  Similarly, we can obtain the fermionic part of the effective action where the indices are contracted by the Minkowski metric andΓ µ := θ µν Γ ν , and S f,other contains irrelevant terms for our main discussions.
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.