Solving the QCD Hamiltonian for bound states

The systematic approach to study bound states in quantum chromodynamics is presented. The method utilizes nonperturbative flow equations in the confining background, that makes possible to perform perturbative renormalization and to bring the QCD Hamiltonian to a block-diagonal form with the number of quasiparticles conserving in each block. The effective block-diagonal Hamiltonian provides constituent description for hadron observables. The renormalized to the second order effective Hamiltonian of gluodynamics in the Coulomb gauge is obtained at low energies. The masses for scalar and pseudoscalar glueballs are predicted.