Unitarily localizable entanglement of Gaussian states

We consider generic $(m\ifmmode\times\else\texttimes\fi{}n)$-mode bipartitions of continuous-variable systems, and study the associated bisymmetric multimode Gaussian states. They are defined as $(m+n)$-mode Gaussian states invariant under local mode permutations on the $m$-mode and $n$-mode subsystems. We prove that such states are equivalent, under local unitary transformations, to the tensor product of a two-mode state and of $m+n\ensuremath{-}2$ uncorrelated single-mode states. The entanglement between the $m$-mode and the $n$-mode blocks can then be completely concentrated on a single pair of modes by means of local unitary operations alone. This result allows us to prove that the PPT (positivity of the partial transpose) condition is necessary and sufficient for the separability of $(m+n)$-mode bisymmetric Gaussian states. We determine exactly their negativity and identify a subset of bisymmetric states whose multimode entanglement of formation can be computed analytically. We consider explicit examples of pure and mixed bisymmetric states and study their entanglement scaling with the number of modes.