  Quantum computers are promising tools for simulating many-body quantum
systems due to their potential scaling advantage over classical computers.
While significant effort has been expended on many-fermion systems, here we
simulate a model entangled many-boson system with the contracted quantum
eigensolver (CQE). We generalize the CQE to many-boson systems by encoding the
bosonic wavefunction on qubits. The CQE provides a compact ansatz for the
bosonic wave function whose gradient is proportional to the residual of a
contracted Schr\"odinger equation. We apply the CQE to a bosonic system, where
$N$ quantum harmonic oscillators are coupled through a pairwise quadratic
repulsion. The model is relevant to the study of coupled vibrations in
molecular systems on quantum devices. Results demonstrate the potential
efficiency of the CQE in simulating bosonic processes such as molecular
vibrations with good accuracy and convergence even in the presence of noise.
