A series of calculations on the self-diffusion behavior of high entropy CoCrFeMnNi alloy were carried out using molecular dynamics methods. By computing both vacancy formation energy and atomic migration energy of the constituent elements in the alloy, the diffusional activation energy of each element is obtained, and the self-diffusion coefficients for bulk diffusion were calculated, with the values exhibiting close to of experiments. A model for structures of symmetrically tilted grain boundary is established, with Σ9 and Σ27 grain boundaries studied based on the coincidence site lattice theory. Measured by the full width at half maxima of the radial distribution function, it is found that the grain boundaries with low index are more ordered than those with high plane index, and the atom fluctuation occurred in the low-indexed grain boundaries is less intensively and sensitively to temperature change. Meanwhile, the diffusion coefficients of ordered grain boundaries are generally smaller than those of disordered grain boundaries. Compared with the experimental values of grain boundary diffusion, the diffusion activation energy of configured grain boundaries from coincidence site lattice is smaller than that of normal large-angle grain boundaries.