Weingarten functions are used for evaluating integrals over the unitary group Ud of products of matrix coefficients of the form
where ∗ {\displaystyle *} denotes complex conjugation. Note that U j i ∗ = ( U † ) i j {\displaystyle U_{ji}^{*}=(U^{\dagger })_{ij}} where U † {\displaystyle U^{\dagger }} is the conjugate transpose of U {\displaystyle U} , so one can interpret the above expression as being for the i 1 j 1 … i q j q j 1 ′ i 1 ′ … j q ′ i q ′ {\displaystyle i_{1}j_{1}\ldots i_{q}j_{q}j'_{1}i'_{1}\ldots j'_{q}i'_{q}} matrix element of U ⊗ ⋯ ⊗ U ⊗ U † ⊗ ⋯ ⊗ U † {\displaystyle U\otimes \cdots \otimes U\otimes U^{\dagger }\otimes \cdots \otimes U^{\dagger }} .
This integral is equal to
where Wg is the Weingarten function, given by
where the sum is over all partitions λ of q (Collins 2003). Here χλ is the character of Sq corresponding to the partition λ and s is the Schur polynomial of λ, so that sλd(1) is the dimension of the representation of Ud corresponding to λ.
The Weingarten functions are rational functions in d. They can have poles for small values of d, which cancel out in the formula above. There is an alternative inequivalent definition of Weingarten functions, where one only sums over partitions with at most d parts. This is no longer a rational function of d, but is finite for all positive integers d. The two sorts of Weingarten functions coincide for d larger than q, and either can be used in the formula for the integral.
The first few Weingarten functions Wg(σ, d) are
where permutations σ are denoted by their cycle shapes.
There exist computer algebra programs to produce these expressions.12
The explicit expressions for the integrals of first- and second-degree polynomials, obtained via the formula above, are: ∫ U d d U U i j U ¯ k ℓ = δ i k δ j ℓ Wg ( 1 , d ) = δ i k δ j ℓ d . {\displaystyle \int _{U_{d}}dUU_{ij}{\bar {U}}_{k\ell }=\delta _{ik}\delta _{j\ell }\operatorname {Wg} (1,d)={\frac {\delta _{ik}\delta _{j\ell }}{d}}.} ∫ U d d U U i j U k ℓ U ¯ m n U ¯ p q = ( δ i m δ j n δ k p δ ℓ q + δ i p δ j q δ k m δ ℓ n ) Wg ( 1 2 , d ) + ( δ i m δ j q δ k p δ ℓ n + δ i p δ j n δ k m δ ℓ q ) Wg ( 2 , d ) . {\displaystyle \int _{U_{d}}dUU_{ij}U_{k\ell }{\bar {U}}_{mn}{\bar {U}}_{pq}=(\delta _{im}\delta _{jn}\delta _{kp}\delta _{\ell q}+\delta _{ip}\delta _{jq}\delta _{km}\delta _{\ell n})\operatorname {Wg} (1^{2},d)+(\delta _{im}\delta _{jq}\delta _{kp}\delta _{\ell n}+\delta _{ip}\delta _{jn}\delta _{km}\delta _{\ell q})\operatorname {Wg} (2,d).}
For large d, the Weingarten function Wg has the asymptotic behavior
where the permutation σ is a product of cycles of lengths Ci, and cn = (2n)!/n!(n + 1)! is a Catalan number, and |σ| is the smallest number of transpositions that σ is a product of. There exists a diagrammatic method3 to systematically calculate the integrals over the unitary group as a power series in 1/d.
For orthogonal and symplectic groups the Weingarten functions were evaluated by Collins & Śniady (2006). Their theory is similar to the case of the unitary group. They are parameterized by partitions such that all parts have even size.
Z. Puchała and J.A. Miszczak, Symbolic integration with respect to the Haar measure on the unitary group in Mathematica., arXiv:1109.4244 (2011). https://arxiv.org/abs/1109.4244 ↩
M. Fukuda, R. König, and I. Nechita, RTNI - A symbolic integrator for Haar-random tensor networks., arXiv:1902.08539 (2019). https://arxiv.org/abs/1902.08539 ↩
P.W. Brouwer and C.W.J. Beenakker, Diagrammatic method of integration over the unitary group, with applications to quantum transport in mesoscopic systems, J. Math. Phys. 37, 4904 (1996), arXiv:cond-mat/9604059. https://arxiv.org/abs/cond-mat/9604059 ↩