Biphonons in the Klein-Gordon lattice

A numerical approach is proposed for studying the quantum optical modes in the Klein-Gordon lattices where the energy contribution of the atomic displacements is non-quadratic. The features of the biphonon excitations are investigated in detail for different non-quadratic contributions to the Hamiltonian. The results are extended to multi-phonon bound states.