arxiv: v2 [hep-th] 2 May 2017

Transcription

1 Prepared for submission to JHEP he Flavoured BFSS Model at High emperature arxiv: v2 [hep-th] 2 May 217 Yuhma Asano, a Veselin G. Filev, a Samuel Kováčik, a,b Denjoe O Connor, a a School of heoretical Physics, Dublin Institute for Advanced Studies, 1 Burlington Road, Dublin 4, Ireland. b Faculty of Mathematics, Physics and Informatics, Comenius University Bratislava, Mlynská dolina, Bratislava, , Slovakia. yuhma@stp.dias.ie, vfilev@stp.dias.ie, samuel.kovacik@fmph.uniba.sk, denjoe@stp.dias.ie Abstract: We study the high-temperature series expansion of the Berkooz-Douglas matrix model, which describes the D/D4 brane system. At high temperature the model is weakly coupled and we develop the series to second order. We check our results against the hightemperature regime of the bosonic model without fermions) and find excellent agreement. We track the temperature dependence of the bosonic model and find backreaction of the fundamental fields lifts the zero-temperature adjoint mass degeneracy. In the low-temperature phase the system is well described by a gaussian model with three masses m t A = 1.964±.3, m l A = 2.1 ±.3 and m f = ±.1, the adjoint longitudinal and transverse masses and the mass of the fundamental fields respectively.

2 Contents 1 Introduction 1 2 he Berkooz-Douglas Model 4 3 High-emperature Expansion Leading order ext-leading order 11 4 High-temperature coefficients from numerical simulations 14 5 he Bosonic Berkooz-Douglas model 17 6 Conclusions 2 A ables for the ω s. 22 B he high-temperature behaviour of energy E, Polyakov loop P, R 2 and mass susceptibility C m for the supersymmetric model Introduction he Berkooz-Douglas model BD model) [1] was introduced as a non-perturbative formulation of M-theory in the presence of a background of longitudinal M5-branes with the M2-brane quantised in light-cone gauge. Its action is written as that of the BFSS model [3] with additional fundamental hypermultiplets to describe the M5-branes. he BFSS model can also be viewed as a many-body system of D-branes of the IIA superstring. In this framework the BD model is a D/D4 system with the massless case being the D/D4 intersection. When the number of D-branes far exceeds that of the D4-branes the dynamics of the D-branes is only weakly affected by that of the D4-branes and is captured by the IIA supergravity background holographically dual to the BFSS model. In this context the D4-branes, representing the fundamental fields of the BD model, are treated as Born-Infeld probe 4-branes. his holographic set up is a tractable realisation of gauge/gravity duality with flavour. Both the BFSS model and the BD model are supersymmetric quantum mechanical models with an SU) gauge symmetry. When they are put in a thermal bath they become strongly coupled at low temperature. At finite temperature their gravity duals involve a black hole whose Hawking-temperature is that of the thermal bath. hese duals can be used to provide non-perturbative predictions at low temperature. he BFSS and BD models can also be 1

3 studied by the standard non-perturbative field theory method of Monte Carlo simulation. hese models therefore provide excellent candidates for testing gauge/gravity duality nonperturbatively and in a broken supersymmetric setting. here are now several non-perturbative studies of the BFSS model [7 11] and several recent reviews[12 14]. Also, the BD model was recently studied non-perturbatively in [15]. In all cases the predictions from the gauge/gravity duals were found to be in excellent agreement with that of Monte Carlo simulations of the finite-temperature models. he situation is conceptually simpler at high temperature as the dimensionless inverse temperature, scaled in terms of the BD-coupling, provides a natural small parameter for the model. In this paper, we obtain the first two terms in the high-temperature expansion of the BD model. In the high-temperature limit only the bosonic Matsubara zero modes survive and the resulting model is a pure potential. his potential, which provides the non-perturbative aspect of our high-temperature study, also plays a role in the ADHM construction [16]. We study the model for adjoint matrix size between 4 and 32 for f = 1 with f the number of D4-branes) and for f between 2 and 16 for from 9 to 2. For f 2 we find that the system has difficulties with ergodicity. In particular, for f = 2 and f = 2 +1 the system failed to thermalise satisfactorily. In contrast the system has no difficulties for f = 2 1. his condition is closely related to the singularity structure of instanton moduli space, where irreducible SU f ) instantons of Chern number exist only for f 2 [17, 18]. he moduli space of such instantons is equivalent to the zero locus of the potential with X a = and D A = see equation 2.4)). his moduli space is in general singular and non-singular only when this bound is satisfied. here is also a natural dimensional analogue of the BD model, which has = 4 supersymmetry, associated with the D1/D5 system of [5], whose BFSS relative was discussed in [19 21]. When the Euclidean finite-temperature version of this dimensional quantum field theory is considered on a torus with the spatial circle of period β and euclidean time 1 of period 1/, then at high temperature the fermions decouple and one is left with the purely bosonic version of the BD model. We refer to this model as the bosonic BD model and study the small period behaviour equivalent for us to our high-temperature regime) of the massless version of this model as a check on our high-temperature series. We find the hightemperature series results are in excellent agreement with Monte Carlo simulations of the bosonic BD model. By fitting the dependence, of the expectation values of our observables, on the number of flavour multiplets, f, we find that extrapolation, to f =, agrees well with the corresponding observables of the BFSS model. As β, the inverse temperature, grows the bosonic BD model undergoes a set of phase transitions. hese are the phase transitions of the bosonic BFSS model. We find the hightemperature series expansion is valid up to β 1/2, which is just below the phase transition region. Above the transition the bosonic BD model is well described by free massive fields, 1 In this paragraph we avoid using β for 1/ for simplicity of the comparison. 2

4 where the backreaction of the fundamental fields has lifted the degeneracy of the longitudinal and transverse masses. he principal results of this paper are: We obtain expansions for observables of the BD model to second order in a hightemperature series. We tabulate the coefficients of this expansion as functions of and f in the range 4 32 and 1 f 16. We measure the expectation values of the composite operator r 2 bos, see equation 2.8)), and the mass susceptibility C m bos, see equation 5.7)), of the bosonic BD model as a function of temperature down to zero and use it to check our coefficients for the high-temperature series of the full BD model. We find that the fundamental fields of the bosonic BD model have mass m f = ±.1. We measure the backreacted mass of the longitudinal adjoint scalars to be m l A = 2.1±.3 and find that the transverse mass is largely unaffected by backreaction being m t A = ±.3, which should be compared with the bosonic BFSS model, where the fields have mass m A = ±.7. We use the measured masses to predict the zero-temperature values of our fundamental field observables r 2 bos and mass susceptibility C m bos and find excellent agreement with direct measurements. he paper is organised as follows: In section 2 we present the finite-temperature BD model and describe our notation and observables. In section 3 we set up and implement the high-temperature series expansion working to second order in the inverse temperature β. Section 4 describes the dependence of our observables on the coefficients in the expansion, which must be determined by numerical simulation of the zero-mode model. In Section 5 we perform lattice simulations of the bosonic BD model and find excellent agreement with the high-temperature expansion. We also find the low-temperature phase of the model is well described by a system of gaussian quantum fields. Section 6 gives our concluding remarks. here are two appendices; appendix A gives tables, for different and f, of the coefficients determined non-perturbatively while appendix B presents graphs of predictions for the hightemperature behaviour of our observables for the supersymmetric model. 3

5 2 he Berkooz-Douglas Model We begin by describing the field content of the model following the notation used in [5]. he action of the BFSS model is given by S BFSS = 1 g 2 dt 9 i=1 { 1 r 2 D X i ) [Xi, X j ] 2 i 2 Ψ C 1 Γ D Ψ + 1 } 2 Ψ C 1 Γ i [X i, Ψ] 2.1) where D = t i[a, ], Ψ is a thirty-two component Majorana Weyl spinor, Γ µ are ten dimensional gamma matrices and C 1 is the charge conjugation matrix satisfying C 1 Γ µ C1 1 = Γ µ. he fields X i and Ψ are in the adjoint representation of the gauge symmetry group SU) and A is the gauge field. o describe the addition of the fundamental fields we break the SO9) vector X i into an SO5) vector X a and an SO4) vector which we re-express as X ρ ρ via 2 X ρ ρ = i 2 4 σ m ρ ρx 1 m, 2.2) m=1 where σ 4 = i1 2 and σ A s A = 1, 2, 3) are the Pauli matrices. he X ρ ρ ρ, ρ = 1, 2) are complex scalars which together transform as a real vector of SO4) which satisfies the reality condition X ρ ρ = ε ρσ ε ρ σ Xσ σ. he indices ρ and ρ are those of SU2) R and SU2) L, respectively, where SO4) = SU2) L SU2) R. he nine BFSS scalar fields, X i, become X a a = 1,, 5) and X ρ ρ. he sixteen adjoint fermions of the BFSS model become λ ρ and θ ρ with λ ρ being SO5, 1) symplectic Majorana-Weyl spinors of positive chirality and satisfying λ ρ = ε ρσ λ c ) σ while θ ρ are symplectic Majorana-Weyl spinors of negative chirality satisfying θ ρ = ε ρ σ θ c ) σ. hey combine together to form an SO9, 1) Majorana-Weyl spinor in the adjoint of SU). his SO9) symmetry is recovered only if the fundamental fields are turned off. o describe the longitudinal M5-branes or D4-branes), we have Φ ρ and χ, which transform in the fundamental representations of both SU) and the global SU f ) flavour symmetry. Φ ρ are complex scalar fields with hermitian conjugates Φ ρ, and χ is an SO5, 1) spinor of negative chirality. After rotating to imaginary time the Euclidean action describing the model at finite temperature = β 1 becomes: [ β 1 S = dτ r 2 D τ X a D τ X a D τ X ρ ρ D τ X ρ ρ λ ρ D τ λ ρ + 1 ) 2 θ ρ D τ θ ρ ) + tr D Φρ τ D τ Φ ρ + χ D τ χ 1 r 4 [Xa, X b ] ) 2 [Xa, X ρ ρ ][X a, X ρ ρ ] 2 Here X 8 of [5] is replaced by X 8., 4

6 r D A D A + tr Φρ X a m a ) 2 ) Φ ρ r A=1 1 2 λ ρ γ a [X a, λ ρ ] θ ρ γ a [X a, θ ρ ] ) 2iε ρσ θ ρ [X σ ρ, λ ρ ] tr χ γ a X a m a )χ + 2iε ρσ χ λ ρ Φ σ + 2iε Φρ ρσ λ χ) ] σ, 2.3) where D A = σ A σ ρ ) 1 2 [ X ρ ρ, X σ ρ ] Φ σ Φρ, 2.4) with D τ the covariant derivative which, for the fields of the fundamental multiplet, Φ ρ and χ, acts as D τ = τ ia). he trace of SU) is written as r while that of SU f ) is denoted by tr. he diagonal matrices, m a, correspond to the transverse positions of the D4-branes. We fix the static gauge: τ A =, so the path integral requires the corresponding ghost fields c and c with the ghost term β dτ r τ cd τ c added to the action 2.3). We will restrict our attention to m a = so that the D4-branes are attached to the D-branes, and the strings between D and D4 are massless, i.e. the fundamental fields are massless. he factor of in front of the integral in 2.3) is the remnant of the t Hooft coupling λ = g 2 which is kept fixed and absorbed into τ and the fields with β = λ 1/3 /. ote that without loss of generality we can set λ = 1. As discussed in the introduction, the BFSS model is the matrix regularization of a supermembrane theory [2], so the BFSS part of this model can be also interpreted as M2-brane dynamics. In this context the D4-branes lift to M5-branes and the model can describe M2- branes ending on longitudinal M5-branes. he BD model is a version of supersymmetric quantum mechanics and could in principle be treated by Hamiltonian methods. he partition function is then Z = re βh ) = [dx][dλ][dθ][dφ][d Φ][dχ][dχ ][da]e S 2.5) with r the trace over the Hilbert space of the Hamiltonian restricted to its gauge invariant subspace and the action S, in the path integral, is given by equation 2.3). he measure in the path integral for the partition function 2.5) has a hidden dependence on temperature due to the presence of the Van Vleck-Morette determinant [22] in the definition of the path integral measure. his determinant arises from the kinetic contribution to the action 2.3) which, as written, is temperature dependent. o remove the temperature dependence from the measure, we rescale the variables in the original action 2.3) so that the kinetic terms, including the gauge potential, are independent of β. For this τ βτ, X i β 1 2 X i, Φ ρ β 1 2 Φ ρ, A β 1 A, c β 1 2 c and c β 1 2 c. he fermions do not need rescaling. he path integral measure is now temperature independent and, when the mass 5

7 is zero, the only temperature dependence is β 3 for the bosonic potential and β 3/2 for the fermionic potential. If the mass term is included it enters as βma in the potential with the overall scales of β 3 and β 3/2 in the bosonic and fermionic contributions respectively. he temperature dependence of the model is now explicit. he principal observable of the model is the energy 3, E = H / 2. Once the temperature dependence of the model has been made explicit, as described above, one can then simply note that 2 E is minus the derivative of logarithm of the partition function with respect to β, returning to the original variables one readily sees that in the path integral formulation: E = ε b + ε f, where ε b = 3 β [ dτ r 14 ) β [Xi, X j ] 2 + tr Φ ρ X a 2 Φ ρ Φ ρ [ X σ ρ, X ρ ρ ]Φ σ 12 ) ] Φ ρ Φ Φσ σ Φ ρ + Φ ρ Φ Φσ ρ Φ σ, ε f = 3 β [ 1 dτ r 2β 2 λ ρ γ a [X a, λ ρ ] 1 2 θ ρ γ a [X a, θ ρ ] + ) 2iε ρσ θ ρ [X σ ρ, λ ρ ] + tr χ γ a X a χ 2iε ρσ χ λ ρ Φ σ 2iε Φρ ρσ λ χ) ] σ. 2.6) We see only the potential contributes and the coefficients 3 and 3/2 of the bosonic and fermionic terms arise from the differentiation. As in [6], there are two other interesting observables: R 2 = 1 β β dτ r X i 2, P = 1 r exp [iβa]). 2.7) Here R 2 is a hermitian operator whose expectation value is a measure of the extent of the eigenvalue distribution of the scalars X i and P is the Polyakov loop. ote: Path-ordering is not needed here for the Polyakov loop as we consider A in the static gauge. Since the model has new degrees of freedom it is important to consider other observables that capture properties of these new fields. he natural candidates are r 2 = 1 β f β β 3 2 dτ tr Φ ρ Φ ρ, 2.8) which is the analogue of R 2 for the fundamental degrees of freedom, and the condensate defined as c a m) = m a 1 ) 1 β { β log Z = dτ tr 2 Φ ρ m a X a )Φ ρ + χ γ χ} a. 2.9) β However, for us, with m a =, c a will be zero. So our focus will be on the mass susceptibility 3 We divide by 2 so that E remains finite in the large- limit. C m := ca ), 2.1) ma 6

8 i.e. the derivative with respect to m a with a fixed not summed over) and evaluated at m a = where C m = 2 β β dτ tr Φ ρ Φ ρ 5β β { } ) 2 dτ tr 2 Φ ρ X a Φ ρ + χ γ a χ. 2.11) Here a in 2.11) is summed over a = 1,, 5 and the same applies hereinafter. 3 High-emperature Expansion In this section, we develop the high-temperature expansion of the BD model. For very high temperatures only the Matsubara zero modes, i.e. the zero modes in a Fourier expansion, survive and the model reduces to a bosonic matrix model for these modes. A high-temperature series expansion is therefore obtained by developing a perturbative expansion of the model in the non-zero modes. he zero-modes must then be treated non-perturbatively and this is done by Monte Carlo simulation. Our strategy is therefore to expand the model in Matsubara modes, show that the temperature can be seen as a coupling constant for these modes and then integrate out the non-zero modes order by order in perturbation theory to obtain an effective action for the zero modes, which can then be treated non-perturbatively. o obtain the series to second order we will only need one loop computations. he nonzero mode integration can be done analytically and yields an effective action and observables in terms of the zero temperature variables. As a final step the integration over these zero modes must then be performed non-perturbatively via Monte Carlo simulation. he Fourier expansion of the fields is given by X i τ) = Xne i 2πinτ/β, λ ρ τ) = λ rρ e 2πirτ/β, θ ρ τ) = θ r ρ e 2πirτ/β, n Z cτ) = n Z,n Φ ρ τ) = n Z r Z+ 1 2 c n e 2πinτ/β, cτ) = n Z,n Φ nρ e 2πinτ/β, χτ) = r Z+ 1 2 c n e 2πinτ/β, r Z+ 1 2 χ r e 2πirτ/β, 3.1) where thermal boundary conditions require that the bosons and ghosts are periodic in τ while the fermions are anti-periodic. he action 2.3) now takes the form of the sum of a zero mode action, the kinetic term for the non-zero modes and an interaction term. As discussed above, in the high-temperature limit, only the zero-modes play a role. As the temperature is lowered one can integrate out the non-zero modes perturbatively with β playing the role of a perturbation parameter. Using this procedure, the first two terms in the high-temperature expansion of E, R 2 and P for the BFSS model were obtained in [6]. We follow the same method here and obtain the corresponding expansion of these observables for the BD model and for them the novel feature 7

9 will be the additional dependence on f, the number of flavour multiplets. In addition we have the new observable r 2 and C m. In order to develop the high-temperature series it is convenient to rescale the scalar fields in 2.3) as follows X i β 1 4 X i, A β 1 4 A, Φ β 1 4 Φ, X i n β 1 2 X i n, Φ n β 1 2 Φn, c n β 1 2 cn, c n β 1 2 cn, 3.2) while the fermions remain unchanged. his rescaling makes the coefficients of the zero-mode terms and the kinetic terms independent of β so that one can concentrate on the β-dependence, which now appears only in the interaction terms. he rescaling 3.2) can be understood as the rescaling of section 2, which was necessary to remove the temperature dependence of the measure, followed by the further zero-mode rescaling X i β 3 4 X i, Φ ρ β 3 4 Φ ρ and A β 3 4 A. his zero-mode rescaling means the partition function becomes Z = β )+4 f ) Z, 3.3) where Z is the partition function in terms of the rescaled fields of 3.2) and the only remaining temperature dependence is in S int. We can now develop the high-temperature series by diagrammatic techniques with β playing the role of a coupling. he action is then written in terms of the variables of 3.2), which we will use for the remainder of the paper, as S = S + S kin + S int, 3.4) where S is a zero-mode action S = 4 r [X i, X j ]2 + 2[A, X i ] 2) + tr S kin is the kinetic part of the action for non-zero modes [ r X nx a n a + S kin = n 2πn) r 2πir 2 Φ ρ A2 Φ ρ + Φ ρ Xa ) 2 Φ ρ Φ ρ σ ρ [ X, X ρ ρ]φ σ 1 2 Φ ρ Φ σ Φ σ Φ ρ + Φ ρ Φ ρ Φ σ Φ σ ), 3.5) ) ρ ρ X n X nρ ρ + 2 c n c n + tr 2 Φ ρ n Φ ) ] nρ [ ) )] r λ ρ r λ rρ + θ ρ r θ r ρ + tr 2χ r χ r, 3.6) and S int is the interaction part of the action. he terms quadratic in non-zero modes present in S int but not present in the BFSS model are where S int = β 3 4 V A) 1 + V B) 1 ) β 3 2 V A) 2 + V B) 2 + V 3 ) + Oβ 9 4 ), 3.7) V A) 1 = 4π n tr Φ ρ n AΦ nρ) + n r triχ r Aχ r), 8

10 V B) 1 = r trχ r γa X a χ r + 2iε ρσ χ r λ rρφ σ + 2iε ρσ Φρ λ σ r χ r), V A) 2 = n tr Φ ρ n A2 Φ nρ ), V B) 2 = tr Φρ n Xa 2 Φ nρ + Φ ρ Xa nxnφ a ρ Φ ρ σ ρ n [ X, X ρ ρ]φ nσ Φ ρ σ ρ [ X n, X nρ ρ]φ σ n Φ ρ n Φ σ Φ σ Φ nρ Φ ρ n Φ nσ Φ σ Φ ρ + 2 Φ ρ n Φ ρ Φ σ Φ nσ + 2 Φ ρ n Φ nρ Φ σ ) Φ σ, V 3 = [ tr Φ ρ n Xa nx a Φ ρ + Φ ρ n Xa XnΦ a ρ + Φ ρ Xa nx a Φ nρ + Φ ρ Xa X nφ a nρ ) n Φ ρ σ ρ n [ X n, X ρ ρ ]Φ σ + Φ ρ σ ρ n [ X, X nρ ρ]φ σ + Φ ρ σ ρ [ X n, X ρ ρ]φ nσ + Φ ρ σ ρ [ X, X nρ ρ]φ nσ ) 1 2 Φ ρ n Φ σ Φ σ nφ ρ + Φ ρ Φ nσ Φ σ Φ nρ ) + Φ ρ n Φ ρ Φ σ nφ σ + Φ ρ Φ nρ Φ σ Φ nσ ) χ r γa X nχ a r+n + 2iε ρσ χ r λ r+n ρφ nσ + ] 2iε ρσ Φρ n λ σ r χ r+n). 3.8) r V 3 does not contribute to the expectation values of operators at next-leading order. wo such vertices would be required and the resultant contribution would therefore be of higher order in β. Similarly, fermionic terms that involve only non-zero modes also scale as β 3 2, and again contribute at a two and higher loop order to the expectation values of observables. he zero-mode action 3.5) corresponds to the bosonic part of the original model 2.3) dimensionally reduced to a point and plays an important role in the ADHM construction as the solutions to S = with D A =, where D A is given in 2.4), provide the ADHM data[16]. his zero-mode model is the flavoured bosonic version of the IKK model[4]. We use the notation DR for the expectation value calculated with this dimensionally reduced model. hus for a generic observable, O which is a function of X i, A and Φ ρ we denote O DR = 1 Z [dx ][d Φ ][dφ ]Oe S. 3.9) Furthermore we denote AB DR,c = AB DR A DR B DR, 3.1) and the subscript c denotes connected part. Identities such as d dx X Oe S ) = for X and similar identities for A and Φ ρ yield the Ward-type identities: 9 2 1) O DR λ ds λx ) dλ O + λ doλx ) λ=1 dλ =, DR λ=1 DR 2 1) O DR λ ds λa) dλ O + λ doλa) λ=1 dλ =, DR λ=1 DR 2 f O DR λ ds λφ ρ ) dλ O + λ doλφ ρ) λ=1 dλ =. 3.11) λ=1 DR DR 9

11 hese identities can be used to simplify various expressions and in particular one can see that one never needs to consider the insertion of S or s = r 1 4 [Xi, X j ]2) + tr Φρ Xa 2 Φ ρ Φ ρ σ ρ [ X, X ρ ρ]φ σ 1 2 Φ ρ Φ σ Φ σ Φ ρ + Φ ρ Φ ρ Φ ) σ Φ σ 3.12) with other correlators as they can be eliminated by use of these identities. identities resulting from 3.11) are that he simplest 4 S DR = 1 2 1) + 4 f and 4 s DR = 8 2 1) + 4 f, 3.13) the latter of which establishes the equivalent leading order expression for the energy using 2.6) as discussed in the next section. We will next present the leading high-temperature expansion of our observables. 3.1 Leading order he expectation values of our observables to leading order are determined solely by the zero modes. o this order the partition function is given by 3.3) with Z a constant. We therefore have E = 3 {8 4 β ) } f + Oβ 3 2 ). 3.14) he direct expression using 2.6) gives E = 3 2 β 1 s DR + Oβ 3/2 ) and using the second identity of 3.13) we see that this agrees with 3.14). Also, the leading terms in the β-expansion of R 2 and the expectation value of the Polyakov loop are R 2 = β r Xi + Oβ), P = 1 β DR 1 r A2 + Oβ 3 ). 3.15) DR For our new observables r 2 and C m we have the leading contributions r 2 = β tr Φ ρ Φ ρ + Oβ) 3.16) f DR and C m = 2β 1 2 tr Φρ Φ ρ 2 DR 5 tr Φρ Xa Φ ρ ) 2 ) + Oβ). 3.17) DR,c ote: All the leading order contributions are purely bosonic, since fermions decouple at high temperature. he necessary expectation values are computed numerically via Monte Carlo simulation with the action S of equation 3.5) and given in the tables in appendix A for different values of and f. 1

12 3.2 ext-leading order he higher order contributions in the high-temperature expansion come from integrating out the non-zero modes in 3.4). he first subleading order is obtained by performing the gaussian integrals over the non-zero modes, where the potential is truncated as in 3.7), and expanding the resulting exponential and ratio of determinants in terms of β. Examining the action 2.3) we see the fermionic terms can be written in the form dτ [ 1 r + r 2 λ ρ D τ + γ a ad X a )λ ρ θ ρ D τ γ a ad X a )θ ρ ) θ ρ J ρ + tr ) ) + tr χ D τ γ a X a )χ χ J + J χ) ], 3.18) and the commutator action of X a is denoted by ad X a. Since J ρ and J are fermionic currents that depend linearly on λ σ, integrating out θ ρ and χ gives the additional contributions 1 2 J ρ G θ J ρ and J G χ J to the quadratic form for λ ρ. Here G λ, G θ and G χ are Green s functions for λ ρ, θ ρ and χ, respectively. hese current-current terms will be of order β 3/2 and contribute at the sub-leading order under consideration. he non-zero modes can now be integrated out and to one loop we obtain S eff = S ln Det[1 + β3/4 ad A 2πn )2 + β 3/2 M X ]Det[1 + β 3/4 A 2πn) 2 2πn )2 + β 3/2 Det[1 + ad A n 2πn ] [ { )} ln Pf ɛ 1 + β 3/4 i ad A + γa ad X a 2πin + 1 n= 2 { ) )} with + n= Pf ɛ 1 + β 3/4 i ad A γa ad X a 2πin ) Det M φ ] 2πn) 2 1 β 3/4 ia + γa X a 2πin ) β 3/2 r[g n λ adx )G n θ ad X )] + 2r[G n λ Φ G n χ Φ ] ), 3.19) G n θ = Gn λ = Gn χ = 1 2πin ). Equations 3.7) with details in 3.8) specify the quadratic forms whose determinants and Pfaffians enter in 3.19), and here r is the operator trace. In detail one has: M X X) a n = 1 2 [Xb, [X, b Xn]] a [Xb, [X a, Xn]] b + 1 ρ ρ [ X 2, [X ρ ρ, Xn]] a Φ ρ Φρ Xa n, M X X) nρ ρ = 1 2 [X σ ρ, [ [Φ ρ Φσ, X nσ ρ ], M Φ Φ) nρ = X a 2 σ ρ Φ nρ + [ X, X ρ ρ]φ nσ σ σ X, X nρ σ ]] [X σ σ ρ σ, [ X, X nσ ρ ]] 1 2 [Xa, [X a, X nρ ρ ]] ) ] 11

13 + Φ σ Φσ Φ nρ + Φ nσ Φσ Φ ρ 2Φ ρ Φσ Φ nσ 2Φ mρ Φσ Φ σ. 3.2) Expanding 3.19) in β we will obtain S eff with S eff = S + β 3/2 S ) Let us look at the individual contributions in more detail. After straightforward computations we find the contribution from the first determinant in 3.19) due to the integration over non-zero X a modes gives S X 1 = 24 he ghost contribution similarly expanded gives and the contribution from the fundamental scalars is 16 r X i 2 18 r A 2) ) tr 12 Φ ρ Φ ρ. 3.22) S g 1 = 2 24 r A2, 3.23) S1 Φ = 1 2f rx a 2 A 2 ) + 3 tr 12 Φ ρ Φ ρ). 3.24) Putting the bosonic contributions together we have S1 bos = 2 { r X i 2 r A ) tr Φ ρ Φ ρ } + f 6 he fermionic contributions can similarly be evaluated to give r X a 2 r A 2). 3.25) S θ 1 = S λ 1 = r X a 2 r A 2 ) 3.26) for the Pfaffian contribution to the integration over θ ρ and λ ρ. he θ current-current contribution gives S J ρ G θ J ρ ρ ρ 1 = 2 r X X ρ ρ). 3.27) hese three contributions together recombine to give an SO9) invariant term, which is the fermionic part of SBF SS 1. ext considering the χ determinant we find S χ 1 = f 2 and its current-current contribution is r Xa 2 r A 2 ), 3.28) S J G χj 1 = 2 2 1) tr Φ ρ Φ ρ. 3.29) Putting all these fermionic contributions together we find { S fer 1 = 2 r X i 2 r A ) } 2 tr Φ ρ Φ ρ f 2 r X a 2 r A 2). 3.3) 12

14 Finally, defining e β3/2 S 1 = 1 + β 3/2 O, so for the supersymmetric model adding 3.25) and 3.3) we have O = 4 { 3 r X i 2 r A ) } 16 2 tr Φ ρ Φ ρ f r X a 2 r A 2). 3.31) Similarly using 3.25) we can define O bos = S bos 1, which we can write O bos = 1 2 O 3 8 tr Φ ρ Φ ρ. 3.32) hese expressions will be useful in section 4. he partition function including the next to leading corrections is now given by Z = β )+4 f ) 1 + β 3/2 O DR ) Z, 3.33) and the temperature dependence is explicit. We immediately have the subleading correction to the energy E by straightforward differentiation of 3.33) to obtain that E = 3 {8 4 β ) } f β1/2 O DR. 3.34) urning to the high-temperature behaviour of R 2 and the Polyakov loop, the resulting expectation values are given by R 2 = β r Xi + β 1 1 ) ) r Xi )O + Oβ 5 2 ), 3.35) DR DR,c and P = 1 β 3 2 [ { 1 1 r A2 β r A4 2 4! DR DR 1 } ] 1 2 r A2 )O + Oβ 3 ). DR,c 3.36) he constant ) is the contribution due to the expectation value of the non-zero 2 modes, which are traceless. Our observable r 2 is similarly given by r 2 = β tr Φ ρ Φ ρ f DR + β tr Φρ Φ ) ρ O f DR,c ) + Oβ 5 2 ), 3.37) and its bosonic version is again obtained by replacing O with O bos. In terms of Fourier modes we have c a m) = tr 2β 1 2 m a Φρ Φ ρ + 2βm a Φ ρ n Φ nρ + χ rγ a χ r 2β 3 4 Φρ Xa Φ ρ n r 13

15 2β 3 4 Φ ρ n Xa Φ nρ + Φ ρ Xa nφ nρ + Φ ) ρ n Xa nφ ρ ) 2β 3 2 Φ ρ n Xa n mφ mρ. n n,m 3.38) However, we will restrict ourselves to the massless case and as discussed SO5) invariance guarantees that this observable is zero so we focus on the mass susceptibility, C m. Calculating C m in the high-temperature expansion to the next to leading order yields C m = 2β 1 2 tr Φ ρ Φ ρ DR 25 ) tr Φ ρ Xa Φ ρ ) 2 DR,c + 2β f 3 + tr Φ ρ Φ ρ)o 2 DR,c 5 tr Φρ Xa Φ ρ ) 2 O ) + Oβ 5 2 DR,c ). 3.39) he contribution f /3 in the second parentheses in 3.39) contains both bosonic and fermionic contributions with the fermionic contribution being f /2, while the bosonic contribution is f /6. herefore, the condensate susceptibility for the bosonic model is obtained from 3.39) by replacing the numerical constant f 3 by f 6 and O by O bos. he resulting expression is C m bos = 2β 1 2 tr Φ ρ Φ ρ DR 25 ) tr Φ ρ Xa Φ ρ ) 2 DR,c f + 2β 6 + tr Φ ρ Φ ρ)o bos DR,c 2 5 ) tr Φρ Xa Φ ρ ) 2 O bos + Oβ DR,c 5 2 ). 3.4) An alternative to the above treatment is to work directly with perturbation theory in β, but we believe the structure of the computations is simpler in the above treatment. he contributions to O, E, R 2 and P from the pure BFSS model were derived in [6] and when f and the fundamental fields are set to zero our results reproduce the results reported there. 4 High-temperature coefficients from numerical simulations In this section we express the coefficients, Ξ i,used in the high-temperature expansion of observables see equations 4.3) and 4.4)), in terms of the primitive observables ω i, ω i,j and ω i,j,k i, j, k = 1,..., 6) defined in equations 4.2) below. We have the following expressions: Ξ 1 = 8ω 1 + ω 3 ) 2 f Ξ 2 = 5ω 1 + 4ω 3, Ξ 3 = f 3 ) f { ω ) } 16 2 ω 4, 5ω 1,1 + 15ω 1,2 + 16ω 1,3 ) ω 1,3 + 4ω 3,3 ) ) 5ω 1,4 + 4ω 3,4 ), 14

16 Ξ 4 = 1 2 ω 1, Ξ 5 = 1 4! ω { 4 + f 3 Ξ 6 = 2 f ω 4 2ω 5,5 ), Ξ 7 = 2 f f Ξ 8 = 4 f { 4 + f 3 Ξ 9 = ω 4, Ξ 1 = 1 4! ω 6, Ξ bos 1 = 1 2 Ξ 1 3 f Ξ bos 3 = f 3 Ξ bos 5 = 1 4! ω ω 4, 2 Ξ bos 7 = f 3 4 f { 4 + f 8 = 2 f 3 Ξ bos { 4 + f 3 5ω 1,2 ω 1,1 ) ω 1,3 + 4 f 1 19 ) } ω 1,4, ω 1, ω 3, ω 2,5,5 + 3ω 1,5,5 ) ω 3,5, ) } 16 2 ω 4,4, ) ω 4,5,5 } ) 4 + f 5ω 1,1 + 15ω 1,2 + 16ω 1,3 ) ω 1,3 + 4ω 3,3 ) 1 5 ) 5ω 1,4 + 4ω 3,4 ), 8 2 { 4 + f 3 { 4 + f 3 5ω 1,2 ω 1,1 ) ω 1,3 + 4 f 3 ω 1, ω 3, ω 2,5,5 + 3ω 1,5,5 ) ω 3,5, ) } ω 4,4, ) ω 1,4 } 1 5 ) } 8 2 ω 4,5,5, 4.1),, where r ω 1 = 1 r A 2 DR, ω 2 = 1 rx a 5 ) 2 = ω DR 1, ω 3 = 1 4 ω 4 = 1 tr Φρ Φ ρ, ω DR 5 = 1 tr Φρ f X1 Φ ρ DR =, ω 6 = 1 ρ ρ X X ρ ρ, DR r A 4, DR ω 1,1 = ω 1,3 = 1 4 r A 2 ) 2, ω 1,2 = 1 DR,c 5 r A 2 r r r A 2 rx a ) 2 DR,c, ρ ρ X X ρ ρ, ω 1,4 = r A 2 tr DR,c Φ ρ Φ ρ f ) 2 DR,c, ω 3,3 = 1 ρ ρ X 16 X ρ ρ, ω 3,4 = ρ ρ r X 4 X ρ ρ tr Φ σ Φ σ DR,c f ω 4,4 = tr Φρ Φ ) 2 ρ, ω 5,5 = tr ) 2 Φρ f DR,c X1 Φ ρ, f DR,c DR,c, 15

17 ω 1,5,5 = 2 r A 2 tr Φ ρ ) 2 X1 Φ ρ, f DR,c ω 2,5,5 = 2 rx 1 ) 2 tr Φ ρ ) 2 X1 Φ ρ, f DR,c ω 3,5,5 = 2 ρ ρ r X 4 X ρ ρ tr Φσ X 1 ) 2 Φ σ f ω 4,5,5 = 2 tr Φρ Φ ρtr Φ σ X 1 ) 2 Φ σ f DR,c DR,c,. 4.2) In terms of the Ξ i the observables of the full BD model 2.3) become: E = 3 {8 4 β ) } f + β 1 2 Ξ1 + Oβ 2 ), R 2 = β 1 2 Ξ2 + βξ 3 + Oβ 5 2 ), [ P = 1 β 3 2 For the bosonic BD model we have Ξ 4 β 3 2 Ξ5 + Oβ 3 ) r 2 = β Ξ9 + β 2 + Ξ ) 7 + Oβ 5 2 ), 2 f ], C m = β 1 2 Ξ6 + βξ 7 + Ξ 8 ) + Oβ 5 2 ). 4.3) E bos = 3 4 β 1 {8 R 2 bos = β 1 2 Ξ2 + βξ bos 3 + Oβ 5 2 ), [ P bos = 1 β ) } f + β 1 2 Ξ bos 1 + Oβ 2 ), Ξ 4 β 3 2 Ξ bos 5 + Oβ 3 ) r 2 bos = β 1 2 Ξ9 + β 1 2 f Ξ bos 7 + Oβ 5 2 ), C m bos = β 1 2 Ξ6 + βξ bos 7 + Ξ bos 8 ) + Oβ 5 2 ). 4.4) he observables of interest for the high-temperature expansion are all expressed in terms of Ξ i and Ξ bos i listed above. As discussed they are temperature independent and depend only on, the matrix dimension of the BFSS fields and f, the number of flavour multiplets. We computed their values for a range of and f by hybrid Monte Carlo simulation with the action S given in 3.5). We tabulate our results for different and f. We choose = 4, 6, 8, 9, 1, 12, 14, 16, 18, 2, 32 for f = 1 and tabulate ω s, the building blocks of Ξ i, in able 1. From the results of able 1 we extrapolate the -dependence of the ω s by fitting them with a function 4, a + b/ + c/ 2 see Figure 3 and 4). he limiting extrapolated values are 4 ote that as expected we find it necessary to include a linear fall of in 1/ for large. his is in contrast to the BFSS model, where the fall off is 1/ 2 ], 16

18 included as the row = in able 1. Ξ i=1,,5,and1 naturally reduce to counterparts in the BFSS model when the fundamental fields are removed. We extrapolate Ξ i=1,,5,and1 for f = 1, 2, 4, 6, 8, 1, 12, 14, 16 for fixed = 12, 14, 16, 18 and 2 to f = and find good agreement, to within the quoted errors, with the measured values for their BFSS counterparts as quoted in [6]. Figure 5 shows plots of the ω s against f for each and we fit the dependence on f basically with a quartic polynomial. However, we find that higher order terms contribute for some ω s and by using the fitting function a + b f + ce d f we capture the dependence on f over the range considered. 5 he Bosonic Berkooz-Douglas model We are in the process of making a direct comparison of both the high-temperature regime of the BD model as determined by the above predictions and the low-temperature regime as predicted by gauge/gravity with results from a rational hybrid Monte Carlo simulation using the code used in [15]. We will present those results in a separate paper as, apart from their value as a check on the code and the computations presented here, they have additional physics that merits a separate discussion. For this paper we restrict our considerations to a comparison of the results obtained here with those obtained from the bosonic Berkooz-Douglas model, whose Euclidean action is given by S bos = β dτ [ 1 r 2 D τ X a D τ X a D τ X ρ ρ D τ X ρ ρ 1 4 [Xa, X b ] ) 2 [Xa, X ρ ρ ][X a, X ρ ρ ] + tr D τ Φρ D τ Φ ρ + Φ ρ X a m a ) 2 Φ ρ ) r 3 A=1 D A D A ]. 5.1) Our comparison is presented in Figure 1, where we restrict our considerations to a high precision test with = 1 and f = 1. As one can see from the figure the agreement is excellent. Furthermore, the high expansion remains valid at temperatures as low as 1.. Below this temperature the figure shows evidence of a phase transition. his is the phase transition of the bosonic BFSS model. From studies of the bosonic BFSS model [1, 2, 21] we know that it undergoes two phase transitions at c2 =.95 ±.2 and c1 =.8761 ±.3. hese are driven by the gauge field A, which at high temperature behaves as one of the X i, while at low temperature it effectively disappears from the system and can be gauged away at zero temperature. As the temperature is increased through c1 there is a deconfining phase transition with the symmetry At) At) + α1 broken and the distribution of eigenvalues of the holonomy 5 becomes non-uniform. When the temperature reaches c2 the spectrum of the holonomy becomes gapped and above this temperature the eigenvalues no longer cover the entire [, 2π] 5 he Polyakov loop, P = 1 ru), where U is the holonomy. 17

19 3 r 2 bos c m Figure 1: Comparison of the high-temperature predictions for the fundamental observable r 2 bos and the derivative of the condensate at zero mass, c/ m) = C m bos, with a Monte Carlo simulation of the bosonic BD model. he simulation is for f = 1 and = 1. range. In the low-temperature phase the bosonic BFSS model becomes a set of massive gaussian matrix models with Euclidean action S eff = β ) 1 dτ r 2 D τ X i ) 2 + m 2 AX i ) 2, 5.2) with the mass m A = ±.7. For the flavoured model the BFSS transition is still present and when the X i become massive they induce a mass for the fundamental scalars and the induced bare mass for these is estimated by integrating out the adjoint fields and expanding it to quadratic order in Φ ρ. his gives a mass m f = 5 2m A However, the fundamental scalars are still strongly interacting as they have a selfcoupling of order one and we expect the bare mass to become significantly dressed. We therefore estimate the physical mass of the scalars at sufficiently low temperature by assuming that they also can be described by a massive gaussian with mass m f, in which case 1 r 2 bos = dτ tr β Φ ρ Φ ρ ) f m f ote that the right-hand side of equation 5.3) is independent of β and from Figure 1 we see that r 2 bos is more or less constant below the transition. A direct measurement of the expectation value 5.3) at =.5 gives m f, which gives the estimate m f However, at zero temperature we can extract the masses for the different fields by measuring their Green s function. o this end we set the holonomy to zero, the parameter β is now just the length of the time circle and not an inverse temperature. Because the SO9) symmetry of the bosonic BFSS model is broken down to SO5) SO4) there are now two bos 18

20 tr Α Α t t r X 9 X 9 t t Figure 2: Plots of the Green s functions equations 5.4) and 5.5) for β = 1, Λ = 144, = 1 and f = 1. he fits correspond to m f = and m l A = 2.1. adjoint masses, a longitudinal mass, m l A, for the four X ρ ρ and a transverse mass, m t A, for the five matrices X a. In figure 2 we present results for the Green s functions: tr Φ ρ )Φ ρ τ) = f m f e mf τ + e mf β τ) 1 e β m f rx 9 )X 9 τ) = 2m l A e ml A τ + e ml A β τ) 1 e β ml A, 5.4), 5.5) where we have chosen the last of the four SO4) adjoint scalars. We have also measured the longitudinal mass m l A by measuring the correlator for X1 defined similarly to 5.5). he results for β = 1, Λ = 144, = 1 and f = 1 are m l A = 2.1 ±.3, mt A = ±.3 and m f = 1.463±.1. he prediction from assuming that the adjoint fields are described by an action of the form 5.2) with different masses for the transverse and longitudinal matrices is: R 2 bos 5 2m t + 4 A 2m l = 2.27 ±.1, 5.6) A which agrees well with the direct measurement where we find R 2 bos = ±.5. Also the measured value of m f using 5.3) predicts that r 2 bos =.6836 ±.6, which is in excellent agreement with the measured value. ote that this estimate of the mass m f is very close to the one obtained from equation 5.3). Also the slightly different values of the adjoint masses m t A and ml A from the purely BFSS case considered in equation 5.2) reflect the presence of backreaction at f / =.1. Observe also the closeness of the transverse mass to the bosonic BFSS mass, which indicates that the backreaction is strongest for the longitudinal modes as one might expect. We can now use this information to estimate the value of C m bos at zero temperature. Assuming that both X a and Φ ρ are well approximated by massive gaussians and using Wick s 19

21 theorem on C m bos = 2 β β to perform the contractions, we obtain dτ tr Φ ρ Φ ρ 4 5β β dτ tr { Φρ X a } ) 2 Φ ρ 5.7) C m bos = = 2 f m f 2 f m 2 f mt A 2m f + m t = 1.27 ± ) A ) Finally, a direct measurement of the measured condensate shown in Figure 1 for =.4,.3,.2 and =.1 extrapolated to = gives C m bos = = ±.3, which is very close to the predicted value and confirms the validity of the gaussian approximation for both the adjoint and fundamental scalars. 6 Conclusions We have obtained the first two terms in the high-temperature series expansion for the Berkooz- Douglas model BD model) for general adjoint matrix size, and fundamental multiplet dimension, f. hese results should prove useful for future studies of this model. he model is an ideal testing ground for many ideas of gauge/gravity duality. he system is strongly coupled at low temperature while at high temperature it is weakly coupled, aside from the Matsubara zero-modes, which remain strongly coupled even at high temperature. It is these modes that provide the residual non-perturbative aspect of the current study. heir effect can be captured in numerical coefficients that depend only on and f. Once the coefficients are determined and tabulated see appendix A) they can be used as input for the high-temperature expansion of the observables of the BD model. We have checked these coefficients by comparing with a high precision simulation of the bosonic version of the BD model which we simulated using the Hybrid Monte Carlo approach. he coefficients in the high-temperature expansion of the bosonic model s observables are similarly determined by the tables presented in appendix A. In fact the observable r 2 bos see equation 2.8)) and mass susceptibility 5.7) of the model, shown in Figure 1, show that the agreement is excellent even down to temperature one. Below this temperature the system undergoes a set of phase transitions. hese are essentially the two phase transitions of the bosonic BFSS model. We found that for f / =.1 our measurements were sensitive to the backreaction of the fundamental fields on the adjoint fields. his backreaction lifted the mass degeneracy of the transverse and longitudinal adjoint fields. he transverse mass was essentially unaffected by the backreaction being m t A = ±.3 while the longitudinal mass was lifted to m l A = 2.1 ±.3, We found that using our understanding of the low-temperature phase of the BFSS model as a system of massive gaussian quantum matrix models we could predict the zero-temperature value of the mass susceptibility 5.7). he additional input that was required was the mass of the fundamental fields which we found by direct measurement to be m f = ±.1. 2

22 he zero-mode model used to obtain the high-temperature coefficients is of independent interest as it is the potential that captures the ADHM data in the theory of Yang-Mills instantons on the four-sphere, S 4. It is the bosonic sector of the dimensional reduction of the BD model to zero dimensions and is equivalent to a flavoured version of the bosonic sector of the IKK model. For this reason we refer to the model as the flavoured bosonic IKK model. he potential is always positive semi-definite and the Higgs branch of its zero-locus is isomorphic to the instanton moduli space [16]. here was some evidence for peculiar behaviour in the zero mode model for f 2. We found that simulations required significant fine tuning for f 2, in that when using the same leapfrog step length which gave 95% acceptance rate for f = 2 1 the acceptance rate for f 2 fell to a fraction of this within a couple of thousand sweeps and Ward identities we use as checks on the simulations were not fulfilled. After tuning the simulation we found the generated configurations had very long auto-correlation time. Also, in fitting the dependence of the observables Ξ i on for a given f we found evidence for a simple pole at = 2 f. Furthermore, one can see from the results tabulated in appendix A that they grow rapidly when the region f = 2 is approached. We expect that these difficulties and the growth of observables as f = 2 is approached are related to the singular structure of the instanton moduli space, i.e. the minimum of the potential in 2.3) with X a =, D A =. We have not pursued this further in the current study as it would take us too far afield, however, we believe it merits further attention. Finally, our preliminary studies of the supersymmetric BD model show [23] that, for some observables, the high-temperature series expansion remains valid to lower temperatures than one might expect. his validity of the high expansion at lower could provide alternative quasi-analytic estimates for observables in the window where gauge/gravity duality is valid. Acknowledgment D.O C. thanks Stefano Kovacs and Charles ash for helpful discussions on the ADHM construction. Support from Action MP145 QSPACE of the COS foundation is gratefully acknowledged. 21

23 A ables for the ω s. In this appendix we gather the numerical data from Monte Carlo simulations for different matrix sizes, and different numbers of flavour multiplets f and present it in tabular form. f = 1 ω 1 ω 3 ω 4 ω 6 ω 1,1 ω 1,2 ω 1,3 ω 1, ) ).89742).9721).131) -.554) -.956) -.292) ).26274).87492).11158).141) -.688) -.978) -.293) ) ).86441) ).1461) -.756) -.976) -.293) ).26872).8661) ).1481) -.787) -.977) -.33) ).26484) ).12889).14884) -.792) -.972) -.291) ) ).85391).12393).1512) -.829) -.979) -.295) ) ).85121) ).1523) -.81) -.11) -.296) ) ).851) ).1542) -.91) -.11) -.297) ) ) ) ).1543) -.91) -.12) -.38) ).25832).84642) ).15539) -.874) -.975) -.293) ) ) ) ).1582) -.919) -.11) -.297).25533) ).83461) ).16192) ) ) ) f = 1 ω 3,3 ω 3,4 ω 4,4 ω 5,5 ω 1,5,5 ω 2,5,5 ω 3,5,5 ω 4,5, ).331).352).7791) -.594).331).25).4836) ).312).2892).7778) -.648).341).218).4726) ).32).282).77269) -.61).332).21).4618) ).32).2772).7719) -.72).342).22).4589) ).31).27526).76993) -.656).3448).27).4553) ).294).2722).76658) -.63).313).23).4519) ).295).2692).76347) -.73).354).23).4419) ).297).2673).76298) -.75).356).25).441) ).298).2663).7627) -.76).387).26).441) ).283).26488).75992) -.62).342).22).4383) ).287).2611).75592) -.75).355).25).436).33526).2721).25421).74796) -.715).362).1886).4162) able 1: Mean values of the ω s for f = 1 were obtained from Monte Carlo samples for all values of but samples) and 2, samples). Errors are estimated with the Jackknife resampling. = values are the one extrapolated by quadratic polynomials: a + b/ + c/ 2. he quoted errors of this extrapolation are the fitting errors of the parameter a. 22

24 In the remaining tables we tabulate fixed = 9, 12, 14, 16, 18 and 2 while we vary f. Mean values of observable ω i, ω i,j and ω i,j,k i, j, k = 1,..., 6) were obtained from samples generated by hybrid Monte Carlo simulation of flavoured bosonic IKK model with the action specified in 3.5). Errors are estimated with the Jackknife resampling. = 9 f ω 1 ω 3 ω 4 ω 6 ω 1,1 ω 1,2 ω 1,3 ω 1, ) ) ).19363).1341) -.627) -.997) -.33) ) ) ).9912).1939) -.355) -.176) -.332) ).31312) ).74382).8878) -.154) ) -.362) ) ) ).59412).796).24) ) -.422) ) ) ).4571).5516).164) ) -.522) ).44727) ).33131).4115).273) -.286) -.72) ).55322) ).21391).2845).374) ) -.183) ).8226) ).1361).1655).454) ) -.234) = 9 f ω 3,3 ω 3,4 ω 4,4 ω 5,5 ω 1,5,5 ω 2,5,5 ω 3,5,5 ω 4,5, ).343).313).8897) -.61).342).22).531) 4.458).443).3975).89718) -.71).332).42).723) 6.541).614).541).1151) -.72).312).63).156) 8.692).916).782).11791) -.82).323).113).161) 1.953).151) 1.256).14232) -.12).33).195).282) ).272) 2.31).18162) -.142).33).438).554) ).666) 5.54).25273) -.212).273).112) 1.31) ) 2.93) 24.2).42759) -.434).15).527) 5.65) = 12 f ω 1 ω 3 ω 4 ω 6 ω 1,1 ω 1,2 ω 1,3 ω 1, ).26611) ) ).1412) -.699) -.11) -.294) ) ).92867).1112).1211) -.487) -.138) -.314) ) ) ).8772).141) -.37) -.199) -.333) ) ) 1.381).75362).891) -.157) ) -.373) ) ) ).63872).751) -.26) -.139) -.413) ) ) ).53161).6318).95) ) -.463) ).39954) ).43141).5187).195) -.178) -.563) ) ) ).33681).4147).275) -.288) -.74) 23

25 = 12 f ω 3,3 ω 3,4 ω 4,4 ω 5,5 ω 1,5,5 ω 2,5,5 ω 3,5,5 ω 4,5, ).324).2953).79447) -.63).333).33).52) 4.411).394).3516).85448) -.73).333).33).624) 6.471).485).431).92699) -.63).314).44).797) 8.542).617).542).11899) -.73).313).74).111) 1.653).829).73).11361) -.83).314).16).152) 12.83).111).975).12921) -.83).283).147).212) ).172) 1.428).1582) -.113).334).239).324) ).273) 2.31).18162) -.143).334).41).546) = 14 f ω 1 ω 3 ω 4 ω 6 ω 1,1 ω 1,2 ω 1,3 ω 1, ) ) ).1182).1432) -.71) -.11) -.295) ).27671).9787).15412).1272) -.548) -.11) -.314) ) ).95276).93632).1112) -.378) -.111) -.324) ).33551) ).82622).981) -.238) -.111) -.354) ).32442) ).72311).851) -.117) ) -.384) ).34312) ).62591).7379).6) ) -.413) ) ) ).53441).6329).16) ) -.473) ) ) ) ).5368).185) ) -.544) = 14 f ω 3,3 ω 3,4 ω 4,4 ω 5,5 ω 1,5,5 ω 2,5,5 ω 3,5,5 ω 4,5, ).315).2884).7878) -.74).354).34).492) 4.41).375).3337).83727) -.63).324).34).584) 6.442).445).391).89589) -.64).325).45).717) 8.52).546).472).96558) -.83).354).55).99) 1.572).678).582).1511) -.74).34).76).122) ).869).733).11561) -.74).35).18).152) 14.83).111).974).12971) -.94).335).158).213) ).162) 1.348).14731) -.14).35).21).35) 24

26 = 16 f ω 1 ω 3 ω 4 ω 6 ω 1,1 ω 1,2 ω 1,3 ω 1, ) ) ) ).1462) -.81) -.11) -.296) ) ).89685).1872).1312) -.61) -.11) -.35) ) ) ).98212).1172) -.41) -.11) -.314) ) ).97446).88261).142) -.39) -.111) -.335) ).39831) 1.266).7891).932) -.188) -.121) -.364) ) ) 1.747).721).821) -.88) -.121) -.385) ) ) ).61591).721).28) -.131) -.424) ) ) ).5361).631).17) -.151) -.474) = 16 f ω 3,3 ω 3,4 ω 4,4 ω 5,5 ω 1,5,5 ω 2,5,5 ω 3,5,5 ω 4,5, ).315).2834).78237) -.75).355).25).482) 4.391).355).3216).82398) -.64).326).35).565) 6.432).416).371).87368) -.75).375).36).668) 8.472).486).431).92958) -.74).335).56).81) 1.522).588).512).9961) -.85).366).68).12) ).71).613).1741) -.76).336).11).122) ).91).764).11751) -.86).317).11).163) ).112).976).12931) -.95).346).11).214) = 18 f ω 1 ω 3 ω 4 ω 6 ω 1,1 ω 1,2 ω 1,3 ω 1, ) ).86767) ).1473) -.81) -.12) -.297) ) ) ) ).1343) -.61) -.11) -.36) ).2869).92725).1182).1212) -.51) -.11) -.325) ).29871) ).92791).112) -.41) -.111) -.335) ).32261).99416).84181).992) -.21) -.111) -.345) ) ) ).7611).892) -.159) -.121) -.365) ) ) ).68251).82) -.69) -.131) -.395) ) ) ).68149).711).28) -.131) -.425) 25