I would use RAxML, Phyml or Paml. The former is the fastest one. These programs allow you to fix the topology (species tree) and optimise only the branch lengths.
Most phylogenetic programs estimate branch lengths for you. You just need to set it in the settings so that it will print/store the branch length information in the resulting tree. You can view the branch lengths in tree visualization programs such as figtree or extract the values using R
But as Frederico said above, it might be good to also explore your data with LG plus gamma, or other models. Do you have DNA data too, or just protein?