How Prof. Zeidler Supported Our Research on Exact Solution of Quantum Field Theory Toy Models

Over many years, we developed the construction of the φ 4 -model on four-dimensional Moyal space. The solution of the related matrix model Z [ E, J ] = (cid:2) d(cid:3) exp ( tr (J(cid:3) − E(cid:3) 2 − λ 4 (cid:3) 4 )) is given in terms of the solution of a non-linear equation for the 2-point function and the eigenvalues of E . The resulting Schwinger functions in position space are symmetric and invariant under the full Euclidean group. Locality is fulfilled. The Schwinger 2-point function is reflection positive in special cases. The quantum field theory is constructed as a perturbative expansion about the free theory, which is solved by the propagator G mn ; kl , the inverse of (cid:14) mn ; kl . After diagonalisation

RW in December 2003. This changed life of all three of us. Vincent jumped to the subject and pushed it enormously, together with his group. They extended the vanishing of the beta-function to all orders in pertubation theory and led us to develop the non-perturbative solution after 2009. On the first version of these equations, we had an intense exchange with Prof. Zeidler in 2011. These discussions inspired our reduction of the equations to a fixed point problem in 2012. In a contribution to a special volume dedicated to Prof. Zeidler's 75th birthday, we were able to prove existence of a solution. We remember the long phone conversation with Prof. Zeidler about this fixed point problem.
In this contribution, we summarise the main developments in our programme since 2002, referring to [8][9][10][11][12][13][14][15][16]. Originally intended as renormalisation of the noncommutative φ 4 4 -model, our work turned into a programme to exactly solve quantum field theory toy models: 1. In December 1999, Minwalla, van Raamsdonk, and Seiberg pointed out [19] that Feynman graphs for scalar fields on noncommutative R 4 show a novel type of singularity (termed UV/IR-mixing) which prevented renormalisation. This came as big surprise, and for several years, no solution was found. Eventually, combining the Wilson-Polchinski programme for noncommutative φ 4 -theory with the harmonic oscillator base of the Moyal plane (which avoids the phase factors appearing in momentum space), we solved in a series of papers [9][10][11] the renormalisation problem. Thereby, we achieved the remarkable balance of proving renormalisability of the φ 4 -model to all orders and reconfirming the UV/IR-duality of [19]. The main steps are given in Section 2. 2. Soon after the renormalisation proof, we showed that the running coupling constant has bounded flow to one-loop order (see [8], reviewed in Section 3). This result led to a close collaboration with Vincent Rivasseau and his group. He emphasised that it should be possible to construct this model non-perturbatively! They established the foundation by proving that at a special self-duality point, the β-function vanishes to all orders [5]. We understood that their method goes beyond the β-function and used it to derive a closed non-perturbative equation for the 2-point function [12] (which we intensely discussed with Prof. Zeidler). 3. In [14], we gave a rigorous derivation of these equations. As reviewed in Section 4, Ward identities for the U(∞) group action lead to an exact solution of the quartic matrix model in terms of the solution of a non-linear equation. As by-product, we find that any renormalisable quartic matrix model has vanishing β-function. 4. Self-dual φ 4 4 -theory on Moyal space [10,11] is of that type. We give a summary of the proof in Section 4.3. The non-perturbative solution leads, for extreme noncommutativity θ → ∞, and after careful discussion of thermodynamic and continuum limit, to a non-linear fixed point equation [14], for which a non-perturbative and non-trivial solution exists for λ < 0 [16]. The key step is the observation that a certain difference function satisfies a linear singular integral equation of Carleman type. 5. Following [13], we identify in Section 5 a limit to Schwinger functions for a scalar field on R 4 . Surprisingly for a highly noncommutative model, these Schwinger functions show full Euclidean symmetry. Otherwise, they have unusual properties such as absent momentum transfer in interaction processes. This seems to suggest triviality, but the numerical investigation [15] of the 2-point function shows scattering remnants from a noncommutative geometrical substructure. Most surprisingly, the Schwinger 2-point function seems to be reflection positive in one of its phases.

Renormalisation of Noncommutative φ 4 -theory to All Orders
The renormalised φ 4 -model corresponds to the classical action withx μ := 2( −1 ) μν x ν and the star product is defined by [6,7]: The appearance of the harmonic oscillator term (1) is a result of the renormalisation proof, as sketched below. It also permits a transformation under Langmann-Szabo duality [17]. Relativistic quantum field theories on noncommutative Minkowski space are much more difficult [1]. Here the UV/IR-mixing problem occurs in different types of graphs [2].

The φ 4 -action in the Matrix Base
We assume for simplicity that θ 12 = −θ 21 = θ 34 = −θ 43 are the only non-vanishing components. A key step is expansion of the fields in the harmonic oscillator base [7] which in two dimensions reads with z ∈ C ≡ R 2 . Collecting them to f mn (x 1 , . . . , where m = m 1 m 2 , these functions satisfy (f mn f kl )(x) = δ nk f ml (x) and R 4 dxf mn (x) = (2πθ) 2 δ mn . Working out the action of Laplace operator and multiplication by x on f mn and expanding φ(x) = m,n mn f mn (x), the action (1) takes the form The quantum field theory is constructed as a perturbative expansion about the free theory, which is solved by the propagator G mn;kl , the inverse of mn;kl . After diagonalisation It is important that the sums in (5) are finite.

Renormalisation Group Approach to Dynamical Matrix Models
The (Euclidean) quantum field theory is defined by the partition function The idea inspired by Polchinski's renormalisation proof of commutative φ 4 -theory is to change the weights of the matrix indices in the kinetic part of S[ ] as a smooth function of an energy scale and to compensate this by a careful adaptation of the effective action L[ , ] such that Z[J ] becomes independent of the scale . If the modification of the weights of a matrix index m ∈ N is described by a function K m θ 2 , then the required -dependence of the effective action is given by the matrix Polchinski equation where In this section, we look for a perturbative solution of the matrix Polchinski equation (6). In terms of the expansion coefficients of the effective action, the matrix Polchinski equation (6) is represented by ribbon graphs (we suppress n → n in ribbon graphs): An internal double-line symbolises the propagator Q mn;kl ( ). In this way, very complicated ribbon graphs can be produced which cannot be drawn any more in a plane. A ribbon graph represents a simplicial complex for a Riemann surface and thus defines the topology of the Riemann surface on which it can be drawn. The Riemann surface is characterised by its genus g computable via the Euler characteristic of the graph, g = We have proven in [9] a power-counting estimation for these coefficients which relates the -scaling of a ribbon graph to the topology of the graph and to two asymptotic scaling dimensions of the differentiated cut-off propagator Q mn;kl ( ). As a result, if these scaling dimensions coincide with the classical momentum space dimensions, then all nonplanar graphs are suppressed by the renormalisation flow. This is a necessary requirement for the renormalisability of a model. On the other hand, as the expansion coefficients A (V ) m 1 n 1 ;...;m N n N [ ] carry an infinite number of matrix indices, the general power-counting estimation proven in [9] leaves, a priory, an infinite number of divergent planar graphs. These planar graphs require a separate analysis.

Power-Counting Behaviour of the Noncommutative φ 4 -model
The key is the integration procedure of the Polchinski equation (9), which involves the entire magic of renormalisation. We consider the example of the planar one-particle irreducible four-point function with two vertices, A (2,1,0)1PI m 1 n 1 ;...;m 4 n 4 . The Polchinski equation (9) provides the -derivative of that function: Performing the -integration of (10) from some initial scale 0 (sent to ∞ at the end) down to , we obtain A ( The second graph in the first line on the right-hand side and the graph in brackets in the last line are identical, because only the indices on the propagators determine the value of the graph. Moreover, the vertex in the last line in front of the bracket equals 1. Thus, differentiating (11) with respect to we obtain indeed (10). As a further check one can consider (11) for m = n = k = l = 0 0 . Finally, the independence of A (2,1,0)1PI mn;nk;kl;lm [ 0 ] on the indices m, n, k, l is built-in. This property is, for 0 → ∞, dynamically generated by the model.
There is a similar 0 -R -mixed integration procedure for the planar 1PI two-point functions A (V ,1,0)1PI and all other A (V ,1,0)1PI mn;nk;kl;lm . These involve in total four different sub-integrations from R up to . We refer to [10] for details. All other graphs are integrated from 0 down to , e.g.,

Theorem 1 The previous integration procedure yields
where P q [X] stands for a polynomial of degree q in X and | m 1 Idea of the proof. For the choice K(x) = 1 for 0 ≤ x ≤ 1 and K(x) = 0 for x ≥ 2 of the cut-off function in (7) one has Thus, the propagator and the volume of a loop summation have the same power-counting dimensions as a commutative φ 4 -model in momentum space, giving the total powercounting degree 4 − N for an N -point function. This is (more or less) correct for planar graphs. The scaling behavior of non-planar graphs is considerably improved by the quasi-locality of the propagator: As a consequence, for given index m of the propagator Q mn;kl ( ) = , the contribution to a graph is strongly suppressed unless the other index l on the trajectory through m is close to m. Thus, the sum over l for given m converges and does not alter (apart from a factor −1 ) the power-counting behaviour of (14): In a non-planar graph like the one in (12), the index n 3 -fixed as an external indexlocalises the summation index p ≈ n 3 . Thus, we save one volume factor θ 2 4 compared with a true loop summation as in (2.3). In general, each hole in the Riemann surface saves one volume factor, and each handle even saves two. A more careful analysis of (5) shows that also planar graphs get suppressed with as the only relevant or marginal ones. In these functions one has to use a discrete version of the Taylor expansion such as which can be traced back to the Meixner polynomials. The discrete Taylor subtractions are used in the integration from 0 down to in prescriptions like (11): This explains the polynomial in fractions like |m| θ 2 in (13).
As the estimation (13) is achieved by a finite number of initial conditions at R (see (11)), the noncommutative φ 4 -model with oscillator term is renormalisable to all orders in perturbation theory. These initial conditions correspond to normalisation experiments for the mass, the field amplitude, the coupling constant, and the oscillator frequency in the bare action related to (1).

Vanishing of the β-Function
Knowing the relevant/marginal couplings, we can compute Feynman graphs with sharp matrix cut-off N . The most important question concerns the β-function appearing in the renormalisation group equation which describes the cut-off dependence of the expansion coefficients m 1 n 1 ;...;m N n N of the effective action when imposing normalisation conditions for the relevant and marginal couplings. We have [8] lim N →∞ Here, Z is the wavefunction renormalisation. To one-loop order, we find [8] From (18) and (19), one finds that λ 2 remains constant under the renormalisation flow. The integration of the resulting differential equation shows that, starting from given small values for R , λ R at N R , the frequency grows in a small region around ln N N R = 48π 2 λ R to ≈ 1.

The coupling constant approaches
, which can be made small for sufficiently small λ R . This leaves the chance of a non-perturbative construction [23] of this model.
The key observation for all subsequent work is the following: The one-loop renormalisation flow has a non-trivial UV fixed point given by the self-dual model = 1 in (2) where β = 0 to all orders. Working exactly at = 1 leads according to (4) to the formulation as a matrix model where H mn = Z 4 θ (|m| + |n|) + (μ 2 bare + 2) and |m| := m 1 + m 2 . Up to an inessential factor 1 64π 2 and a shift μ 2 bare + 2 → μ 2 bare , this action can also be written as where This form, closely related to models studied in [18], will be the starting point of a general treatment of such models in Section 4. Starting from (20), Disertori and Rivasseau were able to prove that the β-function vanishes to three-loop order [4]. This result led to the conjecture of β = 0 to all orders, which was eventually proved in [5] by a combination of Ward identities and Schwinger-Dyson equation.
We understood immediately that the method suggested in [5] has the potential to provide an exact solution of the model. Indeed we proved in [12] that the 2-point function satisfies (in a certain limit) a closed non-linear integral equation, which is essentially equation (30) below, but expressed in other variables. In [12], we also gave a perturbative approximation of the solution. In the following years, we tried to extract non-perturbative information out of (27). We were joined by Prof. Zeidler in 2011 who tried to use techniques from non-linear functional analysis. Although a breakthrough was not achieved in this way, the exchange with Prof. Zeidler stimulated a different strategy via singular integral equations, which will be described in the next section.

Exact Solution of the Quartic Matrix Model
Adding a source term to the action (21), we define the partition function as where D[ ] is the extension of the Lebesgue measure from finite-rank operators to L 2 (H ) and J a test function matrix. For λ = 0 in (21), D[ ] Z[0] would be the Gaussian measure of covariance determined by E. In the presence of interaction λ = 0, a rigorous construction of the measure cannot be expected. Instead we will derive, for finite matrix size N , equations between connected correlation functions formally defined by Then, we prove that, after renormalisation, these equations have a well-defined limit N , V → ∞ which is exact in λ. We are then able to reduce this problem to a fixed point problem where analytic and numerical techniques are applied.

Ward Identity
The first steps apply for actions of the form (21) We can always place ourselves in an orthonormal basis of H where E is diagonal (but J is not). Since E is of compact resolvent, E has eigenvalues E a > 0 of finite multiplicity μ a . We thus label the matrices by an enumeration of the (necessarily discrete) eigenvalues of E and an enumeration of the basis vectors of the finite-dimensional eigenspaces. Writing in {. . . } of (23) as functional derivative ab = ∂ V ∂J ba , we have proved (first obtained in [5]):

Proposition 1 The partition function Z[J ] of the matrix model defined by the external matrix E satisfies the |I | × |I | Ward identities
Without loss of generality, we can assume that the map I m → E m ∈ R + is injective. Namely, correlation functions will only depend on the set of eigenvalues (E m ) of E.

Schwinger-Dyson Equations
We can write the action as where E a are the eigenvalues of E. Functional integration yields, up to an irrelevant constant, Instead of a perturbative expansion of e −V S int [ ∂ V ∂J ] , we apply those J -derivatives to (25) which give rise to a correlation function G ... on the left-hand side. On the right-hand side of (25), these external derivatives combine with internal derivatives from S int [ ∂ V ∂J ] to certain identities for G ... . These Schwinger-Dyson equations are often of little use because they express an N -point function in terms of (N + 2)-point functions.
In the field-theoretical matrix models under consideration, the Ward identity lets this tower of Schwinger-Dyson equation collapse. We prove in [14]:

) satisfies for injective m → E m the Schwinger-Dyson equation
It can be checked [14] that in a genus expansion G ... = ∞ g=0 V −2g G (g) ... (which is probably not convergent but Borel summable), precisely the line (26a) preserves the genus, the lines (26b) increase g → g + 1 and the line (26c) increases g → g + 2. In particular, in a scaling limit V → ∞ with 1 V p∈I finite, the exact Schwinger-Dyson equation for G |ab| coincides with its restriction (26a) to the planar sector g = 0, a closed non-linear equation for G (0) |ab| alone: We have derived in 2007/2008 this self-consistency equation for the Moyal model by the graphical method proposed by [5]. In this form, (27) is meaningless because p∈I diverges. In 2009 we solved the renormalisation problem, namely the renormalisation of infinitely many Feynman graphs at once [12]. This renormalisation increases the non-linearity. In [12] we have solved (27)  |a0| ! It turns out that a real theory with = * admits a short-cut which directly gives the higher N -point functions without any index summation. Since the equations for G ... are real and J ab = J ba , the reality Z = Z implies (in addition to invariance under cyclic permutations) invariance under orientation reversal Whereas empty for the identities (28) lead to many cancellations which result in a universal algebraic recursion formula: The last line of (29) increases the genus and is absent in G The theorem tells us that vanishing of the β-function for the self-dual 4 4 -model on Moyal space (proved in [5] to all orders in perturbation theory) is generic to all quartic matrix models, and the result even holds non-perturbatively! The universal recursion formula (29) computes the planar N -point function G |b 0 ...b N −1 | at B = 1 as a sum of fractions with products of 2-point functions in the numerator and products of differences of eigenvalues of E in the denominator. This structure admits an interesting graphical interpretation. We draw the indices b 0 , . . . b N−1 in cyclic order on the circle S 1 and represent a factor G b i b j as a chord connecting b i with b j and a factor 1 The chords form the non-crossing chord diagrams counted by the Catalan number

Renormalisation and Integral Representation
We return to the noncommutative φ 4 -model at self-duality point = 1, where the action is given by Eq. (20). Our general results on quartic matrix models imply that the planar 2-point function G (0) |ab| satisfies the self-consistency equation (27), where we recall We have introduced a cut-off N 2 N in the matrix size; the index sum diverges for N 2 N → N 2 . As usual, the renormalisation strategy consists in adjusting Z, μ bare in such a way that the limit N 2 N → N 2 exists. This will be achieved by normalisation conditions for the 1PI function ab defined by G (|p|+1)f (|p|). We study a particular scaling limit in which matrix size N and volume V are simultaneously sent to ∞ such that the ratio N √ V μ 4 = 2 (1+Y) is kept fixed. Note that V = θ 4 2 → ∞ is a limit of extreme noncommutativity! The new parameter (1+Y) corresponds to a finite wavefunction renormalisation, identified later to decouple our equations. The parameter 2 represents an ultraviolet cut-off which is sent to → ∞ in the very end (continuum limit). In the scaling limit, functions of In the same way, ren a b converges to a function μ 2 ab with a, b ∈ [0, 2 ], and the discrete sum converges to a Riemann integral This limit makes the restriction to the planar sector (27) of (26a)-(26c) exact. Applying d db | a=b=0 we get Z in terms of ab (and its derivative). Inserted back one gets a highly non-linear integro-differential equation. We can reduce the non-linearity by subtracting from it the same equation taken at b = 0. In terms of G ab := ((a + b)(1 + Y) + 1 − ab ) −1 , this difference equation reads Differentiation d db | a=b=0 of (31) yields Z in terms of G ab and its derivative. The resulting derivative G can be avoided by adjusting This choice leads to Z −1 (1+Y) = 1 − λ 2 0 dpG p0 , which is a perturbatively divergent integral for → ∞. Inserting Z −1 and Y back into (31) we end up in a linear integral equation for the difference function D ab := a b (G ab − G a0 ) to the boundary: The non-linearity restricts to the boundary function G a0 where the second index is put to zero. Assuming a → G ab Hölder-continuous, we can pass to Cauchy principal values. In terms of the finite Hilbert transform the integral equation (32) becomes

The Carleman Solution
Equation ( The possibility of C, C = 0 is due to the fact that the finite Hilbert transform has a kernel, in contrast to the infinite Hilbert transform with integration over R. The two formulae (35a) and (35b) are formally equivalent, but the solutions belong to different function classes and normalisation conditions may (and will) make a choice.
From (34), one introduces an angle function, which leads to a representation first for G a0 : where C is an arbitrary constant.
Recall that G a0 forms the inhomogeneity in the Carleman equation ( for λ < 0, where C is an undetermined constant and bF (b) an undetermined function of b vanishing at b = 0.
Some remarks: -We have proved this theorem in [14] for λ > 0 under the assumption C = 0 in (35b), but knew that non-trivial solutions of the homogeneous Carleman equation parametrised by C = 0 are possible. That no such term arises for λ < 0 (if angles are redefined ϑ → τ ) was proved in [15].
-We expect C, F to be -dependent so that 1 -An important observation is G ab ≥ 0, at least for λ < 0. This is a truly non-perturbative result; individual Feynman graphs show no positivity at all! -As in [12], the equation for G ab can be solved perturbatively. Matching at λ = 0 requires C, F to be flat functions of λ (all derivatives vanish at zero). Because of −→ −∞, the naïve arctan series is dangerous for λ > 0. Unless there are cancellations, we expect zero radius of convergence! -From (37), we deduce the finite wavefunction renormalisation -The partition function Z is undefined for λ < 0. But the Schwinger-Dyson equations for G ab and for higher functions, and with them log Z, extend to λ < 0. These extensions are unique but probably not analytic in a neighbourhood of λ = 0.
Given the boundary function G a0 , the Carleman theory computes the full 2-point function G ab via (37). In particular, we get G 0b as function of G a0 . But the 2-point function is symmetric, G ab = G ba , and the special case b = 0 leads to the following self-consistency equation:

Proposition 4
The limit θ → ∞ of φ 4 4 -theory on Moyal space is determined by the solution of the fixed point equation G = T G, At this point, we can eventually send → ∞.
In [14], we proved via the Schauder fixed point theorem that (39) has a (smooth) solution for λ > 0 (assuming F (b) = 0) bounded by 0 ≤ G b0 ≤ 1 1+b . For the much more complicated case λ < 0, we proved in our contribution [16] to the 75th birthday of Prof. Zeidler the following result: In [15], we solved (39) numerically by approximating G 0b as a piecewise linear function on [0, 2 ] sampled according to a geometric progression and by viewing (39) as iteration G i+1 0b = (T G i ) 0b for some initial function G 0 . We confirmed the convergence of this iteration in Lipschitz norm for a large range λ ∈ R of either sign. It turned out that the required symmetry G ab = G ba does not hold for λ > 0, which is a clear hint that F (b) = 0 for λ > 0. For λ < 0 everything is consistent within small numerical errors. From the solution of (39) we get G ab via (37) and then all higher correlation functions via the universal algebraic recursion formulae. For λ < 0, all these quantities of the model can be obtained with sufficient precision.
We find clear evidence in [15] for a second-order phase transition at λ c ≈ −0.39, which is a common critical value in several independent problems. The most obvious signal is a plot of the derivative 1 + Y := − dG 0b db | b=0 as function of λ shown in Fig. 1. Globally, we found that the numeric solution is close to (but not exactly) a power law G b0 ∝ 1 (1+b) 1−η(λ)/2 , where η has opposite sign as λ. This numerical conjecture was later made precise and proved in [16]. We have indications that the exact critical coupling constant will be λ c = − 1 π . We discuss in the next section how the sign on η relates to reflection positivity of Schwinger functions made from G ... . Reflection positivity requires η ≥ 0 which excludes (unless F (b) = 0 reverses the behaviour) the stable case λ > 0 and prefers λ ≤ 0 where the partition function is meaningless.

Reverting the Matrix Basis
We are interested here in the limit to Schwinger functions in position space. For this end, we revert the matrix representation and take the infinite volume limit V μ 4 → ∞, where we carefully have to pass to densities. Absolute position x ∈ R 4 have no meaning, only μx can be used. This means that we consider where the matrix correlation functions are formally given by (22)   , we symbolise the renormalisation of Section 4.3.
The main question is whether these Schwinger functions satisfy the Osterwalder-Schrader axioms [20,21] which would allow to analytically continue the model to a true Wightman quantum field theory. The first step consists in an explicit evaluation of (40). In [13], we proved the following result: Some comments: -Only a restricted sector of the underlying matrix model contributes to position space: All strands of the same boundary component carry the same matrix index. -Schwinger functions are symmetric and invariant under the full Euclidean group. This comes truly surprising since θ = 0 breaks both translation invariance and manifest rotation invariance. The limit θ → ∞ was expected to make this symmetry violation even worse! -The most interesting sector is the case where every boundary component has N β = 2 indices. It is described by the (2 + · · · + 2)-point functions G p 1 2 .
-This sector describes the propagation and interaction of B particles without any momentum exchange. This is acceptable for a 2D model. In four dimensions, the absence of momentum transfer is a sign of triviality. -However, typical triviality proofs rely on clustering, analyticity in Mandelstam representation or absence of bound states. All this needs verification.

Reflection Positivity
Reflection positivity is the most decisive Osterwalder-Schrader axiom. It gives the spectrum condition and positivity of the reconstructed Hilbert space of the Minkowski model [20,21]. This guarantees a representation as a Laplace transform in the time component of the coordinate ξ 0 , hence analyticity in Re(ξ 0 ) > 0.
We have shown in [13] that the Schwinger 2-point function S(x 1 , x 2 ) given by (41) is reflection positive if and only if a → G aa is a Stieltjes function ( [25]), for a positive and non-decreasing measure ρ. The proof follows from the Källén-Lehmann representation of the two-point function.
The numerical results [15] exclude reflection positivity for any λ > 0 (due to renormalisation). Interestingly, it thus favours the wrong-sign λφ 4 -model studied in [22]. A rigorous proof that G aa satisfies the Stieltjes property for λ < 0 is still missing, although the numerical results of [15] provide strong evidence that this will be true.

Summary
By applying the Wilson-Polchinski ideas to the noncommutative 4 -theory with harmonic oscillator term we were able to solve the renormalisation problem. We obtained renormalisability to all orders of pertubation theory and reconfirmed the UV/IR-duality. Next, we showed that the running coupling constant has bounded flow to one-loop order. This led to a close collaboration with Vincent Rivasseau and his group. They proved the essential result, that at a special self-duality point, the β-function vanishes to all orders in perturbation theory.
Ward identities are the reason behind this result. They allow to decouple the hierarchy of Schwinger-Dyson equations, which allows to solve the model.
We have shown that the φ 4 4 -model on noncommutative Moyal space, considered in the limit θ → ∞ of extreme noncommutativity, is an exactly solvable and non-trivial matrix model. Euclidean symmetry is violated in the beginning, but we identified a limit which projects to diagonal matrices where Euclidean symmetry is restored. Surprisingly, the first consistency checks for OS positivity are passed for the only interesting interval [λ c , 0] of the coupling constant! This model is somewhat strange as 'particles' keep their momenta in interaction processes. Nevertheless, the theory is not completely trivial. We find scattering remnants from the noncommutative geometrical (i.e., matricial) substructure. Only the external matrix indices are put 'on-shell', internally all degrees of freedom contribute. We have seen that clustering is maximally violated. The interaction is insensitive to positions in different boundary components. In particular, 'particles' are never asymptotically free.