The use of crystalline ion beams can increase luminosity in the collider and in experiments with targets for investigation of rare radioactive isotopes. The ordered state of circulating ion beams was observed at several storage rings: NAP-M (Proceedings of the Fourth All Union Conference on Charged Particle Accelerators, Vol. 2, Nauka, Moscow, 1975 (in Russian); Part. Accel. 7 (1976) 197; At. Energy 40 (1976) 49; Preprint CERN/PS/AA 79-41, Geneva, 1979) (Novosibirsk), ESR (Phys. Rev. Lett. 77 (1996) 3803) and SIS (Proceedings of EPAC’2000, 2000) (Darmstadt), CRYRING (Proceedings of PAC’2001, 2001) (Stockholm) and PALLAS (Proceedings of the Conference on Applications of Accelerators in Research and Industry, AIP Conference Proceedings, p. 576, in preparation) (M . unchen). New criteria of the beam orderliness are derived and verified with a new program code. Molecular dynamics technique is inserted in BETACOOL program (Proceedings of Beam Cooling and Related Topics, Bad Honnef, Germany, 2001) and used for numerical simulation of crystalline beams. The sudden reduction of momentum spread in the ESR experiment is described with this code. The simulation shows a good agreement with the experimental results. The code has then been used to calculate characteristics of the ordered state of ion beams for the MUSES Ion Ring (IR) (MUSES Conceptual Design Report, RIKEN, Japan, 2001) in collider mode. A new strategy of the cooling process is proposed which permits to increase significantly the linear density of the ordered ion beam and thereby the luminosity of electron–ion colliding experiments. r 2004 Published by Elsevier B.V.