The equation of state of neutron stars is studied in the newly developed density dependent relativistic Hartree-Fock (DDRHF) theory with the effective interaction PKO1 and applied to describe the properties of neutron stars. The results are compared with the recent observational data of compact stars and those calculated with the relativistic mean field (RMF) effective interactions. The maximum mass of neutron stars calculated with PKO1 is about 2.45 M☉, which consists with high pulsar mass from PSR B1516+02B recently reported. The influence of Fock terms on the cooling of neutron stars is discussed as well.