We calculate the production of exotic electromagnetic bound systems, which ia an consisting of a ($l^{+}l^{-}$) bound state system such as positronium, dimuonium, and ditauonium, by photon-photon interaction with the equivalent photon approximation in ultra-peripheral heavy ion collisions considering the nuclear form factor. The numerical results demonstrate that the experimental study of positronium, dimuonium, and ditauonium in ultra-peripheral collisions is feasible at Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) energies.
Comment: 7 pages, 2 figures