Entangling power of permutations

The notion of entangling power of unitary matrices was introduced by Zanardi et al. Phys. Rev. A 62 030301 (2000)]. We study the entangling power of permutations, given in terms of a combinatorial formula. We show that the permutation matrices with zero entangling power are, up to local unitaries, the identity and the swap. We construct the permutations with the minimum nonzero entangling power for every dimension. With the use of orthogonal latin squares, we construct the permutations with the maximum entangling power for every dimension. Moreover, we show that the value obtained is maximum over all unitaries of the same dimension, with a possible exception for 36. Our result enables us to construct generic examples of 4-qudit maximally entangled states for all dimensions except for 2 and 6. We numerically classify, according to their entangling power, the permutation matrices of dimensions 4 and 9, and we give some estimates for higher dimensions.