In this paper, molecular dynamics simulations (MDS) of symmetric tilt ice grain boundaries are presented. The MDS were carried out using the GROMACS v4.5.5 program, and the water molecules were described using the TIP5P-Ew model. The grain boundary energies, , relative to those of the surface free energies, , were obtained as a function of the misorientation angle Ψ, and compared with the values experimentally obtained. The results show a good correspondence between the experimental and simulated values. The planar density of coincidence sites at the grain boundary planes, Γ, was obtained as a function of ψ. The Γ values were compared with the simulated values and a relation between the minimum of the simulated values and the maximum of the Γ values was observed, suggesting that the CSL theory is a good starting point to detect low energy ice GBs.