|Special Applications of TRIM
Calculating Ion Energy Reducers (Energy Degraders)
TRIM may be used to calculate the energy loss through foils or thick blocks of materials used to lower the energy of a beam of particles. This is often done for light ions, such as protons, to obtain lower energies quickly. For example, if one starts with a beam of monoenergetic protons at 158.6 MeV and introduces a block of plexiglas (called perspex in Europe) 14.1 cm thick, the exiting beam of particles will be at 30 MeV. Degrader blocks made with light atoms can reduce particle energy with a minimum of straggle and little long term radioactivity (plexiglas is C6H4O2).
Two problems occur in planning "energy degraders". First is the calculation of the proper block thickness, and second is the calculation of the final beam energy spread. Below we illustrate both calculations.
The calculation of the required degrader thickness consists of four basic steps. We illustrate these steps below using the example of protons being slowed by a block of plexiglas from an initial energy of 158.6 MeV to a final energy of 30 MeV.
1 Make a table of proton stopping in plexiglas by starting SRIM and going to Stopping and Range Tables.
The default ion is H so this can remain.
Change the Highest Ion Energy to 160000 keV.
Press Compound Dictionary and open the file Common Target Materials. In this list is “Plexiglas / Acrylic”. Select this material and press Add to Target. A pop-up window will ask if you wish to keep the Compound Correction to the stopping of H ions in Plexiglas. Answer Yes.
At the bottom left of the setup window is Stopping Power Units. Open the drop-down menu as specify MeV/mm. This will tell you the energy loss of the beam in convenient units.
Press Calculate Table.
2 Using the table generated in (1), find the range of protons for the initial beam energy and estimate the depth at which the protons will be at 30 MeV by looking at the ranges and working backwards. For example, the table will show that 160 MeV protons have a range of 155 mm in plexiglas, and 30 MeV protons will go 7.7 mm. Hence your first guess for the thickness of plexiglas to reduce 160 MeV protons to 30 MeV will be 155 - 7.7 = 147.3 mm thick. Further, since the initial proton beam is not 160 MeV but 158.6 MeV, one interpolates the distance protons travel between 160 MeV to 158.6 MeV as 2.3 mm, so the initial block size for an incident 158.6 MeV beam will be 147.3 - 2.3 = 145 mm.
3 Close the Table window.
Press the Main Menu button.
Press TRIM to go to the TRIM setup window.
Specify ions of H (158.6 MeV) into a target of plexiglas (145 mm).
At the bottom middle of the Setup window is Output Disk Files. Check the box beside Transmitted Ions.
Double check your inputs, and then press Save Input and Run TRIM.
Calculate about 200 ions, then stop. Close TRIM (do not Save). Quit the input window.
3 You now need to look at the file: SRIM Outputs/TRANSMIT.txt. This will contain the transmitted protons through 145 mm of Plexiglas. You need to find the mean energy of these particles, and their average deviation from this mean (straggling). This can be done in Excel or other statistical packages. There is a BASIC program included with SRIM (/SRIM Tools/TRANSMIT.BAS) if you can run this type of program. It will read the file TRANSMIT.txt and give you a full statistical analysis.
4 Use the difference between your desired final energy (30 MeV) and mean final energy found in step 3 above, along with the stopping power tables, made in step 1 above, to estimate a new target thickness. For example, if the protons exit at 32 MeV and the stopping power for protons at 30 MeV is 2.16 MeV/mm, then you need (2 MeV)/(2.16 MeV/mm) = 0.92 mm more plexiglas in the block to bring the beam down to 30 MeV. You then use TRIM to make up a new plexiglas target with a width of 145 mm (initial estimate) + 0.92 mm (correction) = 145.92 mm. Loop through steps 2-4 until you come get as close as you want to the desired thickness.
The final energy straggle of the beam after it has passed through the energy reducing block is a mixture of two components: (1) the straggle introduced by random collisions in the degrader block and (2) the broadening of any incident beam energy straggle by the block (straggle is the spread of the ion beam energy spectrum). We calculate the straggle of the incident monoenergetic beam with the TRIM calculations above, steps 1-4, and we find, for example, a straggle of about 2.0 MeV when we degrade protons from 160 MeV to 30 MeV using plexiglas (this comes out of the analysis program TRANSMIT.BAS) To this broadening, we have to add the straggle due to expansion of the initial beam energy spread, which may be the dominant factor in the final energy width. Assuming that the incident beam is not monoenergetic, but has a spread of 1 MeV at 160 MeV, we will find that this spread expands to 3.6 MeV by the time the beam degrades to 30 MeV. This expansion occurs because of the non-linearity of stopping powers with particle energy as explained below.
Assume that the energy straggle of the incident beam is . We show below that the final energy straggle of the exiting beam, ', can be estimated as :
' = (S'/S) ,
where S and S' are the stopping powers, dE/dx, in the degrader material at the initial and final beam energy.
Consider two particles of the incident beam with energies E1 and E2, separated by the incident beam energy straggle, . The final energy straggle of these two particles, after transiting the energy-lowering block of material, is defined as ' = E1' – E2'. where E1' and E2' are the final energies of particles 1 and 2.
Considering the incident particles, assuming E1 > E2, , then particle #1 will travel a short distance into the degrader block before its energy will be reduced to that of particle #2. This distance, x, is about x = / S, where S is the particle stopping power (dE/dx) at energy E. After particle #1 reaches depth x, then both particles #1 and #2 will transit through the block similarly until particle #1 reaches the end of the block. Particle #2 will still have more of the block to traverse, an amount x' = ' / S' , where ' = E1' - E2', the final energy difference between the two particles, and S' is the stopping power at the final beam energy, E'. But the block length is the same for both particles, so x = x', and therefore / S = ' / S' . This can be solved for the energy straggle, ' = (S'/S) .
Note that several approximations are made in this argument, especially that S(E1) = S(E2) and S(E1') = S(E2'), and that the total energy loss is much larger than the straggling, (E - E') >> .
For an example, protons at 160 MeV in plexiglas have a stopping power of 0.593 MeV/mm and at 30 MeV have a stopping power of 2.157 MeV/mm. Assuming the initial beam straggle, , is 1 MeV, then the final beam energy straggle, ' = (1 MeV) (2.157 MeV/mm) / (.593 MeV/mm) = 3.6 MeV.
The total beam spread through the energy-degrader block is the rms total value of the two separate quantities : (1) the beam straggle through the block, 2.0 MeV, and (2) the beam energy broadening due to non-linear stopping powers, 3.6 MeV, which, when added in quadrature, gives a total straggle of 4.1 MeV .when protons at 160 MeV are degraded to 30 MeV.
The above comments come from Dr. Bernard Gottschalk, Harvard Cyclotron Laboratory, Harvard University, USA.
SRIM Special – Degrader Blocks.doc Page of 6/26/2016