The azimuthal angle (φ) variation of the Coulomb and nuclear heavy ion (HI) potentials is studied in the framework of the double folding model, which is derived from realistic nuclear density distributions and a nucleon-nucleon (NN) interaction. The present calculation shows that the variation of HI potentials with the azimuthal angle depends strongly on the range of the NN forces. For the long-range Coulomb force, the maximum variation with φ is about 0.9%, and for HI potential derived from zero-range NN interaction the φ-variation can reach up to 90.0%. Our calculations are compared with the recent φ-dependence of the HI potential derived from proximity method. The present realistic φ-dependence calculations of the HI potential is completely different from the results of the proximity calculations. DOI: 10.1103/PhysRevC.75.064610 PACS number(s): 24.10.−iThe calculations of the nucleus-nucleus potential between two deformed, oriented nuclei have been of much interest [1][2][3][4][5][6][7][8][9]. The double folding model [2,4,10] plays a fundamental role in deriving the heavy ion (HI) potential. For two deformed density distributions, the calculation of the nucleus-nucleus potential is a hard task due to the numerical computation of a six-dimensional integral, which is very time consuming. This problem is quite relevant because most nuclei have permanent and/or vibrational deformations. Recently, many authors have considered this problem to study the synthesis of new and superheavy elements [4,5,7,8], since the collision of deformed nuclei is a path to the far side of the proposed island of superheavy nuclei [11]. For the collision between either two spherical [12], or spherical-deformed nuclei [13], the HI potential (including the Coulomb part) had been derived microscopically by different methods. The double folding model [2,4,10] is one of the successful methods to derive the HI potential starting from a finite range NN force. Moreover, this model is the only one used to treat correctly the Coulomb potential between two heavy ions.For two deformed nuclei, the double folding model can be simplified when their symmetry axes are coplanar [1,2]. For arbitrary orientations of the nuclei, the six-dimensional integral of the double folding model is difficult to calculate without making approximations. Due to this difficulty, many authors considered different approximate methods to derive both the nuclear and Coulomb HI potentials [5][6][7]. For example, in Refs. [6,7] the pocket formula of the proximity potential was used to study orientation dependence of the HI potential. In another study [14], the zero-range NN force was used to reduce the six-dimensional integration to three dimensions. The Coulomb potential, in these recent studies, was calculated from a simplified equation derived by Wong [15], or by assuming the nucleus to be of uniform charge distribution with sharp cutoff edge [5]. The two methods for * ali@mailer.eun.eg calculating the Coulomb potential produce a large error in both internal and surface reg...