Calculation of some important parameters in the muon catalyzed fusion by Monte-Carlo simulation

It is determined, how we can take into account the contribution of the nonthermalyzed tµ atoms [1] in the kinetics of muon catalyzed fusion in the steady-state condition. For this purpose we should determine the probability of the dtµ molecular formation by the epithermal tµ atoms during thermal-ization processes (η parameter). We determine the η parameter by the Monte-Carlo simulation [2]. However, to calculate the important parameters in µCF such as fusion yield per muon and cycling rate, solving the kinetics of µCF in the steady-state condition [2, 3] is not completely accurate method. In the present work, the µCF cycle is totally simulated by the Monte-Carlo method. We can approach to the actual situation in the µCF cycle by this method. We have prepared a computer code for this purpose. The simulation starts since a muon enters into the D/T mixture. The collision partner, kinetic energy and spin state of the muonic atom are determined, then we follow the trajectories of the muonic atoms in the collisional processes by using their cross-sections [4]. During the collisions of the tµ atoms, the molecular formation may also take place. We have calculated the fusion yield per muon in the different concentration of hydrogen isotopes. These results are compared with the results due to solve the kinetic equations of µCF cycle[2, 3]. In addition, we can calculate some other useful quantities which increase our knowledge about the processes in the µCF cycle. For example