Entanglement of bosonic modes in symmetric graphs

The ground and thermal states of a quadratic hamiltonian representing the interaction of bosonic modes or particles are always Gaussian states. We investigate the entanglement properties of these states for the case where the interactions are represented by harmonic forces acting along the edges of symmetric graphs, i.e. 1, 2, and 3 dimensional rectangular lattices, mean field clusters and platonic solids. We determine the Entanglement of Formation (EoF) as a function of the interaction strength, calculate the maximum EoF in each case and compare these values with the bounds found in \cite{wolf} which are valid for any quadratic hamiltonian.