Assessing the entanglement of three coupled harmonic oscillators

Quantum entanglement serves as a key phenomenon in understanding correlations in many-body systems, but analytical results remain scarce for coupled three-body oscillators. In this work, we address this gap by introducing a geometrical diagonalization approach that constrains Euler angles, thereby reducing the degrees of freedom in the entanglement analysis. It consists of deriving analytical expressions for linear entropy and purity under the bipartitions $(x|yz)$, $(y|xz)$, and $(xy|z)$ using the Wigner function framework. Our results indicate that excitations in any oscillator basically enhance the redistribution of correlations across the system. The mixing angle $θ$ governs entanglement intensity, ranging from separability to maximal correlation. Moreover, we reveal the symmetry relations, notably $S_{Ly}[(n,m,l),θ]=S_{Lz}[(n,m,l),-θ]$ and an intrinsic symmetry within $(x|yz)$. Hence, we clarify how excitation levels and mixing angles create and enhance entanglement in the three coupled harmonic oscillators.