  We calculate the photoproduction of heavy quarks in proton-proton collisions
at RHIC, Tevatron and LHC energies. The integrated cross section and the
rapidity distributions for open charm and bottom production are computed
employing sound high energy QCD formalisms. For the linear pQCD approaches we
consider both the usual collinear factorization and the k_T-factorization
formalisms, whereas for the nonlinear QCD (saturation) calculations one
considers the Golec-Biernat-Wuesthoff and the Iancu-Itakura-Munier
parameterizations for the dipole cross section within the color dipole picture.
