Classical Monte Carlo code for the determination of atomic cross sections for fusion related collision systems

Complete text of publication follows. Interpretation of the cross sections in multielectron ion-atom collisions is a challenging task for theories. The main difficulty is caused by the many-body feature of the collision, involving the projectile, projectile electron(s), target nucleus, and target electron(s). In the last decade large numbers of non-perturbative studies have been performed to explain experimental total-, single- and double-ionization and charge exchange cross sections. The success of various approaches depends strongly on how far a given theory is capable to describe the many-body character of the collision. Along this line the classical trajectory Monte Carlo (CTMC) method has been quite successful in dealing with the atomic processes in ion-atom collisions. The CTMC method treats the atomic collisions in a non-perturbative manner. Classical equations of motions are solved numerically with randomly selected initial conditions. One of the main advantages of the CTMC method is that the three-body interaction is exactly taken into account during the collision at a classical level. The Integrated Tokamak Modeling Task Force (ITM-TF) was set up in Europe, under EFDA (the European Fusion Development Agreement), in 2004. The main target is to coordinate the European modeling effort and provide a complete European integrated modeling structure, with the highest degree of flexibility, for prediction of actual experiments and -in the long term - of the International Thermonuclear Experimental Reactor (ITER). For the ITM-TF a large numbers of atomic data is necessary. We have experience in theoretical interpretations of processes in ion-atom and ion-solid collisions based on classical descriptions. We developed a classical code for the determination of the cross sections for fundamental scattering processes. To demonstrate the validity and the efficiency of the present code we performed test calculations to determine the state selective charge exchange cross sections for collision system of N7+ + H(1s) at a 50 keV/amu projectile energy (see Fig. 1). The present results are in excellent agreement with the previous results. We also performed systematic calculation to determine charge exchange cross sections with other charge states of N ion. Further works are in progress and will be published soon. Acknowledgements. This work, supported by the European Communities under the contract of Association between EURATOM-HAS, was carried out within the framework of the Task Force on Integrated Tokamak Modelling of the European Fusion Development Agreement.