Programs in Physics & Physical Chemistry
|[Licence| Download | New Version Template] aebd_v2_0.tar.gz(1190 Kbytes)|
|Manuscript Title: TERS v2.0: An improved version of TERS|
|Authors: S. Nath|
|Program title: TERS v2.0|
|Catalogue identifier: AEBD_v2_0|
Distribution format: tar.gz
|Journal reference: Comput. Phys. Commun. 180(2009)2392|
|Programming language: C.|
|Computer: The code has been developed and tested on a PC with Intel Pentium IV processor.|
|Operating system: Linux.|
|RAM: About 8 Mbytes|
|Keywords: Monte Carlo simulation, Recoil separator, Complete fusion reaction, Ion optics, Transmission efficiency.|
|PACS: 02.70.Uu, 29.30.h, 25.70.Jj, 41.85.p..|
External routines: pgplot graphics subroutine library  should be installed in the system for generating residue trajectory plots. (Code included in the CPC distribution file).
Does the new version supersede the previous version?: Yes
Nature of problem:
Recoil separators are employed to select and identify nuclei of interest, produced in a nuclear reaction, rejecting unreacted beam and other undesired reaction products. It is important to know what fraction of the selected nuclei, leaving the target, reach the detection system. This information is crucial for determining absolute cross section of the studied reaction.
Interaction of projectiles with target nuclei is treated event by event, semimicro- scopically. Position and angle (with respect to beam direction), energy and charge state of the reaction products are calculated by Monte Carlo method. Trajectory of each nuclei inside the separator is then calculated by ion optical transfer matrix method. Ratio of the number of trajectories completing their journey up to the detection system to the total number of trajectories is a direct measure of absolute transmission efficiency of the separator.
Reasons for new version:
The method for calculating mean squared scattering angle (<v>2), used earlier, was found to be inadequate particularly for low energy heavy residues. Energy loss of beam in the target-backing foil and energy loss of residues in the charge-reset foil (wherever used) needed to be taken into account for better matching of simulated residue parameters with measurements.
Summary of revisions:
The present version of the code is applicable to complete fusion reactions only. Calculation of transmission efficiency has been illustrated with a specific recoil separator, viz. the Heavy Ion Reaction Analyzer (HIRA)[4,5], at IUAC. One has to make necessary changes in the code, while performing calculations for other recoil separators. Also, atomic number of the residual nucleus should not exceed 92, as the method used for calculating stopping power of ions  is valid for Z ≤ 92. The code can perform energy loss calculation only in elemental targets and foils (i.e. compounds or alloys are not supported). Further, number of events (NEVENT) in ters_tra2.inp should not exceed the same in ters_pti2.inp.
From a few seconds to several minutes depending on the reaction, number of events and separator layout.
|||G.R. Lynch, O.I. Dahl, Nucl. Instr. Methods B 58 (1991) 6.|
|||L. Meyer, Phys. Status Solidi 44 (1971) 253.|
|||A.K. Sinha, N. Madhavan, J.J. Das, P. Sugathan, D.O. Kataria, A.P. Patro, G.K. Mehta, Nucl. Instr. and Methods A 339 (1994) 543.|
|||S. Nath, Nucl. Instr. Methods A 576 (2007) 403.|
|||J.F. Ziegler, J.P. Biersack, U. Littmark, The stopping and range of ions in solids, Vol. I (Pergamon Press, Oxford, 1984).|
|Disclaimer | ScienceDirect | CPC Journal | CPC | QUB|