For entangled polymers, there exist considerable contributions of the cross-correlation between different chains in the relaxation modulus, though the contribution is neglected in the single-chain models. Earlier studies have suggested that the cross-correlation is due to the inter-chain interactions such as the force balance around entanglements and the osmotic force suppressing the density fluctuations. However, the origin of the cross-correlation has been yet to be clarified. In this study, a new multi-chain slip-link model has been developed for the simulations of polymer dynamics without the inter-chain forces while the creation of entanglement is performed according to the geometrical manner and the local equilibration. The simulations reproduced the qualitative features of entangled polymer dynamics for the molecular weight dependence of the longest relaxation time and the diffusion coefficient, just as expected from the success of the earlier single-chain models without the inter-chain forces. On the other hand, the remarkable feature of the model is that the cross-correlation exists between chains, suggesting that the cross-correlation is generated via the local equilibration at the creation of entanglement even without the global equilibration by the inter-chain forces.