Large N Optimization for multi-matrix systems

In this work we revisit the problem of solving multi-matrix systems through numerical large $N$ methods. The framework is a collective, loop space representation which provides a constrained optimization problem, addressed through master-field minimization. This scheme applies both to multi-matrix integrals ($c=0$ systems) and multi-matrix quantum mechanics ($c=1$). The complete fluctuation spectrum is also computable in the above scheme, and is of immediate physical relevance in the later case. The complexity (and the growth of degrees of freedom) at large $N$ have stymied earlier attempts and in the present work we present significant improvements in this regard. The (constrained) minimization and spectrum calculations are easily achieved with close to $10^4$ variables, giving solution to Migdal-Makeenko, and collective field equations. Considering the large number of dynamical (loop) variables and the extreme nonlinearity of the problem, high precision is obtained when confronted with solvable cases. Through numerical results presented, we prove that our scheme solves, by numerical loop space methods, the general two matrix model problem.