Abstract
We study the entanglement entropy of harmonic oscillators in noncommutative phase space. We propose a new definition of quantum R\'enyi entropy based on Wigner functions in noncommutative phase space. Using the R\'enyi entropy, we calculate the entanglement entropy of the ground state of the 2D isotropic harmonic oscillators. We find that for some values of the noncommutative parameters, the harmonic oscillators can be entangled in noncommutative phase space. This is a new entanglement-like effect caused by the noncommutativity of the phase space.