TERS v2.0: An improved version of TERS
We present a new version of the semimicroscopic Monte Carlo code “TERS”. The procedure for calculating multiple small angle Coulomb scattering of the residues in the target has been modified. Target-backing and residue charge-reset foils, which are often used in heavy ion-induced complete fusion rea...
Gespeichert in:
| Veröffentlicht in: | Computer physics communications 2009-11, Vol.180 (11), p.2392-2393 |
|---|---|
| 1. Verfasser: | |
| Format: | Artikel |
| Sprache: | eng |
| Schlagworte: | |
| Online-Zugang: | Volltext |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| container_end_page | 2393 |
|---|---|
| container_issue | 11 |
| container_start_page | 2392 |
| container_title | Computer physics communications |
| container_volume | 180 |
| creator | Nath, S. |
| description | We present a new version of the semimicroscopic Monte Carlo code “TERS”. The procedure for calculating multiple small angle Coulomb scattering of the residues in the target has been modified. Target-backing and residue charge-reset foils, which are often used in heavy ion-induced complete fusion reactions, are included in the code.
Program title: TERS v2.0
Catalogue identifier: AEBD_v2_0
Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBD_v2_0.html
Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland
Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html
No. of lines in distributed program, including test data, etc.: 7309
No. of bytes in distributed program, including test data, etc.: 1 219 555
Distribution format: tar.gz
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
Classification: 17.7
External routines: pgplot graphics subroutine library [1] should be installed in the system for generating residue trajectory plots. (The library is included in the CPC distribution file.)
Catalogue identifier of previous version: AEBD_v1_0
Journal reference of previous version: Comput. Phys. Comm. 179 (2008) 492
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.
Solution method: 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 (〈ϑ〉2), used earlier [2], was found to be inadequate particularly for low energy heavy residues. Energy loss of be |
| doi_str_mv | 10.1016/j.cpc.2009.06.025 |
| format | Article |
| fullrecord | <record><control><sourceid>proquest_cross</sourceid><recordid>TN_cdi_proquest_miscellaneous_1770330017</recordid><sourceformat>XML</sourceformat><sourcesystem>PC</sourcesystem><els_id>S0010465509002069</els_id><sourcerecordid>1770330017</sourcerecordid><originalsourceid>FETCH-LOGICAL-c330t-924c135807a7b3e80e3b72bd6a2c2acd106a447214db69185c6b39dc069a4b413</originalsourceid><addsrcrecordid>eNp9kE1LAzEQhoMoWKs_wNuexMuuk49NNnoqpX5AQdB6DtlkCintbk3aBf-9KevZ01ze552Zh5BbChUFKh82ldu7igHoCmQFrD4jE9ooXTItxDmZAFAohazrS3KV0gYAlNJ8Qu5Wi4_PYmAVPBazrgi7fewH9MWAMYW-K_p1cUpck4u13Sa8-ZtT8vW8WM1fy-X7y9t8tiwd53AoNROO8roBZVXLsQHkrWKtl5Y5Zp2nIK0QilHhW6lpUzvZcu0dSG1FKyifkvuxN5_xfcR0MLuQHG63tsP-mAxVCvImoCpH6Rh1sU8p4trsY9jZ-GMomJMTszHZiTk5MSBNdpKZp5HB_MMQMJrkAnYOfYjoDsb34R_6F1msZZ0</addsrcrecordid><sourcetype>Aggregation Database</sourcetype><iscdi>true</iscdi><recordtype>article</recordtype><pqid>1770330017</pqid></control><display><type>article</type><title>TERS v2.0: An improved version of TERS</title><source>ScienceDirect Journals (5 years ago - present)</source><creator>Nath, S.</creator><creatorcontrib>Nath, S.</creatorcontrib><description>We present a new version of the semimicroscopic Monte Carlo code “TERS”. The procedure for calculating multiple small angle Coulomb scattering of the residues in the target has been modified. Target-backing and residue charge-reset foils, which are often used in heavy ion-induced complete fusion reactions, are included in the code.
Program title: TERS v2.0
Catalogue identifier: AEBD_v2_0
Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBD_v2_0.html
Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland
Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html
No. of lines in distributed program, including test data, etc.: 7309
No. of bytes in distributed program, including test data, etc.: 1 219 555
Distribution format: tar.gz
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
Classification: 17.7
External routines: pgplot graphics subroutine library [1] should be installed in the system for generating residue trajectory plots. (The library is included in the CPC distribution file.)
Catalogue identifier of previous version: AEBD_v1_0
Journal reference of previous version: Comput. Phys. Comm. 179 (2008) 492
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.
Solution method: 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 (〈ϑ〉2), used earlier [2], 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:1.A new method [3] for calculating multiple small angle Coulomb scattering of residues in the target has been adopted. The change is incorporated in function Weibull() in the program ters_pti2.c.2.Isotopically enriched targets are made on a thin backing foil (usually made of carbon) quite often. Energy loss of beam in the backing foil (assuming beam is made to pass through the backing foil first, which is the usual practice) need to be taken into account. This calls for minor changes in the input file ters_pti2.inp. Following is the modified list of input parameters in this file with explanation.Zp, Ap, Zt, At –Atomic no. and mass no. of projectile and target.Ep –Projectile energy [MeV] in laboratory.BeamSpot, TarThick –Dia. [mm] of (circular) beam spot and target thickness [mg/cm2].Zback, Aback, BackThick –Atomic no. and mass no. of the backing material and thickness [mg/cm2] of the backing foil.Qvalue, ILPM –Q value [MeV] for CN formation and inverse level density parameter.AlphaNo, ProtonNo, NeutronNo –Numbers of evaporated alphas, protons and neutrons.Salpha[c] –Alpha separation energies [MeV], to be left blank if no alpha evaporation.Sproton[c] –Proton separation energies [MeV], to be left blank if no proton evaporation.Sneutron[c] –Neutron separation energies [MeV], to be left blank if no neutron evaporation.NEVENT –Number of events i.e. residues to be considered by the program (maximum 5×105). The parameters in input line number 4 are new in this version. If the target is backed by a carbon foil of thickness 125 μg/cm2, the input line would look like “6 12 0.125”. If the target is self-supporting, i.e. there is no backing, value of the last parameter (thickness) should be zero. However, the first two parameters must not be left blank or have 0 values. The input line would look like “6 12 0.0” in this case.3.A new function ThinFoil() has been introduced in the program ters_tra2.c. A thin foil can be inserted anywhere along the path of the residues by calling this function using the following syntax: Status = ThinFoil(argument list); if (Status == 0) continue; The function is particularly useful to place a residue charge-reset foil (usually made of carbon) after the target and is described in Table 1.Table 1Description of the function ThinFoil() used in the program ters_tra2.c. Function name is case sensitive.Name of the functionJob of the functionList of argumentsDescription of argumentsThinFoil()To calculate ion energy loss in a thin foilint arg1, int arg2, int arg3arg1=atomic number of the foil material, arg2=mass number of the foil material, arg3=thickness of the foil [mg/cm2]4.There is a minor change in the input file ters_tra2.inp. Following is the modified list of input parameters in this file with explanation.Z0, A0 –Residue atomic number and mass number.E0, A00, q0 –Energy [MeV], mass no. and charge state of the reference particle.NEVENT –Number of events i.e. trajectories to be calculated.5.Program/input files which have been modified in this version are suffixed by “2” in their names (before the extension), e.g., ters_pti.c has been renamed ters_pti2.c. The complete list of files included in the distributed code can be viewed in the readme file.
Restrictions: 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 [6] 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.
Running time: From few seconds to several minutes depending on the reaction, number of events and separator layout.
References:[1] http://www.astro.caltech.edu/~tjp/pgplot/.[2] G.R. Lynch, O.I. Dahl, Nucl. Instr. Methods B 58 (1991) 6.[3] L. Meyer, Phys. Status Solidi 44 (1971) 253.[4] A.K. Sinha, N. Madhavan, J.J. Das, P. Sugathan, D.O. Kataria, A.P. Patro, G.K. Mehta, Nucl. Instr. Methods A 339 (1994) 543.[5] S. Nath, Nucl. Instr. Methods A 576 (2007) 403.[6] J.F. Ziegler, J.P. Biersack, U. Littmark, The Stopping and Range of Ions in Solids, vol. I, Pergamon Press, Oxford, 1984.</description><identifier>ISSN: 0010-4655</identifier><identifier>EISSN: 1879-2944</identifier><identifier>EISSN: 1386-9485</identifier><identifier>DOI: 10.1016/j.cpc.2009.06.025</identifier><language>eng</language><publisher>Elsevier B.V</publisher><subject>Backing ; Beams (structural) ; Complete fusion reaction ; Evaporation ; Foils ; Ion optics ; Mathematical analysis ; Monte Carlo simulation ; Recoil separator ; Residues ; Separators ; Trajectories ; Transmission efficiency</subject><ispartof>Computer physics communications, 2009-11, Vol.180 (11), p.2392-2393</ispartof><rights>2009 Elsevier B.V.</rights><lds50>peer_reviewed</lds50><woscitedreferencessubscribed>false</woscitedreferencessubscribed><citedby>FETCH-LOGICAL-c330t-924c135807a7b3e80e3b72bd6a2c2acd106a447214db69185c6b39dc069a4b413</citedby></display><links><openurl>$$Topenurl_article</openurl><openurlfulltext>$$Topenurlfull_article</openurlfulltext><thumbnail>$$Tsyndetics_thumb_exl</thumbnail><linktohtml>$$Uhttps://dx.doi.org/10.1016/j.cpc.2009.06.025$$EHTML$$P50$$Gelsevier$$H</linktohtml><link.rule.ids>314,777,781,3537,27905,27906,45976</link.rule.ids></links><search><creatorcontrib>Nath, S.</creatorcontrib><title>TERS v2.0: An improved version of TERS</title><title>Computer physics communications</title><description>We present a new version of the semimicroscopic Monte Carlo code “TERS”. The procedure for calculating multiple small angle Coulomb scattering of the residues in the target has been modified. Target-backing and residue charge-reset foils, which are often used in heavy ion-induced complete fusion reactions, are included in the code.
Program title: TERS v2.0
Catalogue identifier: AEBD_v2_0
Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBD_v2_0.html
Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland
Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html
No. of lines in distributed program, including test data, etc.: 7309
No. of bytes in distributed program, including test data, etc.: 1 219 555
Distribution format: tar.gz
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
Classification: 17.7
External routines: pgplot graphics subroutine library [1] should be installed in the system for generating residue trajectory plots. (The library is included in the CPC distribution file.)
Catalogue identifier of previous version: AEBD_v1_0
Journal reference of previous version: Comput. Phys. Comm. 179 (2008) 492
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.
Solution method: 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 (〈ϑ〉2), used earlier [2], 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:1.A new method [3] for calculating multiple small angle Coulomb scattering of residues in the target has been adopted. The change is incorporated in function Weibull() in the program ters_pti2.c.2.Isotopically enriched targets are made on a thin backing foil (usually made of carbon) quite often. Energy loss of beam in the backing foil (assuming beam is made to pass through the backing foil first, which is the usual practice) need to be taken into account. This calls for minor changes in the input file ters_pti2.inp. Following is the modified list of input parameters in this file with explanation.Zp, Ap, Zt, At –Atomic no. and mass no. of projectile and target.Ep –Projectile energy [MeV] in laboratory.BeamSpot, TarThick –Dia. [mm] of (circular) beam spot and target thickness [mg/cm2].Zback, Aback, BackThick –Atomic no. and mass no. of the backing material and thickness [mg/cm2] of the backing foil.Qvalue, ILPM –Q value [MeV] for CN formation and inverse level density parameter.AlphaNo, ProtonNo, NeutronNo –Numbers of evaporated alphas, protons and neutrons.Salpha[c] –Alpha separation energies [MeV], to be left blank if no alpha evaporation.Sproton[c] –Proton separation energies [MeV], to be left blank if no proton evaporation.Sneutron[c] –Neutron separation energies [MeV], to be left blank if no neutron evaporation.NEVENT –Number of events i.e. residues to be considered by the program (maximum 5×105). The parameters in input line number 4 are new in this version. If the target is backed by a carbon foil of thickness 125 μg/cm2, the input line would look like “6 12 0.125”. If the target is self-supporting, i.e. there is no backing, value of the last parameter (thickness) should be zero. However, the first two parameters must not be left blank or have 0 values. The input line would look like “6 12 0.0” in this case.3.A new function ThinFoil() has been introduced in the program ters_tra2.c. A thin foil can be inserted anywhere along the path of the residues by calling this function using the following syntax: Status = ThinFoil(argument list); if (Status == 0) continue; The function is particularly useful to place a residue charge-reset foil (usually made of carbon) after the target and is described in Table 1.Table 1Description of the function ThinFoil() used in the program ters_tra2.c. Function name is case sensitive.Name of the functionJob of the functionList of argumentsDescription of argumentsThinFoil()To calculate ion energy loss in a thin foilint arg1, int arg2, int arg3arg1=atomic number of the foil material, arg2=mass number of the foil material, arg3=thickness of the foil [mg/cm2]4.There is a minor change in the input file ters_tra2.inp. Following is the modified list of input parameters in this file with explanation.Z0, A0 –Residue atomic number and mass number.E0, A00, q0 –Energy [MeV], mass no. and charge state of the reference particle.NEVENT –Number of events i.e. trajectories to be calculated.5.Program/input files which have been modified in this version are suffixed by “2” in their names (before the extension), e.g., ters_pti.c has been renamed ters_pti2.c. The complete list of files included in the distributed code can be viewed in the readme file.
Restrictions: 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 [6] 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.
Running time: From few seconds to several minutes depending on the reaction, number of events and separator layout.
References:[1] http://www.astro.caltech.edu/~tjp/pgplot/.[2] G.R. Lynch, O.I. Dahl, Nucl. Instr. Methods B 58 (1991) 6.[3] L. Meyer, Phys. Status Solidi 44 (1971) 253.[4] A.K. Sinha, N. Madhavan, J.J. Das, P. Sugathan, D.O. Kataria, A.P. Patro, G.K. Mehta, Nucl. Instr. Methods A 339 (1994) 543.[5] S. Nath, Nucl. Instr. Methods A 576 (2007) 403.[6] J.F. Ziegler, J.P. Biersack, U. Littmark, The Stopping and Range of Ions in Solids, vol. I, Pergamon Press, Oxford, 1984.</description><subject>Backing</subject><subject>Beams (structural)</subject><subject>Complete fusion reaction</subject><subject>Evaporation</subject><subject>Foils</subject><subject>Ion optics</subject><subject>Mathematical analysis</subject><subject>Monte Carlo simulation</subject><subject>Recoil separator</subject><subject>Residues</subject><subject>Separators</subject><subject>Trajectories</subject><subject>Transmission efficiency</subject><issn>0010-4655</issn><issn>1879-2944</issn><issn>1386-9485</issn><fulltext>true</fulltext><rsrctype>article</rsrctype><creationdate>2009</creationdate><recordtype>article</recordtype><recordid>eNp9kE1LAzEQhoMoWKs_wNuexMuuk49NNnoqpX5AQdB6DtlkCintbk3aBf-9KevZ01ze552Zh5BbChUFKh82ldu7igHoCmQFrD4jE9ooXTItxDmZAFAohazrS3KV0gYAlNJ8Qu5Wi4_PYmAVPBazrgi7fewH9MWAMYW-K_p1cUpck4u13Sa8-ZtT8vW8WM1fy-X7y9t8tiwd53AoNROO8roBZVXLsQHkrWKtl5Y5Zp2nIK0QilHhW6lpUzvZcu0dSG1FKyifkvuxN5_xfcR0MLuQHG63tsP-mAxVCvImoCpH6Rh1sU8p4trsY9jZ-GMomJMTszHZiTk5MSBNdpKZp5HB_MMQMJrkAnYOfYjoDsb34R_6F1msZZ0</recordid><startdate>200911</startdate><enddate>200911</enddate><creator>Nath, S.</creator><general>Elsevier B.V</general><scope>AAYXX</scope><scope>CITATION</scope><scope>7SC</scope><scope>7U5</scope><scope>8FD</scope><scope>H8D</scope><scope>JQ2</scope><scope>L7M</scope><scope>L~C</scope><scope>L~D</scope></search><sort><creationdate>200911</creationdate><title>TERS v2.0: An improved version of TERS</title><author>Nath, S.</author></sort><facets><frbrtype>5</frbrtype><frbrgroupid>cdi_FETCH-LOGICAL-c330t-924c135807a7b3e80e3b72bd6a2c2acd106a447214db69185c6b39dc069a4b413</frbrgroupid><rsrctype>articles</rsrctype><prefilter>articles</prefilter><language>eng</language><creationdate>2009</creationdate><topic>Backing</topic><topic>Beams (structural)</topic><topic>Complete fusion reaction</topic><topic>Evaporation</topic><topic>Foils</topic><topic>Ion optics</topic><topic>Mathematical analysis</topic><topic>Monte Carlo simulation</topic><topic>Recoil separator</topic><topic>Residues</topic><topic>Separators</topic><topic>Trajectories</topic><topic>Transmission efficiency</topic><toplevel>peer_reviewed</toplevel><toplevel>online_resources</toplevel><creatorcontrib>Nath, S.</creatorcontrib><collection>CrossRef</collection><collection>Computer and Information Systems Abstracts</collection><collection>Solid State and Superconductivity Abstracts</collection><collection>Technology Research Database</collection><collection>Aerospace Database</collection><collection>ProQuest Computer Science Collection</collection><collection>Advanced Technologies Database with Aerospace</collection><collection>Computer and Information Systems Abstracts Academic</collection><collection>Computer and Information Systems Abstracts Professional</collection><jtitle>Computer physics communications</jtitle></facets><delivery><delcategory>Remote Search Resource</delcategory><fulltext>fulltext</fulltext></delivery><addata><au>Nath, S.</au><format>journal</format><genre>article</genre><ristype>JOUR</ristype><atitle>TERS v2.0: An improved version of TERS</atitle><jtitle>Computer physics communications</jtitle><date>2009-11</date><risdate>2009</risdate><volume>180</volume><issue>11</issue><spage>2392</spage><epage>2393</epage><pages>2392-2393</pages><issn>0010-4655</issn><eissn>1879-2944</eissn><eissn>1386-9485</eissn><abstract>We present a new version of the semimicroscopic Monte Carlo code “TERS”. The procedure for calculating multiple small angle Coulomb scattering of the residues in the target has been modified. Target-backing and residue charge-reset foils, which are often used in heavy ion-induced complete fusion reactions, are included in the code.
Program title: TERS v2.0
Catalogue identifier: AEBD_v2_0
Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEBD_v2_0.html
Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland
Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html
No. of lines in distributed program, including test data, etc.: 7309
No. of bytes in distributed program, including test data, etc.: 1 219 555
Distribution format: tar.gz
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
Classification: 17.7
External routines: pgplot graphics subroutine library [1] should be installed in the system for generating residue trajectory plots. (The library is included in the CPC distribution file.)
Catalogue identifier of previous version: AEBD_v1_0
Journal reference of previous version: Comput. Phys. Comm. 179 (2008) 492
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.
Solution method: 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 (〈ϑ〉2), used earlier [2], 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:1.A new method [3] for calculating multiple small angle Coulomb scattering of residues in the target has been adopted. The change is incorporated in function Weibull() in the program ters_pti2.c.2.Isotopically enriched targets are made on a thin backing foil (usually made of carbon) quite often. Energy loss of beam in the backing foil (assuming beam is made to pass through the backing foil first, which is the usual practice) need to be taken into account. This calls for minor changes in the input file ters_pti2.inp. Following is the modified list of input parameters in this file with explanation.Zp, Ap, Zt, At –Atomic no. and mass no. of projectile and target.Ep –Projectile energy [MeV] in laboratory.BeamSpot, TarThick –Dia. [mm] of (circular) beam spot and target thickness [mg/cm2].Zback, Aback, BackThick –Atomic no. and mass no. of the backing material and thickness [mg/cm2] of the backing foil.Qvalue, ILPM –Q value [MeV] for CN formation and inverse level density parameter.AlphaNo, ProtonNo, NeutronNo –Numbers of evaporated alphas, protons and neutrons.Salpha[c] –Alpha separation energies [MeV], to be left blank if no alpha evaporation.Sproton[c] –Proton separation energies [MeV], to be left blank if no proton evaporation.Sneutron[c] –Neutron separation energies [MeV], to be left blank if no neutron evaporation.NEVENT –Number of events i.e. residues to be considered by the program (maximum 5×105). The parameters in input line number 4 are new in this version. If the target is backed by a carbon foil of thickness 125 μg/cm2, the input line would look like “6 12 0.125”. If the target is self-supporting, i.e. there is no backing, value of the last parameter (thickness) should be zero. However, the first two parameters must not be left blank or have 0 values. The input line would look like “6 12 0.0” in this case.3.A new function ThinFoil() has been introduced in the program ters_tra2.c. A thin foil can be inserted anywhere along the path of the residues by calling this function using the following syntax: Status = ThinFoil(argument list); if (Status == 0) continue; The function is particularly useful to place a residue charge-reset foil (usually made of carbon) after the target and is described in Table 1.Table 1Description of the function ThinFoil() used in the program ters_tra2.c. Function name is case sensitive.Name of the functionJob of the functionList of argumentsDescription of argumentsThinFoil()To calculate ion energy loss in a thin foilint arg1, int arg2, int arg3arg1=atomic number of the foil material, arg2=mass number of the foil material, arg3=thickness of the foil [mg/cm2]4.There is a minor change in the input file ters_tra2.inp. Following is the modified list of input parameters in this file with explanation.Z0, A0 –Residue atomic number and mass number.E0, A00, q0 –Energy [MeV], mass no. and charge state of the reference particle.NEVENT –Number of events i.e. trajectories to be calculated.5.Program/input files which have been modified in this version are suffixed by “2” in their names (before the extension), e.g., ters_pti.c has been renamed ters_pti2.c. The complete list of files included in the distributed code can be viewed in the readme file.
Restrictions: 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 [6] 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.
Running time: From few seconds to several minutes depending on the reaction, number of events and separator layout.
References:[1] http://www.astro.caltech.edu/~tjp/pgplot/.[2] G.R. Lynch, O.I. Dahl, Nucl. Instr. Methods B 58 (1991) 6.[3] L. Meyer, Phys. Status Solidi 44 (1971) 253.[4] A.K. Sinha, N. Madhavan, J.J. Das, P. Sugathan, D.O. Kataria, A.P. Patro, G.K. Mehta, Nucl. Instr. Methods A 339 (1994) 543.[5] S. Nath, Nucl. Instr. Methods A 576 (2007) 403.[6] J.F. Ziegler, J.P. Biersack, U. Littmark, The Stopping and Range of Ions in Solids, vol. I, Pergamon Press, Oxford, 1984.</abstract><pub>Elsevier B.V</pub><doi>10.1016/j.cpc.2009.06.025</doi><tpages>2</tpages></addata></record> |
| fulltext | fulltext |
| identifier | ISSN: 0010-4655 |
| ispartof | Computer physics communications, 2009-11, Vol.180 (11), p.2392-2393 |
| issn | 0010-4655 1879-2944 1386-9485 |
| language | eng |
| recordid | cdi_proquest_miscellaneous_1770330017 |
| source | ScienceDirect Journals (5 years ago - present) |
| subjects | Backing Beams (structural) Complete fusion reaction Evaporation Foils Ion optics Mathematical analysis Monte Carlo simulation Recoil separator Residues Separators Trajectories Transmission efficiency |
| title | TERS v2.0: An improved version of TERS |
| url | https://sfx.bib-bvb.de/sfx_tum?ctx_ver=Z39.88-2004&ctx_enc=info:ofi/enc:UTF-8&ctx_tim=2025-01-19T17%3A55%3A28IST&url_ver=Z39.88-2004&url_ctx_fmt=infofi/fmt:kev:mtx:ctx&rfr_id=info:sid/primo.exlibrisgroup.com:primo3-Article-proquest_cross&rft_val_fmt=info:ofi/fmt:kev:mtx:journal&rft.genre=article&rft.atitle=TERS%20v2.0:%20An%20improved%20version%20of%20TERS&rft.jtitle=Computer%20physics%20communications&rft.au=Nath,%20S.&rft.date=2009-11&rft.volume=180&rft.issue=11&rft.spage=2392&rft.epage=2393&rft.pages=2392-2393&rft.issn=0010-4655&rft.eissn=1879-2944&rft_id=info:doi/10.1016/j.cpc.2009.06.025&rft_dat=%3Cproquest_cross%3E1770330017%3C/proquest_cross%3E%3Curl%3E%3C/url%3E&disable_directlink=true&sfx.directlink=off&sfx.report_link=0&rft_id=info:oai/&rft_pqid=1770330017&rft_id=info:pmid/&rft_els_id=S0010465509002069&rfr_iscdi=true |