The answer is positive, modulo the inherent non-rigor in the path integral quantization. There is a series of papers by Antti Niemi and collaborators where this problem was considered. I'll give you here an introduction to their main reasoning and some references.
The Itzykson-Zuber formula belongs to a special case of integrals where the WKB semiclassical approximation is exact due to the Duistermaat-Heckmann theorem (which is also described in Terry Tao's lecture).
Physically, the integrals satisfying this property can be expressed as classical partition functions:
$$ Z = \int_{\mathcal{M}} e^{-\beta H(x)} d{\mu}(x)$$
Where the integration is over a classical phase space $\mathcal{M}$ which is a symplectic manifold, and $ d{\mu}(x)$ is the Liouville measure on the manifold, which can be locally written as:
$$d{\mu}(x) = det(\omega) d_L(x)$$
Where $ \omega$ is the symplectic form and $d_L$ is the Lebesgue measure.
In the Itzykson-Zuber formula, the integrand is invariant under the action of the maximal torus $T^n$ of $U(n)$. Thus the integration can be effectively performed on the flag manifold $U(n)/T^n$, which is a symplectic manifold.
In physics, the exactness of the WKB approximation is attributed to the existence of supersymmetry.
This property is explained by A Keski-Vakkuri, Niemi, Semenoff and Tirkkonen in the following work. I'll repeat the main arguments (which I'll bring here their main arguments heuristically)
This supersymmetry is obtained by exponentiating the determinant factor of the symplectic form by means of a Berezin Grassmann integral:
$$ Z = \int_{\mathcal{M}} e^{-\beta ( H(x) + c^a \omega_{ab} c^b)} d_L(x)$$
One can easily check that the integrand is invariant under the following supersymmetry operator:
$$ Q = c^a \partial_a + \omega^{ab}\partial_aH \frac{\partial}{\partial c^b}$$,
i.e.,
$$[Q, H(x) + c^a \omega_{ab} c^b] = 0$$
Thus:
$$[Q, e^{-\beta ( H(x) + c^a \omega_{ab} c^b)}] = 0$$
This operator squares to:
$$Q^2 = c^a \partial_a( \omega^{bc} \partial_c H )\frac{\partial}{\partial c^b}$$
Thus this operator vanishes in the stationary points $\partial_a H = 0$.
The Berezin integral picks the top form, the action of Q on the top form has two components, one which reduces the ghost number and the second increases the ghost number, which acts as an exterior derivative. Thus the top form is exact in all points where $Q$ is not zero, that is when $\partial_a H \neq 0$. Thus only points where the Hamiltonian is stationary contribute to the integral. Thus the semiclassical approximation is exact.
Niemi generalized the same reasoning for path integrals in 0+1 dimensions.
Final remarks
The above heuristic reasoning seems to apply to any Hamiltonian, but a deeper analysis shows that it applies only to Hamiltonians which are perfect Morse functions on the phase space.
Witten applied the same reasoning even in field theory computations for example $N=2$ Susy Yang-Mills theory in $4D$.
This post imported from StackExchange Physics at 2017-10-10 21:08 (UTC), posted by SE-user David Bar Moshe