Heavy quark pair correlations in QCD

The azimutal correlations of heavy quarks produced in the high energy $pp$ ($p\bar{p}$) collisions are calculated in the framework of QCD without the usual assumptions of the parton model. The virtual nature of the interacting gluons as well as their transverse motion and different polarizations are taken into account. We give some predictions for the azimutal correlations of charm and beauty hadrons produced at Tevatron-collider and LHC.