The real part of a (nuclear) optical potential in heavy-ion collisions, which could have a phenomenological Woods-Saxon shape, can be determined by the Broglia & Winther parametrization. For instance, you can see the review paper by W. Reisdorf, Journal of Physics G 20 (1994) 1297, where this parametrization is described. It is connected with the double-folding nuclear potential that uses Fermi functions, parametrizing the densities of the interacting nuclei, and the M3Y effective nucleon-nucleon interaction. It is also possible to calculate a single-folding or a cluster-folding nuclear interactions. These methods are based on integrals involving the density of the interacting objects and an effective interaction among them
@Alexis. Thank you very much! I will read this paper now. I am now using Fresco code , written by Prof. I. J. Thompson, to calculating heavy ion 3-CDCC breakup. So I need W-S optical potential. Is it right: we could calculate real part by folding or parametrization, but we have to test imaginary part until it match the elastic scattering data?