The vibrational coupled cluster method in bosonic representation is formulated to describe the molecular anharmonic vibrational spectra. The vibrational coupled cluster formalism is based on Watson Hamiltonian in normal coordinates. The vibrational excited states are described using coupled cluster linear response theory (CCLRT). The quality of the coupled cluster wave function is analyzed. Specifically, the mean displacement values of the normal coordinates and expectation values of the square of the normal coordinates of different vibrational states are calculated. A good agreement between the converged full CI results and coupled cluster results is found for the lower lying vibrational states.